Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,399,978	T→C	T87A (ACG→GCG)	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	0	T	C	100.0%	170.8 / NA	47	T87A (ACG→GCG)	ogt	methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/26);  total (21/26)

GATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGC  >  NC_000913/1399843‑1400114                                                                                                                                        |                                                                                                                                         gATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCcccc                                                                                                                                      <  1:786643/140‑1 (MQ=255)  aTTCGATCCGTTTGCCGTACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGcccccc                                                                                                                                     >  2:299127/1‑140 (MQ=255)     cGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCg                                                                                                                                  <  2:559199/140‑1 (MQ=255)       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTa                                                                                                                                >  1:917502/1‑140 (MQ=255)       aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGCAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTa                                                                                                                                <  1:953995/140‑1 (MQ=255)         ccGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGc                                                                                                                              <  2:469538/140‑1 (MQ=255)               gCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTggg                                                                                                                        <  1:900017/140‑1 (MQ=255)                ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                       <  1:1001202/140‑1 (MQ=255)                ccGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGa                                                                                                                       <  2:60007/140‑1 (MQ=255)                 cGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGaa                                                                                                                      <  2:917502/140‑1 (MQ=255)                              cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCaataat                                                                                                         >  1:599654/1‑140 (MQ=255)                                             cGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGc                                                                                          <  2:834084/140‑1 (MQ=255)                                              ggCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                                         <  1:428811/140‑1 (MQ=255)                                                cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCaaa                                                                                       >  2:100674/1‑140 (MQ=255)                                                 ccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCaaaa                                                                                      <  2:300688/140‑1 (MQ=255)                                                        cTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCAc                                                                               >  2:127945/1‑140 (MQ=255)                                                               aGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                        >  2:677405/1‑140 (MQ=255)                                                               aGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                                        <  2:590547/140‑1 (MQ=255)                                                                   ggCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCg                                                                    >  2:841292/1‑140 (MQ=255)                                                                        tAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaa                                                               <  1:100674/140‑1 (MQ=255)                                                                         aaTGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTaaa                                                              >  2:611175/1‑140 (MQ=255)                                                                               ttACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCt                                                        <  1:299127/140‑1 (MQ=255)                                                                                  ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                     >  1:842723/1‑140 (MQ=255)                                                                                  ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                     <  1:611175/140‑1 (MQ=255)                                                                                  ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                                     >  1:129715/1‑140 (MQ=255)                                                                                      cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                                 >  2:139182/1‑140 (MQ=255)                                                                                       ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggt                                                >  2:648857/1‑140 (MQ=255)                                                                                          cAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGc                                             >  1:52050/1‑140 (MQ=255)                                                                                                 tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                                      <  2:842723/140‑1 (MQ=255)                                                                                                 tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                                      <  2:416036/140‑1 (MQ=255)                                                                                                 tAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                                      <  2:129715/140‑1 (MQ=255)                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTCCCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                   <  1:806629/140‑1 (MQ=255)                                                                                                    tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                   <  2:773459/140‑1 (MQ=255)                                                                                                     gCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCt                                  >  1:148795/1‑140 (MQ=255)                                                                                                       gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCa                                >  2:783977/1‑140 (MQ=255)                                                                                                        tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCAt                               >  1:72559/1‑140 (MQ=255)                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATa                              >  2:139192/1‑140 (MQ=255)                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATa                              >  1:759571/1‑140 (MQ=255)                                                                                                         aGTGTTTTCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATa                              >  2:782915/1‑140 (MQ=255)                                                                                                              tttCCAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCctt                         <  2:544544/140‑1 (MQ=255)                                                                                                                  cAGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCttt                     >  2:51107/1‑140 (MQ=255)                                                                                                                   aGACTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGCATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTg                    >  2:301116/1‑140 (MQ=255)                                                                                                                     aCTTCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCg                  >  2:16889/1‑140 (MQ=255)                                                                                                                      cTTCGCGCTGAAAGGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGa                 <  2:1063991/140‑1 (MQ=255)                                                                                                                       ttCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt                >  2:985407/1‑140 (MQ=255)                                                                                                                       ttCGCGCTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt                <  2:818459/140‑1 (MQ=255)                                                                                                                             cTGAAATGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGa          <  1:165156/140‑1 (MQ=255)                                                                                                                                  aTGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc     <  1:139192/140‑1 (MQ=255)                                                                                                                                  aTGGCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc     <  2:1074256/140‑1 (MQ=255)                                                                                                                                    ggCGCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcag   <  1:16889/140‑1 (MQ=255)                                                                                                                                     gcgcCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCcagc  <  1:677405/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCCAGC  >  NC_000913/1399843‑1400114

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: ‑