Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,228,239	0	G	C	35.3%	62.4 / 14.0	51	intergenic (‑969/‑743)	yohF/yohP	putative oxidoreductase YohF/putative membrane protein YohP
Reads supporting (aligned to +/- strand):  ref base G (13/20);  new base C (8/10);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCT  >  NC_000913/2228001‑2228484                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                  tttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGGCCGTGCTTTa                                                                                                                                                                                                                                                        <  2:516079/250‑1 (MQ=255)  ttcCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑TTCC‑Gcccccccggcttttac                                                                                                                                                                                                                                                      >  1:568984/1‑241 (MQ=255)    ccccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTG‑CTCGC‑G‑TTTC‑GCCGcccccgcttttaccc                                                                                                                                                                                                                                                    >  2:346023/1‑239 (MQ=255)    ccccGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑g                                                                                                                                                                                                                                                   <  2:591315/251‑1 (MQ=255)                   gAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑TTCC‑GCCGCCCCGGGTTTAACC‑CGA‑CCCaggcaggccg                                                                                                                                                                                                                                   >  1:120858/1‑249 (MQ=255)                   gAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCCC‑G‑TTCC‑CCCCCCCCTGCTTTTACC‑CTA‑CCCCAGCGGGccg                                                                                                                                                                                                                                   >  1:422415/1‑249 (MQ=255)                   gAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTGTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCa                                                                                                                                                                                                                                    <  1:492406/250‑1 (MQ=255)                    aGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGGTTT‑CCCGC‑G‑CTCC‑GCCCCCCCCGCGTTTACC‑GTT‑CCCACGCggggagg                                                                                                                                                                                                                                  >  2:18779/1‑247 (MQ=255)                        gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAAc                                                                                                                                                                                                                              <  1:215958/250‑1 (MQ=255)                        gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAAc                                                                                                                                                                                                                              <  1:412750/250‑1 (MQ=255)                         ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTT‑CTCGC‑G‑TTCC‑GCCGCCCCTGCTTTTACC‑GTAACCCAgcggggcataaccc                                                                                                                                                                                                                              >  1:632695/1‑238 (MQ=255)                           gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAtcctcc                                                                                                                                                                                                                          <  1:312862/251‑1 (MQ=255)                              gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTT‑CCCGC‑G‑GTCC‑CCCCCCCGTGGTTTTACC‑GGA‑CCCAGGCggggggtaaaccccgcc                                                                                                                                                                                                                        >  1:202967/1‑237 (MQ=255)                              gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CCTGC‑G‑TTCC‑GCCCCCCGTGCTTTAACC‑GTA‑ACCAGGCCGGCA‑GGAACCTccccc                                                                                                                                                                                                                       >  1:414709/1‑248 (MQ=255)                                gTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑g                                                                                                                                                                                                                    <  1:267506/251‑1 (MQ=255)                                    tGCTGGCATTAATATCGGTTTCCACACCTTCAACAGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GAtt                                                                                                                                                                                                                 <  2:160219/250‑1 (MQ=255)                                     gCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTc                                                                                                                                                                                                                <  2:465984/250‑1 (MQ=255)                                        ggCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCttt                                                                                                                                                                                                             <  1:571621/250‑1 (MQ=255)                                              aaTATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑Ga                                                                                                                                                                                                    <  1:465979/251‑1 (MQ=255)                                                        tCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑At                                                                                                                                                                                          <  1:301229/250‑1 (MQ=255)                                                                     cTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTT‑CTCGC‑G‑TTCC‑GCCGCCCGGGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAAACCCCCCCC‑CATTCTTT‑CCCCC‑GATATCGAG‑GTTCCAGccccccctt                                                                                                                                                                            >  1:139551/1‑245 (MQ=255)                                                                       gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTT‑CTCGC‑G‑TTTC‑CCCGCCCCCGCTTTAACC‑CGA‑CCCCCGCGGGCA‑GTAACCCCCCCCC‑CCTTCTTT‑TCCCC‑GATATCCAG‑ATGCCACCccgccggta                                                                                                                                                                           >  1:144251/1‑247 (MQ=255)                                                                         agaagTAACCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGAAAGCCCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGTTTTAACC‑GTC‑CCCAGGCGGGCA‑GTACCCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGc                                                                                                                                                                        <  1:8413/251‑1 (MQ=255)                                                                          gaagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTC‑CCCAGGCGGGCA‑GTACCCCCCGCCC‑GATTCTTT‑CCCCA‑AATATCGAG‑ATGCCAGCCCGCCGTTAAGcc                                                                                                                                                                       <  1:346023/251‑1 (MQ=255)                                                                                       tCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑GTCC‑CCCGCCCCCGCTTTTACC‑GGA‑CCCATGCGGGGA‑GTAAACTCCGACC‑CATTCTTT‑CCCCC‑GATAACGAG‑AAGCCCGCACCCCGCTAAGGC‑AAAAATACC‑gctg                                                                                                                                                        >  1:111521/1‑248 (MQ=255)                                                                                          gATATAGCTGGCGATCGGTCGCGATAATGGAGTTCTGCTCTTTTTCTATTTGATGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGATACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCCC‑G‑TTCC‑GCCCCCCGTGCCTTTACC‑GCA‑CCCCGGGGGGGA‑GGAACCTCCGCCC‑GAGTCTTT‑CCCCA‑GAAATCGAG‑ATGCCAGCCCCCCGGTAAGCCCAAAAgacccttgggg                                                                                                                                                       >  1:500952/1‑241 (MQ=255)                                                                                              taGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑GGTTCC‑‑CCGCCCGGGCTTTTAAC‑GTA‑CCCCAGCGGGCA‑GTAAACCCCCCCC‑GCTTTTTT‑CCCCC‑GATATCGAG‑ATGCCAGCCCGCCCTTAAACC‑CAAAATACC‑CTTGGGGTA‑cc                                                                                                                                                >  1:611498/1‑251 (MQ=255)                                                                                                          ggtcgcattaatgtggtttggctttttttCTAGTTGCGGCAAGACCG‑GTTTAAACGCCGCC‑GGGGTTG‑GCACCCCCTAGTACGGGCGGGCTACGCGGCCAGCCAGCAAT‑GGGGGGTTTG‑‑GCGCTG‑TTCCCGCCCCCCGTTCTTGAACT‑GCC‑CCCAGG‑GGGAACGTATCCCCCGCCC‑GATT‑TTT‑CCCCAAGTTACGGAG‑AGCCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCg                                                                                                                                    <  2:135951/227‑1 (MQ=255)                                                                                                                 aTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCCCCCGGGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTTTCCCCA‑‑AAAACGAG‑ATGCCAGCCCCCCGTTTAGCC‑CAAAGTACA‑GTTGGGGTA‑CCACCGGCCGGGCGTACCTGc                                                                                                                             >  1:136122/1‑251 (MQ=255)                                                                                                                   aaTGGTTTTCTGCTCTTTTTCAATTGGCGGCAAGCCCGTGTTTAC‑CGCCGCCTGGGTTTG‑CCCCCCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTACCCTCCGCCC‑GATTCTTT‑CCCCA‑GATAGCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCgg                                                                                                                           <  1:364547/251‑1 (MQ=255)                                                                                                                    aTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑TTCC‑GCCCCCCGTGGTTTAAAC‑GTA‑CCCAAGCGGGCA‑GTAACCCCCGCCC‑GATTTTTT‑TCCCC‑GATATCGAG‑ATGCCAACCCGCCGTTAAGCC‑AAAAGTAAC‑GTTGGTGTA‑CCCCGGGCCGGTCGTACCCGCGGc                                                                                                                          >  1:556319/1‑251 (MQ=255)                                                                                                                      ggttttttGCTTTTTTTTTATTTCCTGCAAGCCCGTTTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCAACGCTGCCAGCCAGCAATGGGCGGTTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGGGCTTTACCC‑GTA‑CCCAGGCGGGCA‑GTAACCCCCGCCC‑GATTCTTC‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGcc                                                                                                                         <  1:304554/243‑1 (MQ=255)                                                                                                                         gttcgtctctTTTTCTTTTTGCTGAAAGCCCGTGTTTCC‑CGCCTCGTGGTTTTC‑CCACCCCC‑AGTACAGGTGGGTTCCGCTGCCAGCCAGAATTTGGCGGGTTTG‑CTCCC‑G‑TCCC‑GCCCCCCTTGTTTTAACC‑TCA‑CCCAGGCGGGAA‑GTAACCCCCCCCC‑GATTCTTC‑CCCCA‑GAGTTCGAG‑ATGCCCGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGa                                                                                                                      <  2:383311/244‑1 (MQ=255)                                                                                                                           tCTGCTCTTTTTCTATTTGCTGAAAGACCGTTTTTAA‑CGCCGCC‑GGGTTTGGCCCCCCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑CCCGCCCGTGTTTTAACC‑GTA‑CCCAGGCGGGAA‑GTACCCCCCGCCC‑GATTCTTC‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATc                                                                                                                    <  1:441837/250‑1 (MQ=255)                                                                                                                               ttctTTTTCTATTTGCTGAAAGACCTGGTTTAC‑CGCCGCCTGGGTTTG‑CCCCTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CGCGC‑G‑TTCC‑GCCCCCCGTGCTTTAACC‑GTC‑CCCAGGCGGGCA‑GTACCCCCCGCCC‑GATTCTTC‑CCCCA‑GATATGGAG‑ATGCCGGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTc                                                                                                               <  1:402575/249‑1 (MQ=255)                                                                                                                                                        gtgtTTAC‑CCCCGCCTGGGTTTC‑CCCCCCCC‑AGTACAGGCGGGCTACGCTCCCAGCCAGCATTTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCCCCCGTGCTTTACCC‑GTC‑CCCAGGCGGGAA‑GTACCCTCCCCCC‑GATTCTTC‑CCCAA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTc                                                                                     <  1:41384/251‑1 (MQ=255)                                                                                                                                                               a‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTT‑CTCGC‑G‑TTTC‑GCCCCCCGGGCTTTAAAC‑GGA‑ACCAGGCGGGCC‑GGAACCCCCGGCC‑GATTCTTT‑TCCCA‑GATATCGAG‑GAGCCAACCCGCCGGTAAACC‑AAAAGGACCCGTTGGGGAA‑CCA‑GGGCCGGTCGGACCTGCGGCCGGATCG‑GTCAGAGCAAAAGGCCCCATTACACCTTCTGCC‑gcc                                                                             >  1:490508/1‑249 (MQ=255)                                                                                                                                                                  gccgccGGGTTTTG‑CCACCCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGTTTTG‑CGCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTC‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATACCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACAt                                                                           <  1:469150/251‑1 (MQ=255)                                                                                                                                                                   ccgctggggTTTG‑CCTCCCCC‑AGTACGGGCGGGTTACCCCGCCCGCCAGCAATGGGCGGGTTGG‑CTCGC‑G‑TCCC‑GCCCCCCGTGCTTTAACC‑GTC‑CCAAGGCGGGAA‑GTACCCCCCGCCC‑TATTCTTC‑CCCCA‑GATATCGAG‑TGGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACAtt                                                                          <  2:380196/245‑1 (MQ=255)                                                                                                                                                                    ccgccggggttgccctccccc‑AGTACGGGCGGGCTCCGCTCCCGGCCAGCAATTGGCGGTTTTG‑CCCGC‑T‑TCCC‑CCCCCCCGTGTTTTAACC‑GTC‑CCGGGGGGGGCA‑GTACCCCCCCCCC‑AATTCTCT‑CCCAA‑GATACCGAG‑AGGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACAtt                                                                          <  2:233274/233‑1 (MQ=255)                                                                                                                                                                        tGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑GTTC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCCCCGCCC‑GGTTCTTT‑TCCCA‑GATATCCGG‑ATGCCAGCCCGCCGGTAAGCC‑CAAAGTACC‑GTTGGTGGA‑CCCCCGGCCGGGCGGACCTGCGGCCGGGTCG‑GTCAGAACAAAACGCCCCCTTAAGCCTTCCGCC‑GACAttttttt                                                                     >  1:616046/1‑250 (MQ=255)                                                                                                                                                                               g‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCCTGCTTTTACC‑GTA‑CCCCGGCGGGCA‑GGAAACTCCGCCC‑CGTTCTTT‑CCCCC‑GATATCGAGCATTCAACCCCCCCGT‑AAGCC‑AAAAATACC‑GCTGGTGTA‑CCACGGGCCGGGCGTGGCTCCGGCCGGGGCGGTTC‑GCGCAAACGGCCCCCCTAAGCCCTCTCCCCGAC‑TTTTTTGCCGCTcc                                                              >  2:36879/1‑251 (MQ=255)                                                                                                                                                                                         gataacgggCGGGCAACGCTGCCGGCCAGCACTGGGCGGGTTTG‑CCCGC‑G‑TTCC‑GCCCCCCGTGCTTTAACC‑GTA‑CCCGGGCGGGCT‑GTAACCCCCGCCC‑GATTCTTT‑CCCAA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATccgc                                                      <  1:389343/244‑1 (MQ=255)                                                                                                                                                                                          agaacaggcgggttcccgtGCCCGCCAGCATTGGGCGGTTTTC‑CTCGC‑T‑TCCC‑CCCGCCCGTTCTTTAACC‑GTC‑CCCGGGGGGGCA‑TTACCCCCCGCCC‑GTTTTTTC‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATccgcc                                                     <  2:611498/233‑1 (MQ=255)                                                                                                                                                                                          aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTG‑CTCGC‑G‑GTCC‑CCCGCCCGCGCCTTAAAC‑GGA‑ACCCGGCGGGCA‑GTAACCTCCGCCC‑GGGTCTTT‑CCCCC‑GATATCGAG‑CTGCCAGCCCCCCCTTAAAAC‑CAAAGAACC‑CGAGGGGTA‑ACACCGGGCGGGCCTAACCGCGGCCGGATAG‑GTGAGAGCAAACGGCCCCAATAAGCCCTCTGCC‑GAAATTTTTTTCCGcctcctatccggc                                                     >  1:483337/1‑239 (MQ=255)                                                                                                                                                                                                ggggggCTACGCTGCCACCAAGCATTTGGGGGTTTTC‑CTCGC‑G‑TTCC‑GCCCCCCTTGCTTTACCC‑GTC‑CCCAGGCGGGCA‑GTAACCCCCGCCC‑GTTTCTTC‑CCCCA‑GATACCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGAAAAg                                               <  1:18779/248‑1 (MQ=255)                                                                                                                                                                                                 gCGGGCTACGCGGCCAGCCGGCAATGGGGGGGTTTG‑CCCGC‑T‑TCCC‑GCCGCCTGTG‑TTTAACCCGTC‑CCCGGGCGGGAA‑GTACCCCCCCCCA‑GATTCTTC‑CCCCA‑GATACCGAG‑GTGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGAAAAGt                                              <  2:422415/251‑1 (MQ=255)                                                                                                                                                                                                    ggCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑ATCGC‑G‑TTCC‑CCCGCCCGTGCTTTAAAC‑GTA‑CACAGGCGGGCC‑GTAACCTCCGCCC‑CATTCTTT‑CCCCA‑GGTATTGAG‑ATGCCCGCCCGCCGGTAAGCC‑AAAAATACC‑GTTGG‑GTAACCACGGGCCGGTCGACCCTTCGGCCGGCACG‑TTCAGCGCAAACCACACCATTATGCCTTCTGCC‑GACCTTTTTTTGCGCTCCATCTCCACCGACAAGTCg                                            >  2:530722/1‑250 (MQ=255)                                                                                                                                                                                                         ggCTGCCGGCCAGCAAGGGGCGGTTTTC‑CCCGG‑T‑TTCC‑CCCCCCCGTGTTTTACCC‑GTG‑CCCAGGTGGGCG‑GAACCCTCCCCAC‑GATTCTTC‑CCCCA‑GATTTGGGG‑TTCCCCCGCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGTTCG‑TTCAGAGCAAACGGCCCGATTAAGCCTTCTGCC‑GACATGTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGc                                      <  1:504021/250‑1 (MQ=255)                                                                                                                                                                                                                              gggTTTGTCTCGC‑G‑T‑CC‑CCCCCCCGTTCTTTAACCCGT‑‑CCCAGGGGGGCC‑ATAACCCCCCCCCCGTTTCTTC‑CCCAA‑AATATCGAG‑ATCCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAAc                  <  2:139551/251‑1 (MQ=255)                                                                                                                                                                                                                               ggTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACATCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGATAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGACCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACTGACCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCCTCTCCGCCGTAAAGTCGATCTGCGGCCCATCCTGGGTGGCACCTg                >  2:108903/1‑251 (MQ=255)                                                                                                                                                                                                                                              ttCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGCAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCCCCTgtg    >  2:347828/1‑249 (MQ=255)                                                                                                                                                                                                                                                cc‑gccgCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GAACTTTTTTGCCGCTCCCTATCCGCCAAAAAGTCGATCTGCGGACCCTCCTGAGTGGCAACTGACTTCGCCGGGGCt  >  2:53434/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                  TTTCCCCGCGACATCGTTGAGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAA‑CGCCGCCTGGGTTTG‑CCACTCCC‑AGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTG‑CTCGC‑G‑TTCC‑GCCGCCCGTGCTTTAACC‑GTA‑CCCAGGCGGGCA‑GTAACCTCCGCCC‑GATTCTTT‑CCCCA‑GATATCGAG‑ATGCCAGCCCGCCGTTAAGCC‑AAAAGTACC‑GTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCG‑TTCAGAGCAAACGGCCCCATTAAGCCTTCTGCC‑GACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCT  >  NC_000913/2228001‑2228484

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