Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	360,362	T→G	26.9%	V124G (GTG→GGG)	cynX →	putative cyanate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	360,362	0	T	G	26.9%	70.7 / 2.6	47	V124G (GTG→GGG)	cynX	putative cyanate transporter
Reads supporting (aligned to +/- strand):  ref base T (9/23);  new base G (11/4);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 5.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTACCATGGGCGGGCTGGCGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTCCCCGCGCGTGGACGCTG  >  NC_000913/360127‑360594                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             gtTACCATGGGCGGGCTGGCGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGGGATTAAACGGC‑GGGTTCACCCGCGCCCG‑CCACCGGGGAt                                                                                                                                                                                                                                           >  2:979535/1‑239 (MQ=255)   tacctggggcgggcTGGGGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATg                                                                                                                                                                                                                                          <  2:258897/232‑1 (MQ=255)     ccATGGGCGGGCTGGCGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAAAGGC‑GGTTTTCGCAGGGCACC‑CCACCGGGGATtggg                                                                                                                                                                                                                                       >  2:935523/1‑239 (MQ=255)         gggcgggcTGGCGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGGGATTAAACGGC‑GGTTTCCGCAGCGCACG‑CCACCGGTGATGGGGg‑tgt                                                                                                                                                                                                                                  >  2:888793/1‑239 (MQ=255)                  ggggctggcCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGGGCTTTACCCGCAAAGGGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCgg                                                                                                                                                                                                                        <  2:565423/236‑1 (MQ=255)                            ggAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTTTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATggg                                                                                                                                                                                                              <  2:480962/239‑1 (MQ=255)                                     cTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCCGGAGGGCACG‑CCCCTGGGGATGGGGC‑CGTGGT‑CCGCGGCTCTGGTGGGGGGGGGTgg                                                                                                                                                                                                     >  2:60870/1‑239 (MQ=255)                                        catcaGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGTTGGGGTGGGGATCGGCATCATTCAGGCGGGGCTGCCCTCGGGGATTAAACGGCGGGTTTCA‑CAGCCGACG‑GCCCCGGGGATTGGGC‑TGTGGTGCCGGGGGT‑TGGTGGGCGgggggggggg                                                                                                                                                                                                  >  1:1243552/1‑237 (MQ=255)                                               aTGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGGCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTgc                                                                                                                                                                                           <  2:719058/238‑1 (MQ=255)                                                   cAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCCCG‑CCCCTGGGGATGGGGCTTGTGGC‑CCG‑GGGACTGGTGGGGGGCGGGGGG‑GTTCGGggcgccc                                                                                                                                                                                      >  2:429261/1‑238 (MQ=255)                                                   cAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGGGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGGGATGGGGC‑TGTGGT‑CCGCGGCCCTGATGGGGGGCGGAGGG‑CTTGGTggcgccc                                                                                                                                                                                      >  2:105228/1‑238 (MQ=255)                                                        aaCGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGGTTCAGAAGCGCACG‑CCCCTGGTGATGGGGC‑TGTGGG‑CCCCGGCCCTGTTGGGGGGGGGGGGG‑CTTGTGGCCGCCCTAAc                                                                                                                                                                                  >  2:281436/1‑238 (MQ=255)                                                            tCGCAGTGTCGCCATCAGACTGTTGCTGATTGCCGTCGGTGCATTGAGGCGTGAGCTTTACCCGCAAAGTGCGCTGCGGCTTAGCGGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑Gcc                                                                                                                                                                              <  2:466378/237‑1 (MQ=255)                                                                agtggcggcggccgtttttgggtGATTGCGGTCGGGGGATTGGTGCGTGGGCTTTCGCCGAAAAGCGCGCTGCTGCTTGGGCGCGCACGGCTTGGTGGGGGGGGGAGCGCGAGCCGTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGTTGAGGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGt                                                                                                                                                                         <  1:281436/216‑1 (MQ=255)                                                                    tCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCCCG‑CCACTGGGGATGGGGC‑TGGGGT‑CCCGGGGTCTGGTGGGCGGCGGGGGG‑GTTGTGGCCGCCATaagcacccgggtagt                                                                                                                                                                      >  2:793304/1‑225 (MQ=255)                                                                     cGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGGTGGGCGGCGGGGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAggtt                                                                                                                                                                   >  2:83781/1‑236 (MQ=255)                                                                         ttcctttgtttGCTTATTGCCTTGGGTGCATTGGTGCGGTGGCTTTCCCCGCAAGGGGGGCTGGTGCTTAGCAGGGCACTGTTTGGTGGGGTGGGGATCGGCAGCATTCAGGCGGGGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACa                                                                                                                                                               <  1:390802/230‑1 (MQ=255)                                                                         ctcggttgtttgcctatgcgcGTCGGTGCATTGGTGCGTGAGCTTTACCCGCAAGGGGCGCTGCGGCTTAGAAGCGCACTGCTTGGTGGGGTGGGGATCGGTATCATTCAGGCGGTGATGCCTTGGGGGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAAc                                                                                                                                                                <  2:46843/218‑1 (MQ=255)                                                                                  ttGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGCTTAAACGGC‑GGGTTCAGACCGGCGCGTCCAC‑GGTTGTGGGGG‑CGTGGG‑CCGCGGCTCTGGTGGGGGGCGGTGGG‑GATTGGGCCGCCCTAACGGACCGGGTTAG‑TCAAAAAAGCGaga                                                                                                                                                       >  1:1048366/1‑236 (MQ=255)                                                                                        aTTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAACGCACG‑CCCCTGGGGGTGGGGC‑CGTGGT‑CCCCGGCCCTTATGGGCGGCGGTGGG‑CTTGGGGCCGCCATAAC‑GCCCTGGTTTGTTCAACaactgtctcttttta                                                                                                                                                 >  2:243564/1‑224 (MQ=255)                                                                                                                       aCCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGGTTCACGAGCGCTCC‑CCCCTGGGGCTTGGGC‑CGGGTG‑CCGCTGGCATCGTGGGGGGCGGGGGGCCTTGTTGCCGCGCTAAC‑GCCCGGGTTAGTTCAACAAAGCGAAACCGGGTCTCAAACACCCCCCGGGGGGGGGctgccgcg                                                                                                                  >  1:22968/1‑236 (MQ=255)                                                                                                                              aaGTGCGCTGCTGCTTAGCAGGGCACGGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGttgttg                                                                                                          <  2:1161955/239‑1 (MQ=255)                                                                                                                                                             ggtggggtggGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACAAAGCGAAACCTGGGATCCAACACTCCCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCCCTTTTCCTGGTGGGGGCAAAGCGccc                                                                            >  2:435059/1‑238 (MQ=255)                                                                                                                                                               tgggggggggTCCGGCATCATTCAGGCGGTGATGCTTTCGGTGATTAAACGGC‑GTTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGGGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCcgc                                                                          <  2:703853/238‑1 (MQ=255)                                                                                                                                                                    tggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATa                                                                                                                                                                                    <  2:906993/128‑1 (MQ=255)                                                                                                                                                                    tggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCCCTGATGGGCGGCGGGGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGGTAGTTCAACATAGCGGAACCTGGGATCAAACACTCCCCTGGTGGGCGCTGCCTGCCCGTGTTTCGCCCTTTGCCTGTGGGGGGCAAAACGCCCCCggggg                                                                     >  2:564343/1‑237 (MQ=255)                                                                                                                                                                    gggggggcGGCTTCTTTCGGGCGGGGGGGCCTTGGGTGGTTAACGGGC‑GTTTCAAGGACGGCACG‑CCACTGGTGGGGGGTT‑TGGTGG‑CGGCGGCTCGGATGGGGGGCGGTGGG‑TTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGt                                                                     <  1:979535/231‑1 (MQ=255)                                                                                                                                                                     ggggtccGGAATCATTCAGGGGGTGATGCCTTCGGGGATTAACCGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCg                                                                   <  2:292456/233‑1 (MQ=255)                                                                                                                                                                     ggggACCGGCATCGTTCAGGCGGTGATGCCTTGGGTTATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCg                                                                   <  2:544795/239‑1 (MQ=255)                                                                                                                                                                          tCGGCATGATTCAGGCGGTGATGCCTTGGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCctct                                                              <  2:22968/239‑1 (MQ=255)                                                                                                                                                                              catcatTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGGTTCAGCAGCGCACG‑CCACTGGGGATTGGGC‑TGTGGG‑TCGCGGCTCTGATGGGCGGCGGGGGG‑CTTGGGGCCGCCACAAC‑GCCCTGGt                                                                                                                                                                         >  1:906993/1‑128 (MQ=255)                                                                                                                                                                                                    tgggTGATTAAAGGGC‑GGTTTCAGCAGCGCACG‑CCCTGGGTGATGGGGC‑TGTGGT‑CCGCGGCTCGGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGt                                    <  2:1065498/237‑1 (MQ=255)                                                                                                                                                                                                      ggTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGGGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTTGTGGTGGCAAAGCGCCCGCGAGGGCGCCTCTTCCCCCAAGACAACACCCACTCCCGTTc                                  >  2:674197/1‑239 (MQ=255)                                                                                                                                                                                                           ttAACCGGC‑GGTTTCAGCAGGCCACG‑CCAGTGGTGGTGGGGC‑TGTGGT‑CCGGGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGt                              <  2:1243552/238‑1 (MQ=255)                                                                                                                                                                                                           ttAAAGGGT‑GGTTGCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTTGG‑CTGGGTGCCGCCATAAC‑GCCCTTGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGt                              <  2:686563/238‑1 (MQ=255)                                                                                                                                                                                                           ttAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TTGTGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGt                              <  2:256433/238‑1 (MQ=255)                                                                                                                                                                                                           ttAAACGGC‑GGTTTAAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTTGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACCCTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGt                              <  2:1029177/238‑1 (MQ=255)                                                                                                                                                                                                            ttaacggggggTTTCGGCAGCGCACG‑CCACTGGGGATGGGGC‑TGTGGT‑CCGCGGTTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGt                              <  2:1048366/229‑1 (MQ=255)                                                                                                                                                                                                             taaccggggtgttCAGCAGGGCACG‑CCACTGGGGATGGGGC‑TGGTGT‑CCGCTGCTCTGATTGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTgg                            <  2:1109030/228‑1 (MQ=255)                                                                                                                                                                                                                     ggTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑Gcc                                                                                                                                                                              >  1:518182/1‑85 (MQ=255)                                                                                                                                                                                                                     ggTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑Gcc                                                                                                                                                                              <  2:518182/85‑1 (MQ=255)                                                                                                                                                                                                                        tccggcagcGCCCG‑CCCCTGGGGATGGGGC‑TGGTGT‑CCGCCGCCTGGATGGGCGGCGGTGGG‑CTGGTTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACCCTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTccc                 <  2:143179/235‑1 (MQ=255)                                                                                                                                                                                                                               ggcgaCG‑CCGCTGGTGATGGGGT‑GGTGGT‑CCGCGGTTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTCCCcgcgcg           <  2:56805/234‑1 (MQ=255)                                                                                                                                                                                                                                cgcACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCCCGTGGAATTCACTCCCCGCCCGTg         >  2:138547/1‑239 (MQ=255)                                                                                                                                                                                                                                cgcACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGAATCACACACTCGCCTGGTGTGGGCGGCCTTCCGTTATTGCGCTCTTTGCCTGGGGGGGGAAAAGCGCCCGCGACGGCGCCTCTTCCCACAAGACACCAACCACTCCGGTTCCCGTGGTaagacatcccccccgggg         >  1:906470/1‑222 (MQ=255)                                                                                                                                                                                                                                        cccgggTGAGGGGGC‑TGGGGG‑CAGGGGTTATGAGGGGCGGCGGTGGG‑CTGGGTGCCGCCATAAC‑GCCCTGGTTAGTGAAACAGAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTCCCCGCGCGTGGACGCt   <  2:693116/234‑1 (MQ=255)                                                                                                                                                                                                                                        cccTGGGGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTCCCCGCGCGTGGACGCTg  <  2:846314/237‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             GTTACCATGGGCGGGCTGGCGCTGGCCGGAAGCTGGCTTCATCAGCATGTCAGCGAACGTCGCAGTGTCGCCATCAGTCTGTTGCTGATTGCCGTCGGTGCATTGATGCGTGAGCTTTACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGC‑GGTTTCAGCAGCGCACG‑CCACTGGTGATGGGGC‑TGTGGT‑CCGCGGCTCTGATGGGCGGCGGTGGG‑CTTGGTGCCGCCATAAC‑GCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCACAAGACAACAACCACTCCGGTTCGCGTGGTATTCACTCCCCGCGCGTGGACGCTG  >  NC_000913/360127‑360594

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