Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,031,071	G→A	P416L (CCG→CTG)	dcm ←	DNA cytosine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,031,071	0	G	A	100.0%	109.6 / NA	35	P416L (CCG→CTG)	dcm	DNA cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

CGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGAT  >  NC_000913/2030936‑2031188                                                                                                                                        |                                                                                                                      cGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAgcg                                                                                                                    <  1:872127/139‑1 (MQ=255)       gCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCaa                                                                                                              <  1:423479/139‑1 (MQ=255)                    tGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAgcg                                                                                                 <  2:167223/139‑1 (MQ=255)                      aTTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAgcgcg                                                                                               <  1:1514559/139‑1 (MQ=255)                            ggCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCaaaa                                                                                                            <  2:868401/120‑1 (MQ=255)                            ggCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCaaaa                                                                                                            >  1:868401/1‑120 (MQ=255)                                    cAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGc                                                                                 <  1:1341981/139‑1 (MQ=255)                                                cgcgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGt                                                                     >  1:1914144/1‑139 (MQ=255)                                                 gcgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTg                                                                    <  2:1390815/139‑1 (MQ=255)                                                 gcgGCAAAGACCGGCACGACCACAGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAGTACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTg                                                                    >  1:752430/1‑139 (MQ=255)                                                     cAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACg                                                                <  1:572577/139‑1 (MQ=255)                                                                              aCCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCa                                       >  2:681259/1‑139 (MQ=255)                                                                               ccGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCaa                                      <  2:1914144/139‑1 (MQ=255)                                                                                       gCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTtc                              >  1:1910367/1‑139 (MQ=255)                                                                                               ccTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAgcg                                                                                                 <  1:1716575/64‑1 (MQ=255)                                                                                               ccTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAgcg                                                                                                 >  2:1716575/1‑64 (MQ=255)                                                                                                 tGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTgg                    >  1:470734/1‑139 (MQ=255)                                                                                                   aGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGcc                  >  1:1387887/1‑139 (MQ=255)                                                                                                       tCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCg                                                                                     >  1:417059/1‑68 (MQ=255)                                                                                                       tCCGAAACCGGAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCg                                                                                     <  2:417059/68‑1 (MQ=255)                                                                                                               cGGAATACGGAATTTCGCTTCTCCCAACGATTCAAAAGCCATTAAGGGCGCGCATTACCGAGGCGTTAACCGACGGGGAAGAGGTTGCTGGTTCAGCGGAGCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCcg      >  1:28229/1‑139 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCa                                                                                           >  1:399240/1‑52 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCa                                                                                           <  2:399240/52‑1 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCGCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCcgcg    <  1:333063/139‑1 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa                                                                  >  2:699317/1‑77 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa                                                                  <  1:699317/77‑1 (MQ=255)                                                                                                                 gAATACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCcgcg    >  1:230139/1‑139 (MQ=255)                                                                                                                   aTACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGcgat  <  2:470734/139‑1 (MQ=255)                                                                                                                   aTACGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGcgat  <  2:317758/139‑1 (MQ=255)                                                                                                                      cGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTccc                                                                                      >  1:1152718/1‑52 (MQ=255)                                                                                                                      cGGAATTTCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTccc                                                                                      <  2:1152718/52‑1 (MQ=255)                                                                                                                            ttCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTccc                                                                                      <  2:1401861/46‑1 (MQ=255)                                                                                                                            ttCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTccc                                                                                      >  1:1401861/1‑46 (MQ=255)                                                                                                                            ttCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAg                                                 <  2:1737760/83‑1 (MQ=255)                                                                                                                            ttCGCTTCTCCCAGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAg                                                 >  1:1737760/1‑83 (MQ=255)                                                                                                                                        |                                                                                                                      CGCAACGCCACCGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGATGTTGCTGATTCAGCGGATCGTCAAAGTCTTTCTCACCCGTGGCCATATCCCAGCCGCGAT  >  NC_000913/2030936‑2031188

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