Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,182,601	G→T	17.3%	G18C (GGC→TGC)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,601	0	G	T	17.3%	196.7 / 5.9	72	G18C (GGC→TGC)	zupT	zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (9/49);  new base T (11/2);  total (21/51)
Fisher's exact test for biased strand distribution p-value = 3.80e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAT  >  NC_000913/3182366‑3182838                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  aCCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTTgc‑g                                                                                                                                                                                                                                               >  1:1458590/1‑239 (MQ=255)     tGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑CGGCGGCCACATCTTTTTgc‑gccct                                                                                                                                                                                                                                           >  2:1542755/1‑236 (MQ=255)        ccTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGCTCTGCC‑TGAAGGTGAAGATATATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATCCTGGCGG‑GGGGGGCCACACGTTTTGGC‑GCG‑CTtc                                                                                                                                                                                                                                        >  2:66016/1‑238 (MQ=255)            aaaaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTTGC‑GCG‑TTTCTCGGc                                                                                                                                                                                                                                   >  1:1578241/1‑239 (MQ=255)            aaaaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGAGAAGATATATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGTGGGAGGCTACGTTT‑TTGgc‑ggctctcttcgc                                                                                                                                                                                                                                    >  2:1553171/1‑228 (MQ=255)                 aaTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGACACGTTTATTTGC‑GCG‑TTTTTCGGCGTtc                                                                                                                                                                                                                               >  1:559254/1‑238 (MQ=255)                  aTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTtc‑t‑c                                                                                                                                                                                                                           <  2:1158869/239‑1 (MQ=255)                   tGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCCACTTTATTGGC‑GCG‑TTTCTCGGCGttt                                                                                                                                                                                                                               >  1:1209153/1‑235 (MQ=255)                        gCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGGGCCACGTTTATTTGG‑CCG‑TTTCTTCTCGTTC‑T‑CCGGcca                                                                                                                                                                                                                     >  2:1236046/1‑237 (MQ=255)                                  tGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGGCACGTTTTTTAGC‑GCG‑TGTCTCTGCGTTA‑TATGGGCaaaagactcga                                                                                                                                                                                                            >  2:1730454/1‑232 (MQ=255)                                   gAAGACATGATGATATCCTTAACTAAAGGAGTCAGTGTGATATCCTCACAATACGCTGGTTCGGAGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGaa                                                                                                                                                                                                          <  2:187678/239‑1 (MQ=255)                                          tgatgaTATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTa                                                                                                                                                                                                   <  2:440991/239‑1 (MQ=255)                                           gatgatATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGGTTATTGGC‑GCG‑TTTTTCCGCGTTC‑T‑CGGGCAAAAACCCCCCAACCGCtt                                                                                                                                                                                                    >  1:857048/1‑237 (MQ=255)                                                 atCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGCCCCGTTTTTTGGC‑GCG‑TTTTTCTGCGGTC‑T‑CCGGCAAAAACCCTCGAACCGCTTACTGGCgg                                                                                                                                                                                            >  1:206906/1‑238 (MQ=255)                                                 atCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑GTTTTTCTCGGTT‑T‑CGGGCAAAAAACCCCGAACCGCTTACTGGcct                                                                                                                                                                                            >  1:1541162/1‑237 (MQ=255)                                                      tAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGCCCCGTTTTTTTGT‑GGG‑TTTCTCTGCGTTC‑T‑CGGGGAAAAAACCCCCGACCCCCTACTGGCGttttt                                                                                                                                                                                        >  2:574726/1‑237 (MQ=255)                                                       aaCTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGCGTAAGAACCCGGCTAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCCCGTTTATTGCC‑GCG‑TTTCTCGCCGTGG‑T‑CGGGCTCAAACCCTCGTACCGCCTACTGCCGtcttt                                                                                                                                                                                        >  1:1316174/1‑236 (MQ=255)                                                        aCTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCt                                                                                                                                                                                     <  1:1730454/239‑1 (MQ=255)                                                            aaGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCCCGTTTATTTGC‑GCG‑TTTCTCCGCGTTT‑T‑CGGGCAAAAACCCCCGAACCGCCTTCTGGCGTTTTCGCTTGGt                                                                                                                                                                                 >  1:1000315/1‑239 (MQ=255)                                                            aaGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑AGGGAGCCACGTTTATTTGC‑GCG‑TTTTTCTGCGTTC‑T‑CGGGGAGAAAACCCCGAACCGCTTTCTGGCGGTTTctccctgt                                                                                                                                                                                 >  2:1091390/1‑232 (MQ=255)                                                                 gtcagtgttatatCCTCACAAGTCGCTTGTTCGGCGGAGTAAGAACACGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAATGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGACAGCCAAGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CTGGCAAAATCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTg                                                                                                                                                                             <  2:340326/229‑1 (MQ=255)                                                                 gtCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGCCAAGGTTTTTGGC‑GCG‑TTTCTTGGCGTTCTT‑CGGGGAAAAA‑CCTCCGACCGCTTACTGGCGGTTTTGCTAGGTTTTgc                                                                                                                                                                            >  2:528110/1‑239 (MQ=255)                                                                  tCATTTTTATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGCACCCCGGTTTTTTGT‑GCG‑TTTCTCAGCCGTT‑T‑CGGGCCAGCACCCCCCACCCCCTTTCTTGTCTTTTTGTTAGGttttttg                                                                                                                                                                           >  2:1207725/1‑236 (MQ=255)                                                                   cATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg                                                                                                                                                                          <  1:1651069/239‑1 (MQ=255)                                                                   cATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg                                                                                                                                                                          <  2:761809/239‑1 (MQ=255)                                                                    aTTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCGGCCACGTTTATTTGC‑GCG‑GTTCTCGGCGTTC‑T‑CGGGCAAAAAACCTCCAAACCCCTTCTGGCGTTTTTGCTTGGtttttccgc                                                                                                                                                                         >  1:782431/1‑234 (MQ=255)                                                                    aTTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggc                                                                                                                                                                         <  1:513934/239‑1 (MQ=255)                                                                     ttttGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                                                                                                                                                                        <  1:1091390/239‑1 (MQ=255)                                                                          atatCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGAGCCCCGTTTTTTTGC‑GCGTTTTTTCGCCTTTC‑C‑GCGGCAAAAA‑CCTCCGACCGCTTTCCTGCGTTTTTGCTCGGTTTTGCGGCgggggg                                                                                                                                                                   >  2:638142/1‑237 (MQ=255)                                                                            atCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCa                                                                                                                                                                 <  1:143318/239‑1 (MQ=255)                                                                            atCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGGTTTTTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCCCGAACCGCTTACTGGCGTTTTCGCTAGGGTTTTCGGCggggggc                                                                                                                                                                  >  1:912737/1‑235 (MQ=255)                                                                                 cacaATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTg                                                                                                                                                            <  1:940042/239‑1 (MQ=255)                                                                                   caATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGct                                                                                                                                                          <  2:1578241/239‑1 (MQ=255)                                                                                      tACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTCCCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcat                                                                                                                                                       <  2:619737/239‑1 (MQ=255)                                                                                      tACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcat                                                                                                                                                       <  1:574726/239‑1 (MQ=255)                                                                                            tGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCAtt                                                                                                                                                 <  1:559948/239‑1 (MQ=255)                                                                                               tCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAAt                                                                                                                                              <  1:704272/239‑1 (MQ=255)                                                                                                    ggAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGaaa                                                                                                                                         <  1:1553171/239‑1 (MQ=255)                                                                                                                    tAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACt                                                                                                                         <  2:1007526/239‑1 (MQ=255)                                                                                                                     aaCAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTa                                                                                                                        <  2:990159/239‑1 (MQ=255)                                                                                                                      aCAAGGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTTCCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTag                                                                                                                       <  2:286775/239‑1 (MQ=255)                                                                                                                          tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                   <  1:528110/239‑1 (MQ=255)                                                                                                                          tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                   <  1:335916/239‑1 (MQ=255)                                                                                                                          tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                                   <  1:1558893/239‑1 (MQ=255)                                                                                                                             tgatgatCAGCAGTTATTTTGGCGATCTGCC‑TGAAGGTGAAGATTTATAAGGGGTTGTCGATGTCAGTACCCCTCATTCTGACGATCCTGGCGG‑GGGGAGCCACGGTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCCAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTg                                                                                                                 <  2:1545508/238‑1 (MQ=255)                                                                                                                                   tcatcaGTTATTTTTGCGGTCTGCCGTGAA‑GTGAAGGTTTAAAAGGGGGTGGCGATTGCAGTTCCTCCCATTCCTACCCTTCTGGGGG‑GGGGCGCCCAGTTTATTTGCGGCG‑TTT‑TCCGGGTTC‑T‑CGGGCAAAAACCCCCGAACCCGTTTCCGGCGGTTTCGCTAGGGTTTGGGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAAt                                                                                                          <  2:1476488/239‑1 (MQ=255)                                                                                                                                      tcaGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTc                                                                                                       <  2:1458590/239‑1 (MQ=255)                                                                                                                                           tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                                     <  1:868025/236‑1 (MQ=255)                                                                                                                                           tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                                     <  2:357583/236‑1 (MQ=255)                                                                                                                                           tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                                     <  1:1058814/236‑1 (MQ=255)                                                                                                                                           tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                                     <  1:454714/236‑1 (MQ=255)                                                                                                                                             ttttGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCtgtg                                                                                                <  1:66016/239‑1 (MQ=255)                                                                                                                                               ttGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGtt                                                                                              <  1:690037/239‑1 (MQ=255)                                                                                                                                                       cggCC‑TGAAGGTGAAGATTTATAAGGAGTGGGCGATGTCAGTCCCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTTGC‑GCG‑TTTCTCGGCGTTT‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                                                      <  2:786254/237‑1 (MQ=255)                                                                                                                                                       cTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                                                      <  2:1320802/239‑1 (MQ=255)                                                                                                                                                          cc‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                                     <  1:218225/221‑1 (MQ=255)                                                                                                                                                                        aTTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTag                                                                                                                       <  1:294606/190‑1 (MQ=255)                                                                                                                                                                                      gTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtt                                                        <  1:942751/239‑1 (MQ=255)                                                                                                                                                                                      gTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtt                                                        <  1:263449/239‑1 (MQ=255)                                                                                                                                                                                       tCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                                                             <  1:642212/233‑1 (MQ=255)                                                                                                                                                                                           tGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGc                                                   <  2:1209153/239‑1 (MQ=255)                                                                                                                                                                                              cAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                                                 <  1:1640103/238‑1 (MQ=255)                                                                                                                                                                                                  aCCTGTCATTCTGACCATTCTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc                                            <  1:1259486/239‑1 (MQ=255)                                                                                                                                                                                                    ctctCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc                                          <  1:328327/239‑1 (MQ=255)                                                                                                                                                                                                                    aCTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTc                          <  2:1533264/239‑1 (MQ=255)                                                                                                                                                                                                                          gg‑gggCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg                     <  1:1534643/238‑1 (MQ=255)                                                                                                                                                                                                                             gggCAGCCACGTTTAGTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCa                   <  2:559254/238‑1 (MQ=255)                                                                                                                                                                                                                                    cACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAttt             <  1:1246784/237‑1 (MQ=255)                                                                                                                                                                                                                                       gTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc        <  2:912737/239‑1 (MQ=255)                                                                                                                                                                                                                                        tttATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCa       <  2:49222/239‑1 (MQ=255)                                                                                                                                                                                                                                         ttATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCa       <  2:782431/238‑1 (MQ=255)                                                                                                                                                                                                                                             tGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAt  <  1:1152510/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  ACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAT  >  NC_000913/3182366‑3182838

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