Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,710	0	T	G	43.5%	37.7 / 12.2	48	V43G (GTG→GGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (8/18);  new base G (20/0);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 5.65e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGC  >  NC_000913/729465‑729955                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           atataAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTTAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCCGGTCGGGGTGCCCCCG‑CGGGGGGACgg                                                                                                                                                                                                                                                       >  1:459785/1‑250 (MQ=18)       atatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGGTCAGCCGGGGGACGCCTCGGGGCCCCCACCGGCGTGGCCCGTTCGGTGTGCCCCGG‑CGGGGGGGCGTCCGGcc                                                                                                                                                                                                                                                 >  2:220696/1‑251 (MQ=18)                     gTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGG‑CGGGGGGACGTCCGGCCATCCGGTC‑A‑ATccc                                                                                                                                                                                                                                 >  1:236419/1‑250 (MQ=255)                             cTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGACCCCACCAGCGTGGCCTGTTAGGTGTGCCCCGG‑CGGGGGGGCGTTCGGCCATCCGGGC‑A‑Atacggcgccggg                                                                                                                                                                                                                         >  1:177254/1‑240 (MQ=17)                                          cACTTACGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGAAGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACATTGTCCAGGGGTCAGGCGGGGTACCCCTAGGTGCCCCCCCCCGCGTGGCCGATTCGGGGTGCCCCGT‑CGGGGGGGCGTCCGGCCAACCGG‑C‑ACATCCGCTGTGCGGTG‑C‑AAAGGGCCtt                                                                                                                                                                                                          >  2:488042/1‑251 (MQ=14)                                          cACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCGCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCtt                                                                                                                                                                                                          <  1:4459/251‑1 (MQ=255)                                                     ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGAAGTATGGCGGGAGCATTGTCCAGGGGTCAGCCGGGGTATGCACAGGTGCCCCCACCGGCGGGGGCTGTTCGGTGCGCCCCGG‑CGGGGGGACGTCCGGCCATCCGGGC‑A‑AACCGCCGCGCGGGG‑C‑AAAAGGCCTTTCCggcgtaac                                                                                                                                                                                               >  1:41391/1‑245 (MQ=14)                                                                           cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTtccc‑ggccc                                                                                                                                                                         <  1:170911/250‑1 (MQ=255)                                                                           cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGGGCCCCGG‑CGGGGGGGCGTCCGGCCCTCCGGTC‑A‑AACCCCTGCTCGGGG‑C‑AAAGGGCCTTCCCGGGGAAAACGACATCGcccccccccggccc                                                                                                                                                                         >  1:14801/1‑240 (MQ=14)                                                                             ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCt                                                                                                                                                                      <  1:499243/251‑1 (MQ=255)                                                                                     cagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGGCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCGGGCCTCCCGTC‑A‑ATCCGCCGCTCGGGGGC‑AAAGG‑CCTTCCCGGTGGAACCGCAAACGCCCCGCCC‑GGCCCCCCGCCGtttt                                                                                                                                                              >  1:451309/1‑249 (MQ=14)                                                                                       gACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTAAGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCCGCCCTCCGGTC‑C‑ATCCGCCTCTCGGTG‑C‑CAAGGGCCCTCCCGGGGGAACCGACATCGCCCTGCCC‑CGCCCCCTGcccgtcctt                                                                                                                                                            >  1:317879/1‑244 (MQ=14)                                                                                          aaaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATtctc                                                                                                                                                         <  1:405734/251‑1 (MQ=255)                                                                                          aaaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGTCCGCCGGCCAACCGGGC‑C‑ATCCGGCGGTTGGTG‑G‑AAAGGGCCTTTCCGGGGGAACCGACAGCCCCCCGCCC‑GGCCCCCCGCCGTTCAtttcc                                                                                                                                                         >  2:402487/1‑248 (MQ=14)                                                                                           aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCCGCGTGGCCTGTTCGGCGTGCCCCGG‑CGGGGGGAAGGCCCGCCATCCGGGC‑C‑ATCCGCTGCTCCGGG‑C‑ACAGGGCCCTCCCGGGGAAACCGGCAACGCCCCGccccggcccccgcccttcctttccc                                                                                                                                                         >  2:458019/1‑235 (MQ=14)                                                                                               aCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCGACGGTGCCCCCACCGGCGGGGGCTGTTTGGTGTGCCCCGGGCGGGG‑GACGGCCGGCCCACCGGGC‑C‑ATCCGCTGCGCGGGG‑C‑AAAGGGCCTTCCCGGAGGAACCGACATCGCCCTgccc‑gggcccccgccgcttcttccctcccc                                                                                                                                                    >  2:110936/1‑227 (MQ=14)                                                                                               aCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCg                                                                                                                                                    <  1:209299/251‑1 (MQ=255)                                                                                                cccGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCCG‑CGGGGGGACGGCCGGCCATCCCGGC‑A‑ATCCCCTGCTCGGGG‑C‑AAAGGTCCTTCCCGGGGGAAACGACATCCCCCCGCCC‑CGCCCCCCGCCGCTTATTTTCTccccc                                                                                                                                                   >  1:257929/1‑249 (MQ=14)                                                                                                   gATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCAcc                                                                                                                                                <  2:5829/251‑1 (MQ=255)                                                                                                       aaaCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGG‑CGGGGGGGAGTCCGGCCATTCGGTC‑A‑ATCCGCTGCTCGGGG‑C‑AAAGGTCCTTTCCGGTGAAACCGCCATCGCCCTGCCC‑GGCCCGCTCCCGTTCATTTTCTTCCCCAcccac                                                                                                                                             >  1:486838/1‑247 (MQ=14)                                                                                                        aaCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCCGCCATCCGGTC‑A‑ATCCGCTGCTCGGGG‑C‑AAAGGGCCTTCCCGGTGAAACCCACATCGCCCCTCCC‑GGCCCCCTGCCGTTCCTTCTCTCCCGCACCcacac                                                                                                                                           >  1:377696/1‑250 (MQ=14)                                                                                                                 gagcgaaTGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATtctctc                                                                                                                                                       <  1:146369/224‑1 (MQ=255)                                                                                                                                  aCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGGCCCGG‑CGGGGGGGCGTCCGGCCCTCCGGTCCA‑ATCC‑CTGCTCGGGG‑C‑CAAGGGCCCTTCCGGGGAAACCGGCATCGCCCCGCCC‑CGCCCCCTGGCGTTCATTTTCTCCCCCACCCACCGCAGGTACCCGACAAAAACGCCCccgc                                                                                                                 >  1:268411/1‑251 (MQ=9)                                                                                                                                  aCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCCGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGGCCTTCCCGGGGAAACCGACATCGCCCTGCCC‑GGCCCGCCGCCGTTCATTTTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCgcgc                                                                                                                 >  1:32711/1‑251 (MQ=14)                                                                                                                                     ggcggcGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTCGGTGGGCCCCGG‑CGGGGGGGCGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGGG‑C‑AAAAGTCCTTCCCGGGGAAACCGCCATCGCCCCGCCC‑GGCCCGCCGCCGGTTATTCTCTCCCGCCCCTACAGCCGTTACCGGGCCAAAACCCCCGCCCCgg                                                                                                              >  1:163941/1‑251 (MQ=14)                                                                                                                                           gcgTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGGGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCGG‑CGGGGGGACGGCCGGCCCTCCGGGC‑A‑ATCCGCTGCTCGGTG‑C‑AAAAGGCCTTCCCGGGGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTTATTCTCTCCCGCACCTACAGCAGGTACCGGACAAAAAAGCCCCCGCCGggggggg                                                                                                        >  1:147122/1‑250 (MQ=14)                                                                                                                                                    cGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGATCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGACCGGCCCTCCCGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGGCCTTCCCGGTGAAACCGACATAGCCCTTCCC‑GGCCCCCTGCCCGTCATTCTCTCCCGCGCCCACAGCAGTTACCGGCCAAAAACGCCCGCGCCGGAgggggggcccggcc                                                                                               >  1:267290/1‑241 (MQ=14)                                                                                                                                                     gACATGACGCAGTATGGCGGTAGCCTTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGACCCGG‑GGGGGTGACGTCCGGCCCTCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGGACTTCCCGGCGAAAACCACATCGCCCTGCCC‑CGCCCGCCGCCGGTCATTCTCTCCCGCCCCTCCCGCAGTTACCGGACAAAAAAACCCGCGCCGGTGGGGGGCCTCGGccc                                                                                              >  1:509034/1‑251 (MQ=255)                                                                                                                                                                               tGTCCAGGGTTCAGACGGGGTACGCATCGGTGCCCCCACCGGCGTCACCTGTTCGGTGTGCCCCGG‑CGAGAACAAGACCGGCCCTACGATG‑ACACCCG‑TGTTCGGTC‑CCACGGGACCTTCC‑GGTGACACCGACACCGCCCTGACC‑CGCCCGCACACGTTCATTCTCTTCCGCCCCTACAGCAGTAACCGCACACAAACGCACGCGCCGGTTGCGCGCCTCACGCCCGGCTCGAACATGCCTgtgggtgtcc                                                                    >  1:262141/1‑242 (MQ=255)                                                                                                                                                                                gTCCAGGGTTCAGACGGGGTACGCATCGGTGCCCACACCGGCGTGGCCTGTTCGGTGTGCCTCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑CTCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTAACGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGCAAATGCCTGCGGCTATCCg                                                                   >  2:355292/1‑251 (MQ=14)                                                                                                                                                                                          cAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGGG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCCGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCCCGGCCCCGGCCGGAAAATGCCTGCGGGTATCCGCTTACAGCtc                                                         >  2:248457/1‑250 (MQ=255)                                                                                                                                                                                            gCCGGGGCACGCGTCGGTGCCCCCACCGGCGTGGGCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                                       <  2:114437/251‑1 (MQ=255)                                                                                                                                                                                               ggggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGGGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGGAAATGCCTGCGGGTATCCGCTTACAGCTGcgcg                                                     >  1:301717/1‑250 (MQ=255)                                                                                                                                                                                                     cGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacac                                              <  1:220696/251‑1 (MQ=37)                                                                                                                                                                                                       cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCGGCGCGATAacacac                                            <  1:419766/251‑1 (MQ=37)                                                                                                                                                                                                       cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGACCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGCCAAAAACGCCCGCGCCGGTGGGGGGCCTCGGCCCCGGCTGGAACAAGCCTGCGGATATCCGCTTACACCCCCGCGATAacaaac                                            >  2:236806/1‑248 (MQ=11)                                                                                                                                                                                                          cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTg                                          <  1:9499/250‑1 (MQ=35)                                                                                                                                                                                                           ggTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcgcg                                                                                                                  <  2:162600/177‑1 (MQ=38)                                                                                                                                                                                                           ggTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcgcg                                                                                                                  >  1:162600/1‑177 (MQ=38)                                                                                                                                                                                                           ggTGCCCCCACCGGCGTGGCCTGTGCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAt                                        <  2:14801/251‑1 (MQ=21)                                                                                                                                                                                                                          gTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACCCTGGTACTCAGTGCTAAcgg                         >  2:327445/1‑251 (MQ=21)                                                                                                                                                                                                                                ggggcGGTGTGCCCCGG‑CGGGTTGTCGCCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAGGGTCCTGCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGCTCATTCTCGCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCTCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAGACTCAGTGATAACGGCGGCAg                   <  2:108733/247‑1 (MQ=11)                                                                                                                                                                                                                                 gTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGa                  <  1:58219/251‑1 (MQ=34)                                                                                                                                                                                                                                 gTTCGGTGTGCCCCGG‑CGGGGTGACGGCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCa                    <  1:211026/249‑1 (MQ=25)                                                                                                                                                                                                                                  ttCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑AACCGCTGCTCGGTG‑C‑AAAGGGCCGTCCCGGTGAAACCGACAGCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg                                                            <  1:143090/208‑1 (MQ=255)                                                                                                                                                                                                                                       ggtgCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTg            <  2:236419/250‑1 (MQ=21)                                                                                                                                                                                                                                                  cggGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGc  >  2:276311/1‑251 (MQ=34)                                                                                                                                                                                                                                                   ggggTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCCGAAGCCTGGATTTTGAGc  >  2:234431/1‑250 (MQ=21)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           ATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGG‑CGGGGTGACGTCCGGCCATCCGGTC‑A‑ATCCGCTGCTCGGTG‑C‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCC‑GGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGC  >  NC_000913/729465‑729955

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