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

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

GTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGTTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTA  >  NC_000913/2031687‑2031946                                                                                                                                   |                                                                                                                                  gTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTc                                                                                                                                 >  2:149154/1‑133 (MQ=255)  tACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCa                                                                                                                                >  2:455312/1‑133 (MQ=255)       caTCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAAt                                                                                                                           <  2:552872/133‑1 (MQ=255)        aTCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAAtg                                                                                                                          >  1:91856/1‑133 (MQ=255)        aTCAAAAAACAGCGCGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAAtg                                                                                                                          >  2:308117/1‑133 (MQ=255)                                  cGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCcgccgccgc                                                                                                >  2:36804/1‑133 (MQ=255)                                       gCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTACGCCGCCGCCtcat                                                                                           >  1:339677/1‑133 (MQ=255)                                       gCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGATGAAGAATATGTGCCGCTGCCGCCtcat                                                                                           >  2:194077/1‑133 (MQ=255)                                               gTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACg                                                                                   >  1:536847/1‑133 (MQ=255)                                                 gcgcCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGcc                                                                                 <  2:536847/133‑1 (MQ=255)                                                 gcgcCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGcc                                                                                 <  2:22446/133‑1 (MQ=255)                                                        cccgAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTa                                                                          <  1:308117/133‑1 (MQ=255)                                                             gcgagTTCTTTTTCGGTACGCCAGGCAGCGAAAATGGCTGATAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCt                                                                     <  2:262520/133‑1 (MQ=255)                                                             gcgagTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCt                                                                     <  1:36804/133‑1 (MQ=255)                                                             gcgagTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCt                                                                     <  2:486827/133‑1 (MQ=255)                                                                                                       aaGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGc                           <  2:12783/133‑1 (MQ=255)                                                                                                       aaGGGAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGc                           <  1:149154/133‑1 (MQ=255)                                                                                                         gggAAACCGGCCAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCg                         >  2:432674/1‑133 (MQ=255)                                                                                                                   ccAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCaata               >  2:420274/1‑133 (MQ=255)                                                                                                                   ccAGTAAAACATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCaata               >  1:89615/1‑133 (MQ=255)                                                                                                                        aaaaCATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGt          <  2:91856/133‑1 (MQ=255)                                                                                                                         aaaCATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGtt         <  1:432674/133‑1 (MQ=255)                                                                                                                           aCATCGTGCTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTgg       >  2:438006/1‑133 (MQ=255)                                                                                                                              tCGTGCTCGGGAATGTGTGGACGAATATGTGCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCtt    <  1:250834/133‑1 (MQ=255)                                                                                                                                gtgcTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTta  <  1:420274/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  GTACCACATCAAAAAACAGCGTGCCCTGGGTATCGCAGGCAAAACCGTGCGCCCGCCCGAGCGAGTTCTTTTTCGATACGCCAGCCAGCGAAAATGGCTGACAAGGGAAACCGGCCAGTAAAACATCGTGTTCAGGAATGTGTTGACGAATATGTTCCGCCGCCGCCTCATCACTCACGCCTTCTTTATGGCTGAGGGTGATGTCGCGGATATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTA  >  NC_000913/2031687‑2031946

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