Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,031,160	C→T	E386E (GAG→GAA)	dcm ←	DNA‑cytosine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,031,160	0	C	T	100.0%	78.7 / NA	24	E386E (GAG→GAA)	dcm	DNA‑cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/17);  total (7/17)

GCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTC  >  NC_000913/2030914‑2031398                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               gCCATGTTGGGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCAc                                                                                                                                                                                                                                            <  1:254367/251‑1 (MQ=255)             ctctTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCatat                                                                                                                                                                                                                                >  2:8706/1‑251 (MQ=255)               ctTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGGGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCatat                                                                                                                                                                                                                                <  1:50313/249‑1 (MQ=255)                        gCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCcgcg                                                                                                                                                                                                                      >  1:213634/1‑250 (MQ=255)                               cACCGCCTGTTTGATTTTTGGATCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGAtt                                                                                                                                                                                                               >  1:222160/1‑250 (MQ=255)                                                               ttGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACgcgc                                                                                                                                                                               <  1:228937/250‑1 (MQ=255)                                                                                aCCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGAc                                                                                                                                                             >  1:35948/1‑251 (MQ=255)                                                                                   ggacctaccaccgattgcccGAACTTCCGATAGGCTTGAGTGTCCGACACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCCTTAATCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCTGATCGTCAAATCCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGActct                                                                                                                                                          <  2:183157/233‑4 (MQ=255)                                                                                     cACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTGCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCttt                                                                                                                                                        <  1:139967/251‑1 (MQ=255)                                                                                        gACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCg                                                                                                                                                     >  1:195866/1‑251 (MQ=255)                                                                                                         aCTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAAcc                                                                                                                                    <  1:52902/251‑1 (MQ=255)                                                                                                            ggCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTTCTGATGCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACTCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCa                                                                                                                                   <  1:16031/249‑1 (MQ=255)                                                                                                                          gtgtCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCGTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGcc                                                                                                                   <  2:213634/251‑1 (MQ=255)                                                                                                                                                         cccccGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCata                                                                                     <  2:195866/248‑1 (MQ=255)                                                                                                                                                                    tCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCAGATCGATAGAgg                                                                           <  2:35948/249‑1 (MQ=255)                                                                                                                                                                            cATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCcaca                                                                 <  1:8706/251‑1 (MQ=255)                                                                                                                                                                             aTTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCGTTTTCACCCGTGGCCATATCCCAGCCTCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAg                                                                <  1:169685/251‑1 (MQ=255)                                                                                                                                                                                    gcgcgcATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTCCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGc                                                         <  2:212321/251‑1 (MQ=255)                                                                                                                                                                                     cgcgcATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCg                                                        >  1:113895/1‑251 (MQ=255)                                                                                                                                                                                                        caCCGACGTGAACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCAACCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACTCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCGGATGTTTTTTAGCATATCGATAGAGGCACTTCCTCAGCCCCGGCGTCAGGATATATTTCGCCTc                                     >  2:9552/2‑251 (MQ=255)                                                                                                                                                                                                                   aCGAAGTTGCTGATTCAGCGGAACGTCCAAGACTTTTTCACCCGTGGCCATATCCAAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGCATAGATTTCGCCTCGACCATCGGGt                          <  2:175415/251‑1 (MQ=255)                                                                                                                                                                                                                        gTTGCTGATTCAGCGGATCGTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGGTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAAc                     <  1:118011/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGTCAAAGGCTTTTTCACCCGTGGCCATATCCCAGCCGCGAGCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGATGTACTTCCACAGCACCGGCGTCATGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACt   <  1:292181/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTCAAAGTCTTTTTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTc  <  1:66371/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               GCCATGTTGTGCCTCTTGCTGACGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGATCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGATAAACCATTCCATAACCGAAGCCGTTACCGCGCGCCTGATGTTTTTTCGCATATCGATAGAGGTACTTCCACAGCACCGGCGTCAGGATATATTTCGCCTCGACCATCGGGTCCAACAGCTGCGCCAGCGTCACTC  >  NC_000913/2030914‑2031398

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