Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,438	0	T	G	29.6%	74.2 / 26.3	54	pseudogene (1047/2518 nt)	ydbA	putative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base G (7/9);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAC  >  NC_000913/1466161‑1466704                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                             tGGATGTCAGCGGCGGCGATCACAGTATTGATAGCACCGGTGACAGCGAAACGAGGAATAACAGGGGCACCATGACTATCACCGATCCGGAGTCCATGGGGATCCAGATCGACGGCGACAAGGCCATCGTCAAGAACGAAGGCGAGAGCACTAGCACCAACAGTGGCACCGGCACCCGGATTAATGGTGATGACGCCACGGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                                                                                                                                                                        <  2:79074/300‑1 (MQ=18) tGGACGTCATCGTCGGCGTTCAAGGTATTGACATCACTGTTGATAGTCCGACGGTGGATAACAAGGATACCACGACCGTCACCGGTCCGGAGGCCATGGGTATCCAGATCGACGTTGACGAGTCCGTCGTCAATAACGAAGTCGAGAGCACCATCGCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                                                                                                                                                                        <  2:79691/300‑1 (MQ=18)   gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCANCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGTAATAAAGGAAAGGTGAGTTAGGACGGCGATCTGGATGTc                                                                                                                                                                                                                                                      >  2:237766/1‑300 (MQ=255)                       cGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                                                                                                                                                  >  2:722719/1‑300 (MQ=255)                       cGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                                                                                                                                                  >  1:263457/1‑300 (MQ=255)                          tATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa                                                                                                                                                                                                                               <  1:415132/300‑1 (MQ=255)                                 atCACCGTTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                                                                                                                                                        <  1:755024/300‑1 (MQ=255)                                                  gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCAGCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACgg                                                                                                                                                                                                        <  1:170268/299‑1 (MQ=255)                                                      gTGGATAACAAGGGCACCATGACCGTCACCGCTCCGGAGTCCATAGGTATCCAGATCGACGGTGACAAGGACATCGTCAATAACGAAGGCGAGAACACCATCACCAACGGTGGCACCGGCAACCAGACTAATGGTGATGACGCCAAGGCAAACAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTAAGGAAATCAATGGNAATAACGGCAAAGTTAGTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCActagcgac                                                                                                                                                                                                                  >  2:596180/1‑278 (MQ=17)                                                      gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCNACAACGGTGGCACCGGCACCCAGATTAATGGTNATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                   >  2:571595/1‑300 (MQ=255)                                                        ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATaa                                                                                                                                                                                                 >  2:712254/1‑300 (MQ=11)                                                         gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                                                                                <  1:24501/300‑1 (MQ=255)                                                         gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGAcctttccattataaacacc                                                                                                                                                                                                >  2:502834/1‑282 (MQ=18)                                                           tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                                                                                              <  1:712254/300‑1 (MQ=11)                                                            aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGTTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑g                                                                                                                                                                                            <  1:250651/300‑1 (MQ=255)                                                                    cACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACaatagccccct                                                                                                                                                                                     >  1:683030/1‑291 (MQ=255)                                                                      ccATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATCAATGGTGATGACGCCACGGCAAACACCAACGGCCAAACCACCATTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGACGACAGCGCAACGGTGGATAACAA‑GGGCACCAt                                                                                                                                                                                    >  1:287621/1‑298 (MQ=255)                                                                         tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCg                                                                                                                                                                               <  1:664700/300‑1 (MQ=255)                                                                                         ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAg                                                                                                                                                                <  1:302308/298‑1 (MQ=255)                                                                                            gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCg                                                                                                                                                                         <  1:501508/287‑1 (MQ=11)                                                                                            gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTACTGGTAATAACGGAAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGAtgtttcttatac                                                                                                                                                            >  2:501508/1‑289 (MQ=14)                                                                                            gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAACGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCa                                                                                                                                                            >  1:372671/1‑299 (MQ=11)                                                                                            gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACAGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTNCCGGTAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAACGGCGGCGGTCACGGTATTGATATCCACGGCGAAAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCAcc                                                                                                                                                                          >  2:242851/1‑286 (MQ=14)                                                                                             tCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCCGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATc                                                                                                                                                          <  1:571595/300‑1 (MQ=255)                                                                                                aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCgg                                                                                                                                                        <  1:524022/299‑1 (MQ=255)                                                                                                   ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAAAGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAACGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCCCCGATCC‑GGAGTCCATTGGTATc                                                                                                                                                    >  2:735987/1‑300 (MQ=255)                                                                                                                         aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCANCGTTGACGGCAAGGATTCCACCGGTACGGAAANCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGATTGACGGCGACCAGGccc                                                                                                                              >  2:577821/1‑298 (MQ=255)                                                                                                                         aGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGATTGACGGCGACCAGGcc                                                                                                                               <  1:577821/299‑2 (MQ=255)                                                                                                                                 gTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTAACGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAAGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCGCGGTATTGATATCAACGGTGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCAACGATCC‑GGAGTCCATCGGTATACAGGTTGACAGCGAACAGGCGGTTGTTaac                                                                                                                      >  2:174230/1‑300 (MQ=17)                                                                                                                                            agtcGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                                                                                                                                             <  2:183611/185‑1 (MQ=255)                                                                                                                                             ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGCAATAACGGAAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACTAGGCCgttgtcaataaactgtctagt                                                                                                          >  2:77356/1‑284 (MQ=11)                                                                                                                                                  gagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATANCAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAgcac                                                                                                     >  2:582260/1‑298 (MQ=255)                                                                                                                                                          tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                                                             <  1:582260/300‑1 (MQ=11)                                                                                                                                                            gataactgtggcatcggcatttgtaTTAAGGGTGATGACGCCACGGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTCGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAc                                                                                           <  1:496268/276‑1 (MQ=17)                                                                                                                                                                       caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaat                                                                                           >  2:747731/1‑288 (MQ=11)                                                                                                                                                                           gatccttatattacgggTGATGATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAc                                                                                           <  2:359266/270‑1 (MQ=11)                                                                                                                                                                                   gATTAATGTTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCa                                                                                             <  1:113675/275‑1 (MQ=11)                                                                                                                                                                                   gATTAATGGTGATGACGCCACGGCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                                                                               <  2:667109/241‑1 (MQ=255)                                                                                                                                                                                   gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGAAGGCGAACAGGCAATTGTCAATAACGAAGGCGAGAgcaacatcaaaa                                                                                             >  2:112925/1‑265 (MQ=14)                                                                                                                                                                                   gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAC‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaatgtttcttatacacatctcgca                                                                      >  2:113675/1‑276 (MQ=11)                                                                                                                                                                                    aTTAATGGTGATGACGCCACGGTAAACAACAACGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑ggg                                                                                                                                                                                          <  1:705131/182‑1 (MQ=255)                                                                                                                                                                                      tacgggTGATAATGCCACGGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                 <  2:477100/296‑1 (MQ=14)                                                                                                                                                                                                      aCGGTAAACAACAATGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGc                                                  <  1:685856/299‑1 (MQ=255)                                                                                                                                                                                                                         aaaCCACCGTTGACGGCAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAc                                                                                               <  2:320038/235‑1 (MQ=255)                                                                                                                                                                                                                                    gACGGTAAGGATTCCACCGGTACGGAAATCAACGTTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGTTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                   <  2:96009/300‑1 (MQ=255)                                                                                                                                                                                                                                      cGGCAAGGATTCCAACGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGAACGTCAACGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGACAGTCGTTAATAACGAAGGCGAGAGCGCCATCACCAACGATGGCACCGGCACGCAGATTAACGGTGATGGCGACACGGCAAACAAAAATGACCAAACCACCGttaccggca                  >  2:254900/1‑292 (MQ=17)                                                                                                                                                                                                                                           aGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt            <  1:310760/300‑1 (MQ=255)                                                                                                                                                                                                                                            ggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAATTGTCAACAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACGCAGATCAACGGTAACGACGCAACAGCAAACAACAATGGAAAAACCACTGTTGATGGAAAAGATTc           >  1:255941/1‑300 (MQ=11)                                                                                                                                                                                                                                            ggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCAACGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGTGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAAGGCAAACAACAACGGCAAAACCAACGTTGAAGGCAAGGAtt            >  2:16269/1‑299 (MQ=255)                                                                                                                                                                                                                                            agATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAga                                                                                                                                                >  2:648774/2‑166 (MQ=9)                                                                                                                                                                                                                                            agATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAga                                                                                                                                                <  1:648774/166‑2 (MQ=9)                                                                                                                                                                                                                                             gATTCCACCGGTACGGNAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAGTTGTCAACAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACGCAGATCAACGGTAACGACGCAACAGCAAANAACAACGGAAAAACCACCNTTGACGGAAAAGATTcc          >  1:640045/1‑300 (MQ=14)                                                                                                                                                                                                                                                      gggaaggagataaaCGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc        <  1:379812/281‑1 (MQ=255)                                                                                                                                                                                                                                                       gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATC‑CAACGAGCACCCCGACCGCCAACGATCCAGCACACC‑TCAGCCTCCAAATCCACGGCAAACATGTAGCTTTCAACAACAAAAGCACAAGCACAATCAACAACGGTGGAACAGGCAAGCAGATTAACGGTGATGACGCCAAGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCAACAGGTAc  >  2:720274/1‑298 (MQ=9)                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                             TGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAA‑GGGCACCATGACCGTCACCGATCC‑GGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAC  >  NC_000913/1466161‑1466704

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