BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	887,393	T→A	7.6%	intergenic (‑54/‑30)	ybjJ ← / → rcdA	putative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	887,393	0	T	A	7.6%	614.9 / 27.7	181	intergenic (‑54/‑30)	ybjJ/rcdA	putative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator
Reads supporting (aligned to +/- strand): ref base T (92/75); new base A (6/8); total (98/83)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 197 total aligned reads displayed.

GTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTC  >  NC_000913/887123‑887666

ttttCCCTAACCACCCCTAACCTAGAATACCGCTCATCACACCGCTAGACAGACCAAAGAGAACCCCGGCCATGTCAGCGCTCGAGACAGAGAGGATATCGCGGAGAGCGGGCGTACAGGTGGCCCAGGACGCCATTAACAGGCAGAGAAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGGTGTGTACATTTGTa                                                                                                                                                                                                                                                                                <  4:441651/274‑1 (MQ=255)

 tttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGTTGTGTACATTTGTac                                                                                                                                                                                                                                                                               >  2:114045/1‑275 (MQ=255)

          ccACACCGGGCGAAGAATGCCGCTCATCGNNCCGATCGACGGACCAAAGAGAACACCGCCCNTTTCAGCGATCTAGACAGGGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGCCGCCATTAACAGGCCTNGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCa                                                                                                                                                                                                                                                                      <  3:100679/275‑1 (MQ=255)

              gCCGAGCAGCGAATACCGCTCATCGAACCGATTGACACATTAACGAGAACACCGCCCACTTTAGCGACCGCGAGCGAGAGAATATCGCGGATAGCAGGCGAATGGGTCCCCCAGGCCGTCATTAACACGCCCGGCACAAAGAAGAACATTAACAGCTCACAGGTCTGGAGTTTCAATGCATTACGGCACGGATTTTCGTGCATAGATTACGTCAAAATAAGAAGTGGGCAGACAACATTAGACCGGCTGTTGACCTTTGTAAAAAATTGcgagaa                                                                                                                                                                                                                                                                  >  4:314046/1‑270 (MQ=255)

                              cGCCCATCGAACCGATCGACAGACCAAAGAGAACACAGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAAGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACAAAATTGCAGAGAGAGGAAATGACATGCg                                                                                                                                                                                                                                                  <  2:517233/275‑1 (MQ=255)

                                      gAACCGATCGACAGACCAAAGGGCACACCGCCCCTTTCAGCGATCGAGACAGAGAAAATATCGCGGATAACAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCACAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                                                                                                                                                                                                               <  4:114044/270‑1 (MQ=255)

                                         ccGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAcaca                                                                                                                                                                                                                                                                            <  4:88216/238‑1 (MQ=255)

                                                         cgataaCACCGCCCATTTCAGCGATCGAGACAAAGAGAATATTGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACATGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                                                                                                                                                       <  2:327324/271‑1 (MQ=255)

                                                                  cGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGCGAAAAAATTAt                                                                                                                                                                                                              >  1:136533/1‑275 (MQ=255)

                                                                                                            cAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTGCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGCGAGAGGAAATGAAATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAACTATCCAGGCAACGCGGGAGGCGGGGAAACTTTCCAGAATACAtcc                                                                                                                                                                    >  4:154511/1‑273 (MQ=255)

                                                                                                              ggCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                                                                                                                                                                                                    >  4:182377/1‑161 (MQ=255)

                                                                                                                       gTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGACGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGGGAGAGTAAATGACATGCGTCGCGCTACCTATCAGCAACGGCGAGAAAAAATTACCCAGTCCACGCTGGAGGCGGTGAAACTTACCGAAATACAATCTGTTTCGcct                                                                                                                                                          >  4:537637/1‑272 (MQ=255)

                                                                                                                            ccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                                                                                                                                                  >  8:715116/1‑149 (MQ=255)

                                                                                                                                gACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTac                                                                                                                                                                                                                                                                               <  5:697413/148‑1 (MQ=255)

                                                                                                                                   gCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                            <  6:69864/148‑1 (MQ=255)

                                                                                                                                       ttAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                            <  6:893156/144‑1 (MQ=255)

                                                                                                                                       ttAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                                                                                                                                                                                                        >  6:872461/1‑148 (MQ=255)

                                                                                                                                             ggCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                                                                                                                                                                                                 <  7:247973/149‑1 (MQ=255)

                                                                                                                                              gCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                                                                                                                                                                                                 <  6:537286/148‑1 (MQ=255)

                                                                                                                                               ccTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGACGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCCCGCCGCAGGCGGTGAACCTTTACGGAATACATTCTTTTA‑GCCCCACACAAATTGCTTCCCTTTCCgggg                                                                                                                                >  4:387550/1‑275 (MQ=255)

                                                                                                                                               ccTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                                                                                                                                                <  5:884214/148‑1 (MQ=255)

                                                                                                                                                cTGGCAAAAAGAAGAACATAAACAGAGCCCAGGTTCGGAGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGt                                                                                                                               <  2:298069/275‑1 (MQ=255)

                                                                                                                                                 tGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                                                                                                                                                                                                              >  6:589315/1‑148 (MQ=255)

                                                                                                                                                  ggCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                                <  2:592210/129‑1 (MQ=255)

                                                                                                                                                   gCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                                                                                                                                                                                                           <  7:715116/149‑1 (MQ=255)

                                                                                                                                                         agaagaACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGCAACTTTACGTAATACATGCTGTTACGCACCGCC‑AAATTGCTACCCTTGCTGGGGTACCGTTggg                                                                                                                      >  2:82774/1‑275 (MQ=255)

                                                                                                                                                           aagaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGc                                                                                                                                                                                                                                                   >  8:548820/1‑149 (MQ=255)

                                                                                                                                                             gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                                <  8:824150/118‑1 (MQ=255)

                                                                                                                                                             gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                                >  7:824150/1‑118 (MQ=255)

                                                                                                                                                              aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                                                                                                                                                                                                           <  5:652754/138‑1 (MQ=255)

                                                                                                                                                              aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa                                                                                                                                                                                     >  2:358345/1‑208 (MQ=255)

                                                                                                                                                              aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                            >  6:675472/1‑121 (MQ=255)

                                                                                                                                                               aCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                                                                                                                                                                                                               >  7:799705/1‑149 (MQ=255)

                                                                                                                                                                 aTAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                                                                                                                                                                                                             >  7:93522/1‑149 (MQ=255)

                                                                                                                                                                 aTAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgc                                                                                                                                                                                                                                              >  5:42812/1‑148 (MQ=255)

                                                                                                                                                                  tAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCgcnc                                                                                                                                                                                                                                            <  7:329186/149‑1 (MQ=255)

                                                                                                                                                                  tAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACttt                                                                                                                                                                                 <  4:358343/208‑1 (MQ=255)

                                                                                                                                                                     aCAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                                                                                                                                                <  7:181307/126‑1 (MQ=255)

                                                                                                                                                                      cAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                                                                                                                                                                                                         <  6:535591/148‑1 (MQ=255)

                                                                                                                                                                      cAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                                                                                                                                                                                                         <  5:872461/148‑1 (MQ=255)

                                                                                                                                                                        gcgcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACctacta                                                                                                       >  1:143600/1‑275 (MQ=255)

                                                                                                                                                                             cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                                                                                                                                                                                  <  5:765361/148‑1 (MQ=255)

                                                                                                                                                                                     gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                                                                                                                                                                          <  5:720959/148‑1 (MQ=255)

                                                                                                                                                                                          tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTGTTTTTCAGGAATtgatga                                                                                     >  2:379561/1‑275 (MQ=255)

                                                                                                                                                                                            aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                                                                                                                                                                                                  >  8:169637/1‑149 (MQ=255)

                                                                                                                                                                                            aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaa                                                                                                                                                                                                                   <  6:426099/148‑1 (MQ=255)

                                                                                                                                                                                             aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                                                                                                                                                                                                  >  6:631000/1‑148 (MQ=255)

                                                                                                                                                                                                  ttACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                                                                                                                                                                             >  6:467879/1‑148 (MQ=255)

                                                                                                                                                                                                   tACGTGAAGAATTTACGGCCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGCGAGGAAATGACATGCGTCGCGCTAGCGAACCGCAACGGCGAGAAAAAATTATCCAGGCCACGCGGGAGGCGGTGAAACTTTACGGAATACATGATGTTACGCACCGCA‑AAATTGATACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgattag                                                                            >  4:281081/1‑270 (MQ=255)

                                                                                                                                                                                                    aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                                                                                                                                                                                                          <  8:699899/149‑1 (MQ=255)

                                                                                                                                                                                                         aagaaNTTACNGTCATAGATANTGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCg                                                                                                                                   <  4:770969/214‑1 (MQ=255)

                                                                                                                                                                                                                tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCCAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                                                                                                                                                                             <  8:442883/134‑1 (MQ=255)

                                                                                                                                                                                                                tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                                                                                                                                                                             >  7:442883/1‑134 (MQ=255)

                                                                                                                                                                                                                    gTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCg                                                                                                                                                                                           <  5:589315/148‑1 (MQ=255)

                                                                                                                                                                                                                     tCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg                                                                                                                                                                                          <  5:356948/148‑1 (MQ=255)

                                                                                                                                                                                                                      cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                                                                                                                                                                                                    >  1:203577/1‑57 (MQ=255)

                                                                                                                                                                                                                      cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt                                                                                                                                                                                         >  5:578242/1‑148 (MQ=255)

                                                                                                                                                                                                                       aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa                                                                                                                                                                                       >  7:877285/1‑149 (MQ=255)

                                                                                                                                                                                                                        tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATcc                                                                                                                                                                                                            >  3:46292/1‑127 (MQ=255)

                                                                                                                                                                                                                        tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa                                                                                                                                                                                       >  5:908217/1‑148 (MQ=255)

                                                                                                                                                                                                                          gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg                                                                                                                                                                                          >  6:324229/1‑143 (MQ=255)

                                                                                                                                                                                                                           aTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACt                                                                                                                                                                                   >  7:596719/1‑149 (MQ=255)

                                                                                                                                                                                                                            tCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGc                                                                                                                                                                                                        <  1:46293/127‑1 (MQ=255)

                                                                                                                                                                                                                              aCGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACttt                                                                                                                                                                                 <  6:578242/148‑1 (MQ=255)

                                                                                                                                                                                                                               cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGATTTGTTACTGGAGGCGTTCAGCAGTTTTACTgaatac                                                >  4:217400/1‑271 (MQ=255)

                                                                                                                                                                                                                               cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTgagaac                                                >  2:541415/1‑273 (MQ=255)

                                                                                                                                                                                                                                gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                                                                                                                                                                              >  8:78717/1‑149 (MQ=255)

                                                                                                                                                                                                                                 tCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                                                                                                                                                                              >  5:311566/1‑148 (MQ=255)

                                                                                                                                                                                                                                   aaaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGaa                                                                                                                                                                           <  7:830673/149‑1 (MQ=255)

                                                                                                                                                                                                                                          gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                                                                                                                                                                                                          <  8:838494/111‑1 (MQ=255)

                                                                                                                                                                                                                                          gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGc                                                                                                                                                                    <  8:907007/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGc                                                                                                                                                                    >  5:478201/1‑148 (MQ=255)

                                                                                                                                                                                                                                                  aaGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATcc                                                                                                                                                                                                            >  4:607970/1‑101 (MQ=255)

                                                                                                                                                                                                                                                   aGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGc                                                                                                                                                                                                        <  6:86060/104‑1 (MQ=255)

                                                                                                                                                                                                                                                   aGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCa                                                                                                                                                           >  8:257858/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAc                                                                                                                                                          >  7:6903/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                                                                                                                                                                                                           <  7:229300/81‑1 (MQ=255)

                                                                                                                                                                                                                                                       aacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                                                                                                                                                                                                           >  8:229300/1‑81 (MQ=255)

                                                                                                                                                                                                                                                          aTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtgatg                                                                                      >  2:206560/1‑210 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATa                                                                                                                                                                         >  4:744560/1‑123 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGt                                                                                                                                                                 <  6:815052/131‑1 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTg                                                                                                                                               >  6:506244/1‑148 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTg                                                                                                                                               >  6:171411/1‑148 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGc                                                                                                                                              <  8:803065/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCATCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTg                                                                                                                                               >  6:269569/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                   gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATg                                                                                                                                                                     <  2:744562/123‑1 (MQ=255)

                                                                                                                                                                                                                                                                   gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCg                                                                                                                                   >  2:319152/1‑156 (MQ=255)

                                                                                                                                                                                                                                                                   gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCCATATCAGGCAGTTTTTAGCGATGtt            >  2:436240/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGAACGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTCATGGCTGAGTTCAGCAGTTTTACTGAGATCATGTCCCGGCCATAACAGGACTTTTTTAGCGTTTTTa           >  4:212857/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTAccc                                                                                                                                         >  8:347418/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTGTACAAATGTACACAATCGCAGAGAGAGGAAATGACATGCGTCGCGCTACCGATCCGCAACGGCGGGAAAAAAATATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATTCTGTTACGCACCGCA‑AAATTGCTAccc                                                                                                                                         >  8:293923/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGATAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTg                                                                                                                                      >  7:512430/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTg                                                                                                                                      >  7:387425/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                                                                                                                                                                                                     <  7:267545/53‑1 (MQ=255)

                                                                                                                                                                                                                                                                         tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                                                                                                                                                                                                     >  8:267545/1‑53 (MQ=255)

                                                                                                                                                                                                                                                                         tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGc                                                                                                                                     <  5:536536/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                           cATTTGTACACAATTGCAGAGCGAGGAATTGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                                                                                                                                                       <  2:163811/65‑1 (MQ=255)

                                                                                                                                                                                                                                                                             tttGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGaa                                                                                            >  2:620085/1‑185 (MQ=255)

                                                                                                                                                                                                                                                                             tttGTACACAACTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCggg                                                                                                                                 >  5:83360/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                              ttGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCgg                                                                                                                                  <  5:172540/146‑1 (MQ=255)

GTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCA‑AAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAGTGATGCTC  >  NC_000913/887123‑887666

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`