Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,102,762	T→C	T31A (ACC→GCC)	trmI ←	tRNA m(7)G46 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,102,762	0	T	C	100.0%	110.8 / NA	32	T31A (ACC→GCC)	trmI	tRNA m(7)G46 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/15);  total (17/15)

CCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAAACGGGCATT  >  NC_000913/3102628‑3102894                                                                                                                                       |                                                                                                                                     ccccATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTc                                                                                                                                 <  2:588102/140‑1 (MQ=255)       gCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCt                                                                                                                           >  1:740655/1‑140 (MQ=255)           aaaCCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTggcg                                                                                                                       <  2:385886/140‑1 (MQ=255)           aaaCCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTggcg                                                                                                                       >  1:601004/1‑140 (MQ=255)                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGt                                                                                                                                  >  1:207303/1‑123 (MQ=255)                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGt                                                                                                                                  <  2:207303/123‑1 (MQ=255)                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGcac                                                                                                                 >  2:418888/1‑140 (MQ=255)                  tctcAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGcaca                                                                                                                <  2:473661/140‑1 (MQ=255)                  tctcAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGcaca                                                                                                                <  2:24838/140‑1 (MQ=255)                       aGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACt                                                                                                           >  2:642605/1‑140 (MQ=255)                             aCCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGAt                                                                                                     >  1:216090/1‑140 (MQ=255)                             aCCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGAt                                                                                                     >  1:587366/1‑140 (MQ=255)                                      tCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAg                                                                                            >  1:819250/1‑140 (MQ=255)                                                  aGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGtt                                                                                <  2:89727/140‑1 (MQ=255)                                                   gcgcGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGttt                                                                               <  1:418888/140‑1 (MQ=255)                                                    cgcgGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGtttt                                                                              >  1:579733/1‑140 (MQ=255)                                                      cgGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTtca                                                                            >  2:543862/1‑140 (MQ=255)                                                                     tatCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGaa                                                             >  1:294379/1‑140 (MQ=255)                                                                        cTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATg                                                          >  2:222417/1‑140 (MQ=255)                                                                           cGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGAcg                                                       <  2:740655/140‑1 (MQ=255)                                                                                aaCTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGTAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGtt                                                  <  2:587366/140‑1 (MQ=255)                                                                                aaCTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGtt                                                  >  2:550177/1‑140 (MQ=255)                                                                                                  ggCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTc                                <  1:26524/140‑1 (MQ=255)                                                                                                    ccAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGc                              <  1:222417/140‑1 (MQ=255)                                                                                                            tttCCAGCGCATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGaaa                      <  2:294379/140‑1 (MQ=255)                                                                                                                  gcgcATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGc                              >  2:914786/1‑126 (MQ=255)                                                                                                                  gcgcATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGc                              <  1:914786/126‑1 (MQ=255)                                                                                                                  gcgcATGTTCCTGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTc                >  1:276659/1‑140 (MQ=255)                                                                                                                             tGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGaaa                      >  2:570258/1‑123 (MQ=255)                                                                                                                             tGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGaaa                      <  1:570258/123‑1 (MQ=255)                                                                                                                             tGGCCTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAAACGGGc     <  2:216090/140‑1 (MQ=255)                                                                                                                                ccTTTGGCCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAAACGGGCAtt  >  1:540178/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                     CCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAAACGGGCATT  >  NC_000913/3102628‑3102894

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