Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,710	0	T	G	28.8%	94.0 / 10.3	60	V43G (GTG→GGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (11/31);  new base G (16/1);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCAC  >  NC_000913/729467‑729957                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            ataAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGATCGGAAAACCGGCGGGGCGTCTGTGCGACATGACGCAGTATGGCGGTAGCATTGTCGAGGTGTCAGCCGGGGGACGCATCGTTGCTCCCACCGGCGTTGCCTGTT‑CTGTGT‑GCCCCGGCGCGGTGACGTct                                                                                                                                                                                                                                                      <  2:353587/251‑2 (MQ=18) ataAATATTTTCCTAATTGTTCATATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTcc                                                                                                                                                                                                                                                      <  1:217562/251‑1 (MQ=255)       atTTTCATAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCCGGGCCCCCACCGGCGTGGCCTGTT‑AGGGGT‑GGCCCCGCGGGGGGAGGTCCgggccc                                                                                                                                                                                                                                                >  1:419490/1‑247 (MQ=17)                aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAt                                                                                                                                                                                                                                       <  1:334484/251‑1 (MQ=255)                aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCCTCGGCGCCCCCACCCGGGTGGGCCGTT‑CGGGGT‑GGCCCCGCGGGGGGGCGTCCCGCCAACCgggcca                                                                                                                                                                                                                                        >  1:268468/1‑246 (MQ=17)                           cTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACGTCCGGCCCTCCGGTCCATCCCCTGCCCgg                                                                                                                                                                                                                            >  1:182271/1‑251 (MQ=255)                             gACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTc                                                                                                                                                                                                                              <  2:78406/247‑1 (MQ=255)                                   tatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGAGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGGCGTCCCGCCCTCCCGGCAACCCGCTGCTCGGTG‑G‑AAAgg                                                                                                                                                                                                                  >  1:328554/1‑251 (MQ=14)                                                                gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCt                                                                                                                                                                                    <  1:378738/250‑1 (MQ=255)                                                                gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGAGAAAACAAGGATGAGCTGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATGGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCGTCCACCGGCGGGGCCTGTT‑CGGTGT‑GCGCCGGCGGGTTGACGTCCGGCCATCCGGTCAAGCCGCGGCGCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTg                                                                                                                                                                                   <  1:522074/251‑1 (MQ=14)                                                                        ccGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACGGCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACTTCGCCCTgcccggccc                                                                                                                                                                           <  2:501417/250‑1 (MQ=255)                                                                              gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCGTTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGG‑GCCCCCGCGGGGGGACGTCCGGCCCTCCCGTCAAACCCCTGCTCGGGG‑C‑AAAGG‑GCC‑TTTCCCGTGAAACCGACATCCCCCCtcccggccccccggc                                                                                                                                                                     >  1:11833/1‑245 (MQ=14)                                                                              gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGcc                                                                                                                                                                     <  1:305239/251‑1 (MQ=255)                                                                                cagcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGG‑GGCCCCGCGGGGGGACGGCCGGCCCTCCCGTCAAACCGCTTCCCGGGG‑CCAAAGG‑‑CC‑CTCCCGGGGGAACCCACATCCCCCCgcccggccccccgcccg                                                                                                                                                                   >  1:197477/1‑243 (MQ=9)                                                                                   cagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGGGGCCTGTT‑CGGTGT‑GCCCCCGCGGGGGGACGGCCCGCCACCCGGGCAATCCCCTGCTCGGGG‑C‑AAAGG‑GCC‑TTCCCGGGGAAACCGACAACGCCCTGCCCGGCCCGCCGCCGTTcc                                                                                                                                                                >  1:147694/1‑250 (MQ=14)                                                                                        aaaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAAACGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGATCAACCGGGGTACGCATCGGGGAACCCCCCGGCGTGGCCTGTT‑GCGTGTTGCCCCGGGGGGGGGA‑GTCCGGGCACACGGTCAACTCGCCGCCCCGCG‑C‑AAAGG‑TCCCTTCCCGG‑GAAACAGACCTCCCCCCccccggccccccgcccgttcttctt                                                                                                                                                           >  1:361038/1‑235 (MQ=9)                                                                                         aaTAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATcc                                                                                                                                                                                                                                              <  1:287144/171‑1 (MQ=255)                                                                                           tAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGGGT‑GCCCCCGCGGGGGGGCGTCCGGCCCTCCGGTCAATCCCCTGCTCGGGG‑C‑AAAGG‑GCC‑CTCCCGGGGAAAACGGCATCGCCCCGCCCCGGCCGCCGCCGGTCCTTCTCTcc                                                                                                                                                        >  1:472012/1‑251 (MQ=14)                                                                                                       aCAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCCTTGTCCAGGGTTAAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGCGT‑GGCCCCGCGGGGGGTCGGCCGGCCATCCGGTCACTCCGCTGCTCCGTGGC‑AAGGG‑CCT‑TTCCCGG‑GGAACCCCCCTCGCCCTGCCCCGCCCGCCGCCCGTCCTTCTCTCCCCCACCCAcagc                                                                                                                                            >  1:235741/1‑251 (MQ=9)                                                                                                        cAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACGTCCCGCCATCCGGTCCATCCCCTGCTCGGGG‑C‑AAAGGGCCC‑TCCCCGGGGAACCCG‑CATCGCCCCTCCCCGCCCCCCGCCGCTCATTCTCTCCCCCACCCAcagcc                                                                                                                                           >  1:426890/1‑250 (MQ=14)                                                                                                                    gATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGGGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACCTCCGGCCCTCCGGGCCATCCCCTGCTCGGGG‑C‑AAAGGGCCC‑TTCCCGG‑GAAAACGCCATCCCCCCGCCCCGCCCGCCGCCGTTCAATTTCTccccccccccaaacagttacccggca                                                                                                                                >  1:453091/1‑227 (MQ=2)                                                                                                                        aGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAAc                                                                                                                           <  1:264402/251‑1 (MQ=255)                                                                                                                         gCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                                                                          <  2:459976/251‑1 (MQ=255)                                                                                                                            gAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAAGCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGccc                                                                                                                       <  1:244258/251‑1 (MQ=255)                                                                                                                             aaaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcg                                                                                                                      <  1:103412/251‑1 (MQ=255)                                                                                                                             aaaaCCGGCGGCGCGTCAGGGCGACAGGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCTCCTACAGCAGTTACCGGACAAAAACGCCcg                                                                                                                      <  1:288213/251‑1 (MQ=255)                                                                                                                              aaaCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcgc                                                                                                                     <  1:152333/251‑1 (MQ=255)                                                                                                                                   ggcggcGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TGCCCGGTGAAACCGACAGCGCCCTGCCCTGGCCGCTGCCGTGCAGTCTCGCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCgg                                                                                                                <  1:217321/251‑1 (MQ=14)                                                                                                                                   ggcggcGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACGGCCGGGCCTCCGGGCAATCCCCTGCTCGGGG‑C‑AAAGG‑GCC‑CTCCCCGGGAAACCGACATCGCCCTGCCCGGCCCCCCGCCGGTCATTCTCCCCCCCCCCCACAGCCGTTTCCGGGCAAAAACGCCCCCGCCgg                                                                                                                >  1:435180/1‑251 (MQ=14)                                                                                                                                          cgTCAGGGCGACACGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAg                                                                                                         <  1:426970/251‑1 (MQ=255)                                                                                                                                            gcAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGcc                                                                                                       <  2:183568/250‑1 (MQ=40)                                                                                                                                             cAGGGCGCCATGGCGCAGTATGGCGGTAGCATTGGCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCt                                                                                                      <  2:478676/251‑1 (MQ=255)                                                                                                                                               gggCGAAATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCg                                                                                                    <  2:58/251‑1 (MQ=255)                                                                                                                                                 gCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGGGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCGCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCGTCGGc                                                                                                  <  2:472012/251‑1 (MQ=255)                                                                                                                                                   gACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑TGGTGT‑GCCCCCGCGGGGGGACGTCCCGCCCTCCGGTCAATCCGCCGCTCCGGG‑C‑AAAGG‑GCC‑CTCCCCGGGAAACCCACATCGCCCCGCCCCGCCCCCCTCCGTTCCTTCTCTCCCCCCCCTCCAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGGGCCTCCGccc                                                                                                >  1:284225/1‑251 (MQ=9)                                                                                                                                                        gACGCAGTATGGCGGTAGCATTGGCCAGGGTGCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCGCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCt                                                                                           <  1:520348/251‑1 (MQ=255)                                                                                                                                                           gCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGGCGTCCCGCCCTCCGGTCAATCCGCCGCCCGGTG‑C‑AAAGG‑GCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTTTCCCGCACCTACAGCAGCTACCGGACAAAAAAGCCCGCGCCGGGGGGGGGCCTCGGGCCCGGGTggg                                                                                        >  1:20721/1‑250 (MQ=14)                                                                                                                                                           gCAGGATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGa                                                                                        <  2:147694/251‑1 (MQ=255)                                                                                                                                                                      gTAGCATTGTCCAGGGTTCAGCCGGGGAACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGAGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGGG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGAACCTACAGCCGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCg                                                                             >  1:189794/1‑251 (MQ=255)                                                                                                                                                                       tAGCATTGTCCAGGGTTCAGCCGGGGTACGAATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACCGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCCGGAAACTGCCTGCgg                                                                            >  2:245710/1‑251 (MQ=255)                                                                                                                                                                         gCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGGT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGat                                                                          <  2:138777/251‑1 (MQ=255)                                                                                                                                                                           aTTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCCCTCCCGCACCTACAGCAGTTACCGGCCAAAAACGCCCCCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCCGCGGatat                                                                        >  2:102165/1‑251 (MQ=255)                                                                                                                                                                                 cAGGGTTCAGCCGGGGAACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGGGACGTCCGGCCCTCCGGGCAATCCGCTGCTCCGTG‑C‑AAAGG‑TCC‑TTCCCGGCGAAACCGACATCCCCCCGCCCGGCCCGCCGCCGTTCATTCTCTCCCCCCCCCTCAGCAGTTACCGGACAAAAACGCCCGGGCCGGTGGGGAGCCCCGGCCCCGGCCGGaacaagaccgcgggtatcccatt                                                                  >  2:191465/1‑230 (MQ=14)                                                                                                                                                                                  ggggtccggCCGGGGTACGCGGCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGGCGGCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCGGCCGTTCATTCTCTCCCGCGCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATAGCCGCTTa                                                                 <  2:197477/243‑1 (MQ=11)                                                                                                                                                                                             ggggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGcgcg                                                       >  1:150091/1‑250 (MQ=39)                                                                                                                                                                                                   cgaaccGGTGCCCCCACCGGCGTGGGC‑GTTGCGGTGG‑GCCCCGGCGGGGTGACGGCCGGCCACCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGGGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATGCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGGGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacac                                                <  2:235741/246‑1 (MQ=11)                                                                                                                                                                                                   cGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGGGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGCTAacac                                                >  1:159602/1‑251 (MQ=37)                                                                                                                                                                                                           tGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGGCAATCCCCTGCTCGGTG‑C‑AAAGG‑GCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGACGTTACCGGACAAACACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCCGCGGATATCCGCTTACAGCTGCGCGAAAACACACTGaaac                                        >  1:345462/1‑248 (MQ=11)                                                                                                                                                                                                             cccccACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTc                                      >  1:9606/1‑251 (MQ=34)                                                                                                                                                                                                               cccACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTc                                      >  1:122540/1‑249 (MQ=35)                                                                                                                                                                                                                 cccctacacgacgcgcttccgatcggCCCCGGCGGGGTGACGGCCGGCCATCCGGTCAATCCGCGGCGCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACACCGCCCTGCCCGGCCCGCTGCCGTTCATTCGCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCTGCCCCGGCTTGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAg                                    <  2:235804/226‑1 (MQ=11)                                                                                                                                                                                                                   ccGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGAt                                <  2:182271/251‑1 (MQ=34)                                                                                                                                                                                                                     ggCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGGTATCCGCTTACAGCTGCGCGATAACACACTGATACTCCGTGATa                               >  2:158179/1‑250 (MQ=25)                                                                                                                                                                                                                     ggCGTGGCCTGTT‑CGGTGG‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATaa                              <  2:11833/251‑1 (MQ=21)                                                                                                                                                                                                                         tGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGGCGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGCAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGCCAACAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGGCTGGGGGTATCCGATTACAGCTGCGCGACACCACACTGTTACTCAGTGATAAcggc                          >  2:146638/1‑251 (MQ=11)                                                                                                                                                                                                                          ggCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggc                          <  1:389425/250‑1 (MQ=35)                                                                                                                                                                                                                               gTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCa                      <  1:158179/249‑1 (MQ=35)                                                                                                                                                                                                                                          gCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCAGCTGCCGTTCATTCTCTCCCGCACCTACAGCCGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAg                                    >  1:235804/1‑226 (MQ=25)                                                                                                                                                                                                                                              cggcggGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTg       <  1:131192/251‑1 (MQ=34)                                                                                                                                                                                                                                                   gggTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑ACC‑TTCCCGGTGAAACCGACATCGCCCTAACCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCAGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTTAGCAc  >  1:346167/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                            ATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTT‑CGGTGT‑GCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTG‑C‑AAAGG‑TCC‑TTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTTTGAGCAC  >  NC_000913/729467‑729957

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