BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	314,258	G→A	5.0%	intergenic (‑453/‑99)	ykgP ← / → eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	314,258	0	G	A	5.0%	595.7 / 10.6	176	intergenic (‑453/‑99)	ykgP/eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand): ref base G (90/77); new base A (5/4); total (95/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 178 total aligned reads displayed.

CTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTT  >  NC_000913/314006‑314520

                                   gCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTgtgt                                                                                                                                                                                                               <  4:21416/275‑1 (MQ=255)

                                            tGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCATTTTCTGCGGAAAAGGAAATCGGGAAATCACCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTAt                                                                                                                                                                                                       >  1:6324/1‑274 (MQ=255)

                                                      tctGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGACAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCAATTATTTGTGTCTGCCTATGTTCGTTAAtt                                                                                                                                                                                            >  2:231471/1‑275 (MQ=255)

                                                                           tGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtctc                                                                                                                                                                       >  2:131886/1‑275 (MQ=255)

                                                                                             tataACTTAGGTATCTATTTAATAACTTGACCAAAAAACTAAATTTGTCCCCATAAATAAAAATATAAGCCCGGGCCCAACCACCGGATAAGTGGCTATAGGAACTGGATATATTAGGTGGCGGTGCACTGTCGTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                     <  4:149355/275‑1 (MQ=255)

                                                                                                            tATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                                                                                                                                                    >  5:868110/1‑149 (MQ=255)

                                                                                                            tATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaa                                                                                                                                                                                                                                                                    >  6:875720/1‑149 (MQ=255)

                                                                                                             aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                                   <  5:1429098/149‑1 (MQ=255)

                                                                                                             aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                                   <  6:1461954/149‑1 (MQ=255)

                                                                                                                tAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTCTCCGCAGAAAAGGa                                                                                                                                                                                                                                                                <  5:603793/149‑1 (MQ=255)

                                                                                                                tAATCAATTACAAAAAAAACTTNNTTTTCNCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGa                                                                                                                                  <  3:6324/275‑1 (MQ=255)

                                                                                                                  aTGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaaa                                                                                                                                                                                                                                                              <  7:338580/149‑1 (MQ=255)

                                                                                                                      cTTACACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                          <  5:470179/149‑1 (MQ=255)

                                                                                                                            caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATc                                                                                                                                                                                                                                                    >  5:1331481/1‑149 (MQ=255)

                                                                                                                             aaaaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGAAAtacc                                                                                                                                                                                                                                                          <  7:185465/142‑4 (MQ=255)

                                                                                                                               aaaGCTACATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTCTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                                                                                                                                                 >  8:918117/1‑149 (MQ=255)

                                                                                                                               aaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                                                                                                                                                 <  7:1098187/149‑1 (MQ=255)

                                                                                                                               aaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcccc                                                                                                                                                                                                                                                 <  6:989672/149‑1 (MQ=255)

                                                                                                                                   cTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGtt                                                                                                                                                                                                                                             >  6:770060/1‑149 (MQ=255)

                                                                                                                                    tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttt                                                                                                                                                                                                                                            >  8:393911/1‑149 (MQ=255)

                                                                                                                                     aaaTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGtttt                                                                                                                                                                                                                                           >  6:403940/1‑149 (MQ=255)

                                                                                                                                      aaTTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttttt                                                                                                                                                                                                                                          >  8:1253139/1‑149 (MQ=255)

                                                                                                                                        ttttCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCt                                                                                                                                                                                                                                        <  6:429786/149‑1 (MQ=255)

                                                                                                                                            cccccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTgcg                                                                                                      <  2:319256/275‑1 (MQ=255)

                                                                                                                                             ccccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATACCCGTTTTTTCTGACaa                                                                                                                                                                                                                                   >  5:536462/1‑149 (MQ=255)

                                                                                                                                               ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                                                                                                                                                 >  5:386658/1‑149 (MQ=255)

                                                                                                                                                cATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                          >  7:582982/1‑123 (MQ=255)

                                                                                                                                                 ataaataaAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGa                                                                                                                                                                                                                                      <  8:236953/142‑1 (MQ=255)

                                                                                                                                                 ataaataaAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                                                                                                                                               >  8:423917/1‑149 (MQ=255)

                                                                                                                                                  taaataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGa                                                                                                                                                                                                                              <  6:585719/149‑1 (MQ=255)

                                                                                                                                                     ataaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGc                                                                                                                                                                                                                           >  5:443545/1‑149 (MQ=255)

                                                                                                                                                      taaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                                                                                                                                                                                          >  5:567943/1‑149 (MQ=255)

                                                                                                                                                      taaaAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGcc                                                                                                                                                                                                                          >  5:985924/1‑149 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                                                                                                                                       >  8:1452581/1‑149 (MQ=255)

                                                                                                                                                         aaaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatt                                                                                                                                                                                                                       >  6:138230/1‑149 (MQ=255)

                                                                                                                                                             ataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTAttt                                                                                                                                                                                                                   <  7:651767/149‑1 (MQ=255)

                                                                                                                                                               aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                                                                                                                                                 <  8:815441/149‑1 (MQ=255)

                                                                                                                                                               aaTCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTtgt                                                                                                                                                                                                                 >  7:486282/1‑149 (MQ=255)

                                                                                                                                                                 tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                    >  4:359171/1‑192 (MQ=255)

                                                                                                                                                                   ccGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCt                                                                                                                                                                                                             >  8:571467/1‑149 (MQ=255)

                                                                                                                                                                     gcgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                                <  2:359171/192‑1 (MQ=255)

                                                                                                                                                                         ccAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTAt                                                                                                                                                                                                       <  6:518071/149‑1 (MQ=255)

                                                                                                                                                                           aaCCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGt                                                                                                                                                                                                     <  6:805372/149‑1 (MQ=255)

                                                                                                                                                                            accaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTATGCCTATGtt                                                                                                                                                                                                    >  5:798218/1‑149 (MQ=255)

                                                                                                                                                                             ccaccTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTc                                                                                                                                                                                                   <  5:138230/149‑1 (MQ=255)

                                                                                                                                                                               accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                                                                                                                                                                                                 <  8:1284038/149‑1 (MQ=255)

                                                                                                                                                                               accTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGt                                                                                                                                                                                                 >  8:830166/1‑149 (MQ=255)

                                                                                                                                                                                   gatgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAAt                                                                                                                                                                                             <  6:1094341/149‑1 (MQ=255)

                                                                                                                                                                                     tgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTc                                                                                                                                                                                           <  7:571467/149‑1 (MQ=255)

                                                                                                                                                                                      gaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGt                                                                                                                                                                                                     >  6:1280925/1‑138 (MQ=255)

                                                                                                                                                                                        gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                          >  6:1086207/1‑83 (MQ=255)

                                                                                                                                                                                        gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCgg                                                                                                                                                                                                                                                          <  5:1086207/83‑1 (MQ=255)

                                                                                                                                                                                            cTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGGCTATGTTCGTTAATTCGTtcatc                                                                                                                                                                                    <  7:220705/149‑1 (MQ=255)

                                                                                                                                                                                              ataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCt                                                                                                                                                                                                                                        >  7:115763/1‑95 (MQ=255)

                                                                                                                                                                                              ataGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAg                                                                                                                                                                                  <  5:875720/149‑1 (MQ=255)

                                                                                                                                                                                                    aCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTGCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACCAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGt                                                                                                                                                                                         <  6:876372/136‑1 (MQ=255)

                                                                                                                                                                                                     cTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTa                                                                                                                                                                           >  7:721991/1‑149 (MQ=255)

                                                                                                                                                                                                      tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCa                                                      >  2:262550/1‑265 (MQ=255)

                                                                                                                                                                                                       ggATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc                                                                                                                                                                         <  6:619357/149‑1 (MQ=255)

                                                                                                                                                                                                          tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCa                                                                                                                                                                      >  8:567678/1‑149 (MQ=255)

                                                                                                                                                                                                          tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCa                                                                                                                                                                      >  7:983166/1‑149 (MQ=255)

                                                                                                                                                                                                          tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCtt                                                  <  4:262550/265‑1 (MQ=255)

                                                                                                                                                                                                           atatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                                                                                                                                     <  5:777240/149‑1 (MQ=255)

                                                                                                                                                                                                            tatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                    <  8:983166/149‑1 (MQ=255)

                                                                                                                                                                                                            tatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                    >  7:1226763/1‑149 (MQ=255)

                                                                                                                                                                                                            tatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGGTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                                    >  7:448652/1‑149 (MQ=255)

                                                                                                                                                                                                             atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                                   <  8:67810/53‑1 (MQ=255)

                                                                                                                                                                                                             atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGa                                                                                                                                                                                                                                                                <  6:33805/56‑1 (MQ=255)

                                                                                                                                                                                                             atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGa                                                                                                                                                                                                                                                                >  5:33805/1‑56 (MQ=255)

                                                                                                                                                                                                             atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                                                                                                                                              >  7:1027446/1‑122 (MQ=255)

                                                                                                                                                                                                             atTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTaa                                                                                                                                                                                              <  8:1027446/122‑1 (MQ=255)

                                                                                                                                                                                                                 ggTTGCGGTGCACTTTCTTACCAAAAGGTATTTCCGTTTCTGCGGAAAAGGAATTCGGGAAATCCCCGGTTTTTCGGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTACTTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATggc                                   <  2:9871/273‑1 (MQ=255)

                                                                                                                                                                                                                  gTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCTACTCGCTGTCACCTTTACCCCAGTAATggcgg                                 >  1:40503/1‑274 (MQ=255)

                                                                                                                                                                                                                      cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCCCATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACCCGCTGTCCCCTTTACCCCAGTAATGGCGGCACGt                            >  4:48412/1‑275 (MQ=255)

                                                                                                                                                                                                                        gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                                                                                                                                        <  7:361497/149‑1 (MQ=255)

                                                                                                                                                                                                                        gTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaa                                                                                                                                                        >  6:468704/1‑149 (MQ=255)

                                                                                                                                                                                                                            aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                                                                                                                                    <  6:1022527/149‑1 (MQ=255)

                                                                                                                                                                                                                            aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                                                                                                                                    >  8:509797/1‑149 (MQ=255)

                                                                                                                                                                                                                            aCTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACa                                                                                                                                                    <  7:1408456/149‑1 (MQ=255)

                                                                                                                                                                                                                              tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTtcatc                                                                                                                                                                                    <  5:769813/115‑1 (MQ=255)

                                                                                                                                                                                                                              tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                                  >  6:226288/1‑149 (MQ=255)

                                                                                                                                                                                                                              tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                                  >  7:1035256/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAaacaac                                                                                                                                       >  5:145483/1‑149 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAAcc                                                                                                                                      <  6:1471236/149‑1 (MQ=255)

                                                                                                                                                                                                                                          aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTTGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCAACTCGCTGTCACCTTTCCCTCAGTAATGGCGGAACGTGCGCAGCATGCGGTTCAGcc        >  2:238678/1‑275 (MQ=255)

                                                                                                                                                                                                                                             gaTTTCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACg                                                                                                                                   >  8:120985/2‑149 (MQ=255)

                                                                                                                                                                                                                                              aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGa                                                                                                                                  >  5:1257557/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAt                                                                                                                                 >  6:764785/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAt                                                                                                                                 >  7:527144/1‑149 (MQ=255)

                                                                                                                                                                                                                                               tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCCCTGTTCAGGTTCTTTTTCCACTTGTTGTCCCCTTTACCCCAGTAATGGCGGCACGTGCCGAGCATGCGGTTCATCCACGGt   >  4:334582/1‑275 (MQ=255)

                                                                                                                                                                                                                                                ttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGc                                                                                                                                                                                                           >  2:161900/1‑74 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGt                                                                                                                                                 >  8:1003647/1‑131 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAttt                                                                                                                               >  6:837633/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tCCTTTTCCGCAGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGAttt                                                                                                                               >  6:1222332/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTGTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt                                                                                                                            >  7:1376185/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt                                                                                                                            >  8:1291717/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGt                                                                                                                            >  5:701186/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtat                                                                                                                                                          >  2:79670/1‑116 (MQ=255)

                                                                                                                                                                                                                                                            ggAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGt                                                                                                                    >  7:1276289/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            ggAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGt                                                                                                                    >  5:100504/1‑149 (MQ=255)

CTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCATAAACAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTT  >  NC_000913/314006‑314520

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