Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,276	0	T	C	23.5%	127.7 / 12.0	72	pseudogene (885/2518 nt)	ydbA	putative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base C (13/3);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 9.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGC  >  NC_000913/1466026‑1466522                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt                                                                                                                                                                                                                                                         <  1:90664/251‑1 (MQ=255)      aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGACCCGGAGTCCATGGGCATCCAGAACGACGGCG‑Acc                                                                                                                                                                                                                                                    >  2:93037/1‑249 (MQ=18)       aCGGTGAGCTTGTTATTACGGGTGATTATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCG‑Accag                                                                                                                                                                                                                                                  <  2:56594/251‑4 (MQ=11)         ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGGATCCAGGTCGACGGCG‑ACCAGGc                                                                                                                                                                                                                                                >  2:420299/1‑251 (MQ=255)         ccacgctcttccgataagggTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCG‑ACCAGGc                                                                                                                                                                                                                                                <  2:123756/234‑1 (MQ=11)                     attaCGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGGGATGAAAGACGGCGAGCTGGATGGCGGAGGCCGCGGTCGCGGTATTGTTATTACAGACGACAGCGCGGAGGGGGATAACAAGGGCACAATGACCGTCACCGCTACGGAGGCCATCGGAATCCAGATCGACGGTA‑ACAAggacaacgtgaaca                                                                                                                                                                                                                                     >  2:387547/1‑238 (MQ=12)                       taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATTGACGGCG‑ACAAGGCCATCGTCAATAACg                                                                                                                                                                                                                                  >  1:221319/1‑251 (MQ=255)                         cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGCACCGGAGTCACTGGGTATCCAGATCGACGGTG‑ACCAGGCCATCGTCAATAACgc                                                                                                                                                                                                                                 >  1:263511/1‑249 (MQ=18)                          gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCG‑ACCAGGccgaccgcaagagcgcac                                                                                                                                                                                                                               >  1:123756/1‑235 (MQ=11)                               aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCG‑ACAAGGCCGTCGTCAATAACGAAGGCg                                                                                                                                                                                                                            >  1:263738/1‑249 (MQ=255)                                 aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACCAGGCCGTCGTCAATAACGAAGGCgag                                                                                                                                                                                                                          >  1:323655/1‑249 (MQ=255)                                          gaaaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAGCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                                                                                                                                                                                                               <  2:441502/249‑1 (MQ=255)                                                aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAGTAACGGGAAAGGGATTCAGGACGGCGACCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGCCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGG‑ACAAGGCCATCGTAAATGACGAAGGCGAAAGCCACATCACCAAc                                                                                                                                                                                                           >  2:67778/1‑249 (MQ=18)                                                            accaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcac                                                                                                                                                                                             <  2:356625/251‑1 (MQ=255)                                                            accaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGAATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcac                                                                                                                                                                                             <  1:157353/251‑1 (MQ=255)                                                            accaccGTGTACGGTAAAGATGCCACCGGTACGGAAATCAACTGTAATAACGGGAAAGTGATTCAGGGCGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCGCCATGACCGGCACCGATCCGGAGTCCATGGGTAGCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcac                                                                                                                                                                                             <  2:14973/251‑1 (MQ=18)                                                             ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                                                                                            <  1:316790/251‑1 (MQ=255)                                                                ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                                                         <  1:334248/251‑1 (MQ=255)                                                                 cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGCTGTGGATGTCAGCGGCGGCAGTCACGGTATTGATATCAGCGGGGACAGCCCTACGGAGGATAACAAGGGCACCCAGACCGTCACCGATCCGGAGTCCATAGGTATCACGACTGACGGGGGACA‑GGCGGTCGTAAACAACGAAGGCGAGCGCACCATCACCAACGGTggcaccggcaacc                                                                                                                                                                                           >  2:454790/1‑245 (MQ=2)                                                                    tGACGGTAAAGAGTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                                                                                                                                                      <  2:291173/250‑1 (MQ=255)                                                                     gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTCGACGGCG‑ACACGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCGAGATTa                                                                                                                                                                                     >  1:162219/1‑250 (MQ=255)                                                                           aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACGGAGCGCCGGCGGCTGACATCGGGACCATGACCGTCACCGATCCGGGGTCCATGGGTGTCCAAATCGACGGCG‑ACAAGACCATCGTCAATAACGCATGCGAGAGAACCCTCACCAACGGTGGGACCGGCACCCCGGTTAAtag                                                                                                                                                                                 >  1:136274/1‑246 (MQ=17)                                                                             aGATTCCGCCGGTACGGAAATCAACGGTAATAACGGGTAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGTCGAGAGTACCATCACCAACGGTGGCTCCGGCACCCAGATTAATGGtgatg                                                                                                                                                                            <  1:253200/251‑1 (MQ=255)                                                                                 tCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGCG‑ACAAGGCCGTCGTTAATAACGAAGGCGAGAGCGCCATCACCAAAGGTGGGAACGGCACCCAGATTAACGCTAATGACGc                                                                                                                                                                        >  1:85275/1‑251 (MQ=14)                                                                                   cACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGcc                                                                                                                                                                       >  1:357686/1‑250 (MQ=255)                                                                                     ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACg                                                                                                                                                                    <  2:192945/251‑1 (MQ=255)                                                                                            ggAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGGCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                                                              >  1:52795/1‑250 (MQ=255)                                                                                             gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCGCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCATCGGTTGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                                                                                                                                            <  1:384871/251‑1 (MQ=255)                                                                                                    aCGGTAATAACGGGAAAGTGATTCAGGACGGCGACCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                                                                     >  1:48264/1‑251 (MQ=255)                                                                                                       gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAaacaa                                                                                                                                                            >  1:113370/1‑241 (MQ=255)                                                                                                            aaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGCCGGGG‑ACCAGGCCGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCCCGGCCCACAACAACGGCAAAacacc                                                                                                                                             >  2:415294/1‑248 (MQ=17)                                                                                                                gtaacgTGATGCAGGACGTCGATGTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGGGACGGTGGATAACAAGGGCACCATGACCGTCACCGGGCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAATAATGGTGATGACGCCACGGCAAACAACAACGGCAAAacctcagtt                                                                                                                                         <  2:72037/246‑7 (MQ=18)                                                                                                                ggAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaaatccggtg                                                                                                                                         >  1:104612/1‑243 (MQ=255)                                                                                                                 gAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACCAGGGCACCATGACCGTCACCGTTCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAACGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCCAACAACAACGGCAAAACCACCGTTg                                                                                                                                        >  1:422328/1‑251 (MQ=255)                                                                                                                            aGGACGGCGAGCTGGGTGGCAGCGGCGGCGGTCACGGTATTGATAGCACCGGTGACAGCGCGACGGTGGATAACAAGGGTACCAGGACCGTCACCGATCCGGAGTCCATAGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGTCAAAACCACCGTTGACGGCAAGGAt                                                                                                                             <  2:69851/251‑1 (MQ=21)                                                                                                                                     aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGt                                                                                                                    >  1:307778/1‑251 (MQ=255)                                                                                                                                         ggATGTCAGCGGCGGCGGTCACGGTATTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCGGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                                                                                    <  2:48264/247‑1 (MQ=255)                                                                                                                                         ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGCGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                                                <  1:91270/251‑1 (MQ=255)                                                                                                                                          gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                                                               <  1:67778/251‑1 (MQ=255)                                                                                                                                              gcAGCGGCGGTGGTGACGGTATTGATAGCACCGGTGCCCGCGAGACGGTGGATAACAAGGGGACCATGACCGTCACCGATCCGGAGTCCATGGGTAAACAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                                                           <  2:357777/250‑1 (MQ=18)                                                                                                                                                 gcggcggcgGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCGCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                                                        <  1:130797/251‑1 (MQ=255)                                                                                                                                                        cggTCACGGTATTGATATCACCGGTGACATCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaata                                                                                                 <  1:359680/251‑1 (MQ=255)                                                                                                                                                           tCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGCG‑ACAAGGCCGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGTCGGCAGCGGGTCGCAGATTGAAGGGGAGGACGCCACGGCGGACAACGACGCAAAAAACACCGCTGACGGCAAGGACTCCACCGGGCCGGAAATCAACGGTAATAACg                                                                                              >  2:253627/1‑251 (MQ=14)                                                                                                                                                                    tGATATCTCCGGAGAACGCGCGAAGATGGATGGCAAGGGCACCACGACTATCCCCGATCCGGAGGACATGGGGATCAAGATCGATGGCC‑ACAAGGCCATCGTCAATGACGAAGGCGAGAGCACCCTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtt                                                                                      <  1:252791/250‑1 (MQ=18)                                                                                                                                                                        atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCG‑ACCAGGCCATCGTCCATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACgggaaggtgattc                                                                                  >  2:194193/1‑239 (MQ=14)                                                                                                                                                                         tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCG‑ACAAGGCAGTGGTCACTAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGTCCGCTGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                  >  1:272049/1‑249 (MQ=17)                                                                                                                                                                              ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGCTTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                                                           >  1:23947/1‑251 (MQ=255)                                                                                                                                                                                   cAGCGCGACTGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGCTATTCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCAACACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                                                      <  1:181695/251‑1 (MQ=255)                                                                                                                                                                                         gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGCATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                >  1:209358/1‑251 (MQ=255)                                                                                                                                                                                            ggTGGATAACAAGGGCACCATGACCGTCACCGGGCCGGAGGCTATGGGTATCCAGATCGACGGTG‑ACAAAGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCGCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                                                             <  1:279020/251‑1 (MQ=21)                                                                                                                                                                                             gTGGATAGCAAGGGCACCATGACCGTCTCTGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGGAGGCGAGAGGACCCACACCAACGGTGGCAGCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAATTCAACGGTAATAACGAAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                                                            <  1:243535/251‑1 (MQ=21)                                                                                                                                                                                                   aaCAAGGTCACCATGACCGCCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcg                                                       <  2:307778/250‑1 (MQ=255)                                                                                                                                                                                                      aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                   <  1:439669/251‑1 (MQ=255)                                                                                                                                                                                                       aGGGCACCATGACCGACCCCGATCCGGAGTACATGGGTATCCAGACCGACGGTG‑ACAAGGGCATCGTCAATAACGAAGGCGAGAGCTCCAACTCCAACCGTGGCACCGGCACCCAGATTAATGGTGATAACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                  <  1:199820/251‑1 (MQ=21)                                                                                                                                                                                                        gggCCCCATGACCGTCACCGATCCGGAGTCCATGGGTAACCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGTCACCCAGATTAATTGTGATTACGGCACCGCCAACAACAGAGACAAACCCTCCTTCCCCACCCAGCGCTCCACCGATCTGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc                                                 <  2:52795/251‑1 (MQ=19)                                                                                                                                                                                                           cACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAcc                                              >  1:321807/1‑250 (MQ=255)                                                                                                                                                                                                                tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGAGAAGGATTCAGCCGACACGGGAATCCACGTTAATAAGGCAACAGGTATTCAGGACGCCGATCTGAAAGTCAGCGGCGGCGGTAACGGTAtt                                         >  1:5415/1‑251 (MQ=18)                                                                                                                                                                                                                 gtccGTCCCCGATCCGGAGTACATGGGTATCCAGAGCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCCCCCAGATTACTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                        <  1:34266/249‑1 (MQ=21)                                                                                                                                                                                                                 gACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                        >  1:336579/1‑251 (MQ=255)                                                                                                                                                                                                                   ccGTCTCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAAGAACGAAGGCGAGAGCACGATCACCAAGGGTGGCACCGGCGCCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                      <  1:210018/251‑1 (MQ=255)                                                                                                                                                                                                                     gTCACCGCTCCGGAGTCCATGGGTATCCAGATCGACGGAG‑ACAAGGCCGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCCAATAACAGCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGAAAATCAACGGCAATAACGGGAATGTAATTCAGGAGGGCGGTCTGNGTGTAAgcggcggcgagcctggtattgat                                      >  2:215097/1‑235 (MQ=17)                                                                                                                                                                                                                      tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGCAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGtagagcggaaggacgc                                   >  2:267579/1‑236 (MQ=18)                                                                                                                                                                                                                      tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTAATGACGCCACCGCAAACAACAACGGCAAAACCACCGTCGACGGCAAGGCTTCTCCCGGTACGGAAATCCATGGCAATAACGGCAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGAGGTCACGGCATTg                                        >  1:36469/1‑246 (MQ=18)                                                                                                                                                                                                                              ccGGAGGCCATGGGTATCCAGATCGACGGTG‑ACAAGGGCATCGTCAATAACGAAGGCGAGAGCACCATCTCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCg                            <  1:279959/250‑1 (MQ=255)                                                                                                                                                                                                                                 gAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCNGCGGTCACGGTATTGATATCACCGGCggc                          >  2:231935/1‑247 (MQ=255)                                                                                                                                                                                                                                  aGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                       >  1:115679/1‑251 (MQ=255)                                                                                                                                                                                                                                         gggtgtcaggagcGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCAAGAGCACCATCGCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCGCGGCAAACAACAACGAAAAAACCACCGTTGACGGCAAGGAACCCTCCGGTACGGAAATCAGCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACg                 <  2:422328/238‑1 (MQ=18)                                                                                                                                                                                                                                                  gATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATACCGGAAAGGTTCTTCAGGACGGCGATCTGGTTGTCAGCGGCGGCGGGCACGGCATTGATATCACCCGCGACAGCGCAACGGTggctcac        >  1:66757/1‑245 (MQ=18)                                                                                                                                                                                                                                                  gATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGCGACGGTGCTGATATCACCgcccgcagcacataggcagatagc        >  1:73752/1‑227 (MQ=18)                                                                                                                                                                                                                                                      gACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCAACCAGATTAATGGAGATGTCGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGAGCCCGCGCAACGGGGGATAACAAggg   >  2:282179/1‑251 (MQ=17)                                                                                                                                                                                                                                                       gcGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGc  <  2:336579/250‑1 (MQ=255)                                                                                                                                                                                                                                                       aCGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAAGACGACCTGGATGTAAGCGACGGCGGTCAGCGTATTGATATCACCGGTGACAGCGCAACGAAGGATAACAAgggg  >  1:18591/1‑250 (MQ=18)                                                                                                                                                                                                                                                       aCGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCAACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGGCAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCCCGGTATTGATATCACCGGCGACAGCGCAAAGGTGGATAACAAGGGc  >  2:414294/1‑251 (MQ=14)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        ATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTG‑ACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGC  >  NC_000913/1466026‑1466522

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑