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

ATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAA‑GCCTGTATTTTGAGC  >  NC_000913/729473‑729955                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             atTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTTGGTGTGC‑CCCGGCGGGGTGGCGGCCGGCCAt                                                                                                                                                                                                                                                 >  1:279717/1‑251 (MQ=255)                  tATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGGGCCCCCACCGGCGTGGCCCGTTCGGTGTGG‑CCCGGCGGGGGGACGGCCGGCCCTcccggtaaacccctgcc                                                                                                                                                                                                                                >  1:360268/1‑236 (MQ=17)                      tGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGACGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGGGACGTCCGGCCATCCGGGCAAT‑CCGCTGCTCggg                                                                                                                                                                                                                           >  1:563820/1‑250 (MQ=255)                                  cACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAACATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCCGCGTGGGCCGTTCGGGGTGC‑CCCGGCGGGGGGCAGGCCGGCCATCCGGGCCAA‑CCGCTGCTccgggacaagggcctt                                                                                                                                                                                                               >  2:379483/1‑236 (MQ=14)                                                     aaTCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGGGACGTCCCGCCCTCCCGGCCAT‑CCCCTGCTCGGTGC‑AAAGGTCCTTCCCCGTGGAACCCACATcc                                                                                                                                                                                           >  1:252979/1‑250 (MQ=14)                                                      aTCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGAATTGTCCAGGGTTCCGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGC‑CCCGGCGGGGTGACGTCCGGGCATCCGGTCAAT‑CCGCTGCTCGGGGG‑AAAAGGCCTTCCCGGTGAAACCCACACCGc                                                                                                                                                                                          >  1:443772/1‑251 (MQ=14)                                                          gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGC‑CCCCGCGGGGGGGCGTTCGGCCATCCGGTCAAA‑CCGGTGCTCGGGGCCAAAGG‑CCTTTCCGGTGGAACCGACATCGCCCTg                                                                                                                                                                                      >  1:259402/1‑251 (MQ=14)                                                                          cagcagACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGCATCGGGGCCCCCACCGGCGTGGGCCGGTCGGGGTGC‑CCCGGCGGGGTGACGTCCGGCCCTCCGGTCAAT‑CCCCTGCTGGGGGG‑CAAGGTCCTTCCCCGTGCAACCCACATccacccgccccgcccaccgcccg                                                                                                                                                                      >  2:29025/1‑229 (MQ=14)                                                                                     tAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTTGGTGTGC‑CCCCGCGGGGGGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTTGGGGG‑AAAGGGCCTTCCCGGGGAAACCGACCTCGCCCTGCCCGGCCCCCTGCCC‑TTCATTCTCTcc                                                                                                                                                          >  1:422421/1‑251 (MQ=14)                                                                                     tAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAACATTGTCCAGGGGTCATCCGGGGTACGCATCGGGGCCCCCACCGGCGTGGGCTGTTCGGGGTGC‑CCCGGCGGGGGGGAGGCCGGCCCTCCCGTACAT‑CCGCCGCCCGGGGC‑GAAGGTCCTTCCCGGTGGAACCGACCTCGCCCCGCCCCGCCCGCGGCCG‑TGCATtctcccc                                                                                                                                                          >  1:547193/1‑248 (MQ=9)                                                                                          cGATAAAACACGGATGAGCAGATGAGCGGAAAACAGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCCGCGTGGCCCGTTCGGTGTGC‑CCCGGCGGGGGGGCGTCCGGCCCTCCGGTCAAT‑CCCCTGCCCCGGGC‑CAAGGGCCTTCCCCGGGAAACCCACATCGCCCCGCCCCGCCCCCCCCCG‑GTCATTCTCTCCCCCAc                                                                                                                                                     >  1:228727/1‑251 (MQ=14)                                                                                                      gattaGCAGATGAGGGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACTCAGTATGGCGGTAGCATGGTCCAGGGTCCAGCCGGGGTGCGCATCGGTGCCCCCGCCGGCGTGGCCTGGTCGGTGTGC‑CCCGGCGGGGTTACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTa                                                                                                                                        <  2:52469/247‑1 (MQ=14)                                                                                                      gATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGGTCAGCCGGGGGACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGGGGGGGGACGTCCGGCCATCCGGTCAAT‑CAGCTGCTCGGGGC‑AAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCCGCCCGCTGCCG‑TTCATTCTCTCCCCCACCTAcagcaggt                                                                                                                                          >  1:220117/1‑248 (MQ=14)                                                                                                         gAGCAGATGAGCGGAAAACCGGCGGAGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAACATTTTCCCGGGGTCAGCCGGGGGGCGCAGCGGTGCCCCCCCCGGGGTGGCCCGTTCGGTGCGCACCCGG‑GGGGTGGAGTCCGGGCAACCGCACCATACCG‑TGCTCCGGGC‑AACGGGCCTTCCCGGGGGAACCGACATACCACTGCGCGGCCCGCCGGCG‑CTCACGCTCTCCCGAACCTaacgccgggtccc                                                                                                                                      >  2:215005/1‑239 (MQ=2)                                                                                                             aGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGATCAGCCGGGGGAAGCACCGGTGCCCCCACCGGCGTGGCCTGTTCCGTGTGG‑CCCGGCGGGGGGACGTTCGGCCCTCCGGGCAAT‑CCGCCGCGCGGGGC‑AAAAGTCCTCCCCGGGGAAACCGACATCGCCCCGCCCCGCCCCCCGCCG‑CTCATTCTCTCCCCCACCAACCGCCG‑TTACCGGAc                                                                                                                                  >  1:406915/1‑250 (MQ=9)                                                                                                                tGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGGGCCCCCACAGGCGTGGCCTGTTTGGTGTGC‑CCCGGGGGGGGGACGTCAGGCCATCCGGGAAAT‑CCGCTGGTCGGAGC‑AAAGGGCCTTCCCGGTGAAACCGGCATCGCCCTGCCCCGCCCGCTGCCG‑CTCATTCTTTCCCGCCCCTACCGCAG‑TTCCCGGACaaa                                                                                                                               >  1:95117/1‑250 (MQ=14)                                                                                                                       aaaaCCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGCAAGCATCGGCGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGGCATCCGGTCAAT‑ACCCTGCTCGGTGC‑AAAGGGCCTTCCCCGTGAAAACGACATCCCCCTGCCCGGCCCCCCGCCC‑TTCATTCTCTCCCGCACCTACAACAG‑TTACCGGGCAAAAAAGcccc                                                                                                                       >  1:328759/1‑250 (MQ=14)                                                                                                                            cggcggCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGTCCCGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCg                                                                                                                  <  2:387527/251‑1 (MQ=255)                                                                                                                             ggcggcGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAACAATGTCCCGGGTTCCGGCGGGGGACGCATCGGTGCCCCCACCGGCGGGGCCTGTTCGGAGGGC‑CCCGGCGGGGGGACGGCCGGCCATCCGGTCAAA‑CCGCTGCTCGGGGC‑AAAGGGCCTTCCCGGTGGAACCGACATCGCCCCGCCCCGCCCGCCGCCG‑TTCCTTCTCTCCCGCACCTCCCACAG‑GTAACGGACAAAAACGCCCCCGCCgg                                                                                                                 >  2:261568/1‑251 (MQ=9)                                                                                                                                ggcgCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGGGGACGCATCGGTGCCCCCAACGGCGGGGCCTGTTCGGTGGGC‑CCCGGCGGGGGGAAGTGCGGCCATCCGGTCAAT‑CCGCTGACCGGGGC‑CAAGGTCCTCCCCCGGGAAACCGACCTCGCCCCTCCCCGCCCCCCGCCC‑TTCATTTTCTCCCCCACCCACAGCAGGTTCCCGG‑CAAAAACGCCCGCGCCCGTgg                                                                                                              >  2:410915/1‑251 (MQ=9)                                                                                                                                  cgcgTCAGGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGC‑CCCGGCGGGGGGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAAGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCCCCGCCG‑TTCATTCTCTCCCGCCCCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCgggggg                                                                                                             >  1:74902/1‑251 (MQ=14)                                                                                                                                    cgTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAACCGGGGTACGCATCGGTGCCCCCACCGGCGGGGGCTGTTCGGTGTGC‑CCCGGCGGGGTGGCGGTCGGCCCTTCGGTCCAT‑CCCCTGCTCGGTGG‑AAAGGGCCTTTCCGGTGAAACCCACATCGCCCCGGCCGGCCCGCCTCCG‑TTCATTTTCTTCCGCACCTACAGCAG‑TTACCGGAAAAAAACCCCCCCGCCGgggggg                                                                                                            >  1:704374/1‑249 (MQ=11)                                                                                                                                         ggg‑CGACAGGACGCAGTATGTCTGTAGCATTGTCCAGTGTTCAGCCGGTGTACGCATCGTTGCCCCCACCGGCGGGGCCTGGTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCGGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑g                                                                                                    <  2:259402/251‑1 (MQ=14)                                                                                                                                               aCATGACGCAGTATGGCGGATGAATTGTCCCGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGC‑CCCGGAGGGGGGACCTCCGGCCATCCGGGCAAT‑CCCCCGCCCGGTGA‑AAAGGGCCTTTCCGGGGAAAACGCCATCGCCCCGCCCGGCCCGCCTCCG‑TTCCTTATCTCCCGCACACCCAGCAG‑TTACCGGGCAAAAAAGCCCCCGCCGGGGGGGGGCCTCCGGccc                                                                                                >  2:88227/1‑251 (MQ=255)                                                                                                                                                 atgactcggTATGGCGGTAGCAGTTTCCAGGGTTCAGCCGGGGTACGCATCGGGGGCCCCACCGGCGTGGCCTGTTCGGTGTGT‑CCCGGCGGGGTGACGTCCGGCCAGCCGGTCAAT‑CCGCTGCTCGGTGG‑AAAGGTCCTTCCCGGTGAAGCCGACAGCGCGCTGCCTGGCCCTCTGCCG‑TTTATTCTCTCCCGCACCGACGGTAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCgg                                                                                             <  1:356485/243‑1 (MQ=14)                                                                                                                                                   gACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTTCGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGTGAGCCTC‑GGCCCCGGCt                                                                                           <  1:555213/251‑1 (MQ=255)                                                                                                                                                   gACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTCCGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCt                                                                                           <  2:175047/251‑1 (MQ=38)                                                                                                                                                     cGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGGGACGTCCGGCCCTCCGGTCCAT‑CCGCTGCTCGGTGC‑AAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCCGCCCCCTGCCGGTTC‑TTCTCTCCCGCACCCACAGCCG‑TTTCCGGACAAAAACCCCCGCGCCGGGGGGGAGCCTC‑GGCCCCCGCTgg                                                                                         >  1:269235/1‑251 (MQ=255)                                                                                                                                                             ggCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGGGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCCTTCTCTCCCCCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGcc                                                                                 >  2:555865/1‑251 (MQ=255)                                                                                                                                                                         gTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCg                                                                     <  2:55285/251‑1 (MQ=255)                                                                                                                                                                           ccAGGGTTCAGCCGGGGTACGCATAGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTTGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAACATGCCTGCGGATATCCGc                                                                    <  2:126232/250‑1 (MQ=255)                                                                                                                                                                             aGGGTTCAGCCGGGGCCCGCAACGGTGACCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTTCGGCACTCAGGTCAAT‑CCGATGCTCGGGGC‑AAAGGTCCTTCCCGGTGAAAACGACCTCCCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAACCCCCCGCCCCGGTCGGGAGCCCC‑GGCCCCCGCTGGGAAATGGCTGCCGATATcccatt                                                                  >  2:348546/1‑246 (MQ=11)                                                                                                                                                                                   cAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTg                                                          <  1:379483/251‑1 (MQ=38)                                                                                                                                                                                        ggggTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCcgcgcg                                                      >  1:94415/1‑250 (MQ=39)                                                                                                                                                                                        ggggTACGCATCGGGGCCCCCACCGGCGTGGCCTGTTAGGTGTGC‑CCAGGCGGGGGGACGTCCGGCAATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGACCTTCCCGGTGAAACCGACATCGCACTGCCCGGCCCGCTGCCG‑TTCCTTCTCTCCCGCACCTACCGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGCGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTCCAGCTGcgcg                                                      >  2:117349/1‑250 (MQ=14)                                                                                                                                                                                         gggTACGCATCAGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATACGGTCAAT‑CCGCTGCCCGGCGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTCCAGCAG‑TTACCGGACACAAACGCCCCGGCCGGTGGGGAGCCCC‑CGCCCCGGCAGGAAGATGCCTGCGGATATCCGC‑TTCCCGCTGCGCGAt                                                    >  1:69785/1‑251 (MQ=14)                                                                                                                                                                                          gggaCGCAGCGGTGCCCCCACCGGCGTGGCGTGTGCGGTGTGC‑CCGGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCGCCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCGGCGGATATCCGC‑TTACAGCTGCGCGATa                                                   <  2:228727/248‑1 (MQ=11)                                                                                                                                                                                           gTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGCCATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGAGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGAACTCCGC‑TTACAGCTGCGCGATac                                                  >  2:386782/1‑250 (MQ=255)                                                                                                                                                                                               gCATCG‑TGCCCGCACCGGCGTGGCCTGTTCGGTGTGC‑CGCGGCGGGGGGACGTCCGGCCATCCGGTCGAT‑GCGCTGCGCGGTGC‑AAAGGGCGTTCCCGGTGAAGCCGACAGCGCCCTGCCCGGTCCGGTGCCG‑TTCATTCTCGCCCGCACCTACAGCAG‑ATACCGTACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCTCCGGCGGGAAAATGCCTGCGGATATCCGC‑GTACAGCTGCGCGATAacacac                                             <  1:356617/251‑1 (MQ=11)                                                                                                                                                                                                cATCGGTGCCCCCACCGGCGTTGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTTCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAacactc                                             <  1:658248/251‑3 (MQ=255)                                                                                                                                                                                                cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACCAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAacacac                                             >  1:680190/1‑251 (MQ=37)                                                                                                                                                                                                cATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAacacac                                             >  1:533521/1‑251 (MQ=37)                                                                                                                                                                                                 aTCGGCGCCACCAACGGCGTGGCCTGTACGGTGTGC‑CCAAGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑CACGGTCCTTCCCGGTGAAACCGACATCCCCCTGCCCGCGCCGCTGCCG‑TTACTTCTCTCCCGCCCCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGGGCCTC‑GGCCCCGGCGGGAGAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAacaacc                                             >  2:480519/1‑247 (MQ=11)                                                                                                                                                                                                   cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGa                                          <  1:29025/251‑1 (MQ=35)                                                                                                                                                                                                   cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGa                                          <  1:587013/251‑1 (MQ=35)                                                                                                                                                                                                             aCCGGCGTGGCCTGTTCGGGGTGC‑CCCGGCGGGGTGACATCCGGCCATCCGGTCCAT‑CCGCTGCTCGGGGA‑CAAGGGCCTTCCCGGTGAAACCGACCTCGCCCTGCCCGGCCCGCTGCCC‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCCCGCCGGTGGGGCGCCCC‑GGCCCCCGCCGGGTAATGCCCGCGGGTCGCCGCCTTAC‑GCTGCGCGCTACCACCCTGCTACTCAgcgg                                >  2:270341/1‑248 (MQ=2)                                                                                                                                                                                                                   gTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcgg                          <  1:322673/251‑1 (MQ=34)                                                                                                                                                                                                                   gTGGCCTGCTCGGCCTGC‑CCCGGCGGGGTGACGTCCGGCCCCCCGATCACT‑CCGCTCCTCCGCGC‑AAAGGCCACTCCCGGTGAAACCGCCCTCGCCCCACCCGGCCCGCTGCAG‑TTCATTCTCTCACGCCCCTACAGCAG‑TTACCGGACAAAAAAGCCCGCGCCGGGGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGCATATCCGC‑CTCCAACTGCGCCGTAACACCCTGATaatccgggataaa                            >  2:658587/1‑237 (MQ=2)                                                                                                                                                                                                                               gtgtGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAA‑GCCt             <  1:367176/251‑1 (MQ=34)                                                                                                                                                                                                                                    c‑cccGGCGGGGTGACGTCCGGCCAGCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAA‑GCCTGTAtt        <  2:328759/251‑1 (MQ=21)                                                                                                                                                                                                                                       ccGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGAGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGCTATCCGC‑TTACAGATGCGCGATAACACACTGATACTCAGTGATAACGGCGGGAGAACGCC‑GTAtttt      >  2:532246/1‑251 (MQ=11)                                                                                                                                                                                                                                            ggggTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAA‑GCCTGTATTTTGAGc  >  2:241580/1‑250 (MQ=34)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             ATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGG‑CGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGCATCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGC‑CCCGGCGGGGTGACGTCCGGCCATCCGGTCAAT‑CCGCTGCTCGGTGC‑AAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCG‑TTCATTCTCTCCCGCACCTACAGCAG‑TTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTC‑GGCCCCGGCTGGAAAATGCCTGCGGATATCCGC‑TTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAA‑GCCTGTATTTTGAGC  >  NC_000913/729473‑729955

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