Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,102,751	Δ1 bp	coding (102/720 nt)	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,750	0	C	.	100.0%	97.4 / NA	24	coding (103/720 nt)	trmI	tRNA m(7)G46 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (13/11);  total (13/11)

TTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAA  >  NC_000913/3102612‑3102885                                                                                                                                           |                                                                                                                                        ttGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTct                                                                                                                                        >  1:580501/1‑139 (MQ=255)  tGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTctg                                                                                                                                       >  1:173219/1‑138 (MQ=255)          gCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTgg                                                                                                                              >  1:278994/1‑140 (MQ=255)                           aaaCCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGGAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTggcgg                                                                                                             >  2:55344/1‑140 (MQ=255)                                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGgcgc                                                                                                          >  2:391224/1‑137 (MQ=255)                                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGgcgc                                                                                                          <  1:391224/137‑1 (MQ=255)                                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGcaca                                                                                                       >  2:614817/1‑140 (MQ=255)                                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGcaca                                                                                                       >  1:493852/1‑140 (MQ=255)                                 aTCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGcaca                                                                                                       >  1:378290/1‑140 (MQ=255)                                  tctcAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACaa                                                                                                      <  2:580501/140‑1 (MQ=255)                                       aGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTa                                                                                                 >  2:59719/1‑140 (MQ=255)                                                  cgcTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGgcg                                                                                      >  1:234070/1‑140 (MQ=255)                                                         cGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTggg                                                                               >  2:40174/1‑140 (MQ=255)                                                                      cgGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTtcat                                                                  >  1:245917/1‑140 (MQ=255)                                                                          gAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCaaa                                                              <  2:245917/140‑1 (MQ=255)                                                                             aTCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGTCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTc                                                           <  1:55344/140‑1 (MQ=255)                                                                                   cATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGa                                                     >  1:708214/1‑140 (MQ=255)                                                                                               gAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGtt                                         <  1:83612/140‑1 (MQ=255)                                                                                                aaCTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGttt                                        <  2:378290/140‑1 (MQ=255)                                                                                                      aCGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATa                                  <  2:234070/140‑1 (MQ=255)                                                                                                      aCGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATa                                  <  2:493852/140‑1 (MQ=255)                                                                                                               aCCGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAg                         <  2:278994/140‑1 (MQ=255)                                                                                                                 cGGCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTc                       <  1:161851/140‑1 (MQ=255)                                                                                                                   gCCAGTAGTTTTCCAGCGCATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGc                     >  1:251872/1‑140 (MQ=255)                                                                                                                                     cATGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGaa   >  1:230120/1‑140 (MQ=255)                                                                                                                                      aTGTT‑CTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGaaa  <  1:59719/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        TTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCTGAACTCAACGCCCATCACCGGCCAGTAGTTTTCCAGCGCATGTTCCTGGCCTTTGGTCAGTCGCCCCTGGCGGCGCACAAAACTACGGATACGGCGCAGTGGGCGGCCGTTTTCATCAAATTCCGGTGAAATGACGTCGTTTTTCATAAAGGTTTAGTCGCTTGTGAAAGTGTTCTGAAA  >  NC_000913/3102612‑3102885

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