Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,182,585	C→G	16.4%	A13G (GCA→GGA)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,585	0	C	G	16.4%	177.1 / 2.6	75	A13G (GCA→GGA)	zupT	zinc transporter
Reads supporting (aligned to +/- strand):  ref base C (14/46);  new base G (12/2);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 2.50e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTCATCGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCA  >  NC_000913/3182351‑3182819                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                cGTTCATCGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑Gcc                                                                                                                                                                                                                                            <  2:396006/239‑1 (MQ=255)        cGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTa                                                                                                                                                                                                                                     <  1:155210/239‑1 (MQ=255)            cTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCG‑GGCACGtcgtatag                                                                                                                                                                                                                                 >  1:584447/1‑232 (MQ=255)              aCCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTTgc‑g                                                                                                                                                                                                                              >  2:1458590/1‑239 (MQ=255)                  tGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCGCGGTTAttttggcgtt                                                                                                                                                                                                                           >  2:840143/1‑231 (MQ=255)                  tGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑CGGCG‑GCCACATCTTTTTgc‑gccct                                                                                                                                                                                                                          >  1:1542755/1‑236 (MQ=255)                     ccTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGCTCTGCC‑TGAAGGTGAAGATATATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATCCTGGCGG‑GGGGG‑GCCACACGTTTTGGC‑GCG‑CTtc                                                                                                                                                                                                                       >  1:66016/1‑238 (MQ=255)                         aaaaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCAGGCCCCGTTTTTTgcgcctgtttttg                                                                                                                                                                                                                      >  1:362633/1‑226 (MQ=255)                         aaaaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTTGC‑GCG‑TTTCTCGGc                                                                                                                                                                                                                  >  2:1578241/1‑239 (MQ=255)                         aaaaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGAGAAGATATATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGGTGGGAG‑GCTACGTTT‑TTGgc‑ggctctcttcgc                                                                                                                                                                                                                   >  1:1553171/1‑228 (MQ=255)                              aaTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GACACGTTTATTTGC‑GCG‑TTTTTCGGCGTtc                                                                                                                                                                                                              >  2:559254/1‑238 (MQ=255)                               aTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTtc‑t‑c                                                                                                                                                                                                          <  1:1158869/239‑1 (MQ=255)                                tGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCCACTTTATTGGC‑GCG‑TTTCTCGGCGttt                                                                                                                                                                                                              >  2:1209153/1‑235 (MQ=255)                                     gCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGG‑GCCACGTTTATTTGG‑CCG‑TTTCTTCTCGTTC‑T‑CCGGcca                                                                                                                                                                                                    >  1:1236046/1‑237 (MQ=255)                                               tGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GGCACGTTTTTTAGC‑GCG‑TGTCTCTGCGTTA‑TATGGGCaaaagactcga                                                                                                                                                                                           >  1:1730454/1‑232 (MQ=255)                                                gAAGACATGATGATATCCTTAACTAAAGGAGTCAGTGTGATATCCTCACAATACGCTGGTTCGGAGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGaa                                                                                                                                                                                         <  1:187678/239‑1 (MQ=255)                                                       tgatgaTATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTa                                                                                                                                                                                  <  1:440991/239‑1 (MQ=255)                                                        gatgatATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGGTTATTGGC‑GCG‑TTTTTCCGCGTTC‑T‑CGGGCAAAAACCCCCCAACCGCtt                                                                                                                                                                                   >  2:857048/1‑237 (MQ=255)                                                              atCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GCCCCGTTTTTTGGC‑GCG‑TTTTTCTGCGGTC‑T‑CCGGCAAAAACCCTCGAACCGCTTACTGGCgg                                                                                                                                                                           >  2:206906/1‑238 (MQ=255)                                                              atCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑GTTTTTCTCGGTT‑T‑CGGGCAAAAAACCCCGAACCGCTTACTGGcct                                                                                                                                                                           >  2:1541162/1‑237 (MQ=255)                                                                   tAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GCCCCGTTTTTTTGT‑GGG‑TTTCTCTGCGTTC‑T‑CGGGGAAAAAACCCCCGACCCCCTACTGGCGttttt                                                                                                                                                                       >  1:574726/1‑237 (MQ=255)                                                                    aaCTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGCGTAAGAACCCGGCTAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCCCGTTTATTGCC‑GCG‑TTTCTCGCCGTGG‑T‑CGGGCTCAAACCCTCGTACCGCCTACTGCCGtcttt                                                                                                                                                                       >  2:1316174/1‑236 (MQ=255)                                                                     aCTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCt                                                                                                                                                                    <  2:1730454/239‑1 (MQ=255)                                                                         aaGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCCCGTTTATTTGC‑GCG‑TTTCTCCGCGTTT‑T‑CGGGCAAAAACCCCCGAACCGCCTTCTGGCGTTTTCGCTTGGt                                                                                                                                                                >  2:1000315/1‑239 (MQ=255)                                                                         aaGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑AGGGA‑GCCACGTTTATTTGC‑GCG‑TTTTTCTGCGTTC‑T‑CGGGGAGAAAACCCCGAACCGCTTTCTGGCGGTTTctccctgt                                                                                                                                                                >  1:1091390/1‑232 (MQ=255)                                                                              gtcagtgttatatCCTCACAAGTCGCTTGTTCGGCGGAGTAAGAACACGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAATGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGACA‑GCCAAGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CTGGCAAAATCCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTg                                                                                                                                                            <  1:340326/229‑1 (MQ=255)                                                                              gtCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GCCAAGGTTTTTGGC‑GCG‑TTTCTTGGCGTTCTT‑CGGGGAAAAA‑CCTCCGACCGCTTACTGGCGGTTTTGCTAGGTTTTgc                                                                                                                                                           >  1:528110/1‑239 (MQ=255)                                                                              gtCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCCCGTTTTTTTGC‑CCG‑GTTTTCGGCGTTC‑T‑CCGGCAAAAACCCCCCAACCCCTTTCCTGGGTTTTctgttcgtttttg                                                                                                                                                           >  1:1558893/1‑227 (MQ=255)                                                                               tCATTTTTATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGC‑ACCCCGGTTTTTTGT‑GCG‑TTTCTCAGCCGTT‑T‑CGGGCCAGCACCCCCCACCCCCTTTCTTGTCTTTTTGTTAGGttttttg                                                                                                                                                          >  1:1207725/1‑236 (MQ=255)                                                                                cATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg                                                                                                                                                         <  2:1651069/239‑1 (MQ=255)                                                                                cATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcgg                                                                                                                                                         <  1:761809/239‑1 (MQ=255)                                                                                 aTTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCG‑GCCACGTTTATTTGC‑GCG‑GTTCTCGGCGTTC‑T‑CGGGCAAAAAACCTCCAAACCCCTTCTGGCGTTTTTGCTTGGtttttccgc                                                                                                                                                        >  2:782431/1‑234 (MQ=255)                                                                                 aTTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggc                                                                                                                                                        <  2:513934/239‑1 (MQ=255)                                                                                  ttttGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTgcggcg                                                                                                                                                       <  2:1091390/239‑1 (MQ=255)                                                                                       atatCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GCCCCGTTTTTTTGC‑GCGTTTTTTCGCCTTTC‑C‑GCGGCAAAAA‑CCTCCGACCGCTTTCCTGCGTTTTTGCTCGGTTTTGCGGCgggggg                                                                                                                                                  >  1:638142/1‑237 (MQ=255)                                                                                         atCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCa                                                                                                                                                <  2:143318/239‑1 (MQ=255)                                                                                         atCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGGTTTTTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCCCGAACCGCTTACTGGCGTTTTCGCTAGGGTTTTCGGCggggggc                                                                                                                                                 >  2:912737/1‑235 (MQ=255)                                                                                              cacaATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTg                                                                                                                                           <  2:940042/239‑1 (MQ=255)                                                                                                caATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGct                                                                                                                                         <  1:1578241/239‑1 (MQ=255)                                                                                                   tACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTCCCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcat                                                                                                                                      <  1:619737/239‑1 (MQ=255)                                                                                                   tACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGctcat                                                                                                                                      <  2:574726/239‑1 (MQ=255)                                                                                                     cGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGGA‑GCCACGGTTATTTGT‑GGG‑TTTCTCGGCGGTT‑T‑CGGGGAAAAACACCCGAACCGCTTTCTGGCGCTTTTGTTTGGGTTTGcggcggcggggagcttgttcttt                                                                                                                                     >  1:263449/1‑221 (MQ=255)                                                                                                        ttGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑AGGGA‑GCCCCGTTTATTTGC‑GCG‑GTTCTCGGCGTTC‑T‑CTCGCAAAAAACCCCGCACCGCTTCTTGGCGGTTTTGCTAGGTTTTGCGGCgggggggctgtttgtcctctc                                                                                                                                   >  1:794202/1‑220 (MQ=255)                                                                                                         tGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCAtt                                                                                                                                <  2:559948/239‑1 (MQ=255)                                                                                                            tCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAAt                                                                                                                             <  2:704272/239‑1 (MQ=255)                                                                                                                 ggAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGaaa                                                                                                                        <  2:1553171/239‑1 (MQ=255)                                                                                                                                 tAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACt                                                                                                        <  1:1007526/239‑1 (MQ=255)                                                                                                                                  aaCAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTa                                                                                                       <  1:990159/239‑1 (MQ=255)                                                                                                                                   aCAAGGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTTCCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTag                                                                                                      <  1:286775/239‑1 (MQ=255)                                                                                                                                       tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                  <  2:1558893/239‑1 (MQ=255)                                                                                                                                       tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                  <  2:335916/239‑1 (MQ=255)                                                                                                                                       tgatgatgaTCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTagcagc                                                                                                  <  2:528110/239‑1 (MQ=255)                                                                                                                                          tgatgatCAGCAGTTATTTTGGCGATCTGCC‑TGAAGGTGAAGATTTATAAGGGGTTGTCGATGTCAGTACCCCTCATTCTGACGATCCTGGCGG‑GGGGA‑GCCACGGTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCCAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTg                                                                                                <  1:1545508/238‑1 (MQ=255)                                                                                                                                                tcatcaGTTATTTTTGCGGTCTGCCGTGAA‑GTGAAGGTTTAAAAGGGGGTGGCGATTGCAGTTCCTCCCATTCCTACCCTTCTGGGGG‑GGGGC‑GCCCAGTTTATTTGCGGCG‑TTT‑TCCGGGTTC‑T‑CGGGCAAAAACCCCCGAACCCGTTTCCGGCGGTTTCGCTAGGGTTTGGGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAAt                                                                                         <  1:1476488/239‑1 (MQ=255)                                                                                                                                                   tcaGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTc                                                                                      <  1:1458590/239‑1 (MQ=255)                                                                                                                                                        tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                    <  1:357583/236‑1 (MQ=255)                                                                                                                                                        tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                    <  2:868025/236‑1 (MQ=255)                                                                                                                                                        tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                    <  2:454714/236‑1 (MQ=255)                                                                                                                                                        tATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                    <  2:1058814/236‑1 (MQ=255)                                                                                                                                                          ttttGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCtgtg                                                                               <  2:66016/239‑1 (MQ=255)                                                                                                                                                            ttGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGtt                                                                             <  2:690037/239‑1 (MQ=255)                                                                                                                                                                    cggCC‑TGAAGGTGAAGATTTATAAGGAGTGGGCGATGTCAGTCCCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTTGC‑GCG‑TTTCTCGGCGTTT‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                                     <  1:786254/237‑1 (MQ=255)                                                                                                                                                                    cTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATg                                                                     <  1:1320802/239‑1 (MQ=255)                                                                                                                                                                       cc‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGc                                                                                    <  2:218225/221‑1 (MQ=255)                                                                                                                                                                                     aTTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTag                                                                                                      <  2:294606/190‑1 (MQ=255)                                                                                                                                                                                                   gTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtt                                       <  2:263449/239‑1 (MQ=255)                                                                                                                                                                                                   gTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAtt                                       <  2:942751/239‑1 (MQ=255)                                                                                                                                                                                                    tCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTgg                                            <  2:642212/233‑1 (MQ=255)                                                                                                                                                                                                        tGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGc                                  <  1:1209153/239‑1 (MQ=255)                                                                                                                                                                                                           cAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCt                                <  2:1640103/238‑1 (MQ=255)                                                                                                                                                                                                               aCCTGTCATTCTGACCATTCTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGAcc                           <  2:1259486/239‑1 (MQ=255)                                                                                                                                                                                                                 ctctCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGc                         <  2:328327/239‑1 (MQ=255)                                                                                                                                                                                                                                 aCTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTc         <  1:1533264/239‑1 (MQ=255)                                                                                                                                                                                                                                       gg‑gggCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg    <  2:1534643/238‑1 (MQ=255)                                                                                                                                                                                                                                          gggCA‑GCCACGTTTAGTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCa  <  1:559254/238‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                CGTTCATCGGACTACCGTGACCTAAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCC‑TGAAGGTGAAGATTTATAAGGAGTTGTCGATGTCAGTACCTCTCATTCTGACCATACTGGCGG‑GGGCA‑GCCACGTTTATTGGC‑GCG‑TTTCTCGGCGTTC‑T‑CGGGCAAAAACCCTCGAACCGCTTACTGGCGTTTTCGCTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCA  >  NC_000913/3182351‑3182819

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