Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,399,985	(C)7→8	coding (252/516 nt)	ogt ←	methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,399,978	1	.	C	100.0%	77.1 / NA	22	T87A (ACG→GCG)	ogt	methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (13/9);  total (13/9)

GATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGG  >  NC_000913/1399843‑1400105                                                                                                                                          |                                                                                                                                 gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑cccc                                                                                                                             <  2:273423/140‑1 (MQ=255)  aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑ccccc                                                                                                                            <  1:169335/140‑1 (MQ=255)    tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑ccccccc                                                                                                                          >  2:515215/1‑140 (MQ=255)     cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑cc                                                                                                                               <  2:183256/134‑1 (MQ=255)     cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑cc                                                                                                                               >  1:183256/1‑134 (MQ=255)            tttGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAg                                                                                                                   >  2:123536/1‑140 (MQ=255)                        aCGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                                       <  2:50101/140‑1 (MQ=255)                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCaataa                                                                                                 >  1:233262/1‑140 (MQ=255)                                     cgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                    <  2:328571/130‑1 (MQ=255)                                     cgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                    >  1:328571/1‑130 (MQ=255)                                                           aGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGa                                                                    <  2:233262/140‑1 (MQ=255)                                                                 ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAata                                                                           <  1:468524/127‑1 (MQ=255)                                                                 ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAata                                                                           >  2:468524/1‑127 (MQ=255)                                                                    gCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCg                                                           <  1:383497/140‑1 (MQ=255)                                                                          aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaaa                                                     >  1:29347/1‑140 (MQ=255)                                                                          aTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaaa                                                     <  2:29347/140‑1 (MQ=255)                                                                                     gCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGAtt                                          >  2:433699/1‑140 (MQ=255)                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                         >  2:194395/1‑140 (MQ=255)                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                        >  1:30564/1‑140 (MQ=255)                                                                                                  aGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                             <  1:313472/140‑1 (MQ=255)                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TCCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgc                            >  2:67703/1‑140 (MQ=255)                                                                                                     gCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                             >  1:213338/1‑137 (MQ=255)                                                                                                     gCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                             <  2:213338/137‑1 (MQ=255)                                                                                                             ttttCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTGCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc                  >  2:376115/1‑140 (MQ=255)                                                                                                             ttttCCAGACTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc                  >  1:296824/1‑140 (MQ=255)                                                                                                                 ccAGACTTCGCGCTGAAATGGCGCC‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCcttct              <  2:296824/140‑1 (MQ=255)                                                                                                                     aCTTCGCGCTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGc          >  2:164087/1‑140 (MQ=255)                                                                                                                             cTGAAATGGCG‑TC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATgg  <  1:164087/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                 GATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCG‑T‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGG  >  NC_000913/1399843‑1400105

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