Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,102,340	(C)6→5	coding (513/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,335	0	C	.	100.0%	200.1 / NA	47	coding (518/720 nt)	trmI	tRNA m(7)G46 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (24/23);  total (24/23)

TCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAA  >  NC_000913/3102201‑3102464                                                                                                                                       |                                                                                                                                  tCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGccccca                                                                                                                              >  2:96886/1‑139 (MQ=255)                  gCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCttttt                                                                                                            >  1:122632/1‑140 (MQ=255)                     ggACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTAc                                                                                                          <  1:138605/139‑1 (MQ=255)                     ggACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTAc                                                                                                          >  2:138605/1‑139 (MQ=255)                         gcgGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTCTTTACCAGtt                                                                                                     >  1:489711/1‑140 (MQ=255)                            gTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCgg                                                                                                  <  2:320208/140‑1 (MQ=255)                                 tAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                                                                                             <  1:924004/140‑1 (MQ=255)                                 tAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                                                                                             <  1:760720/140‑1 (MQ=255)                                     cATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCa                                                                                         >  2:270385/1‑140 (MQ=255)                                      aTTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCAc                                                                                        <  2:342412/140‑1 (MQ=255)                                       ttGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCAcc                                                                                       >  1:167547/1‑140 (MQ=255)                                        tGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCt                                                                                      >  2:52666/1‑140 (MQ=255)                                         gCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                                                                                             <  1:573935/132‑1 (MQ=255)                                         gCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                                                                                             >  2:573935/1‑132 (MQ=255)                                           tctctGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGaa                                                                                   <  1:649471/140‑1 (MQ=255)                                                gACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATa                                                                              <  1:803417/140‑1 (MQ=255)                                                gACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATa                                                                              <  2:260763/140‑1 (MQ=255)                                                    ggTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGc                                                                          <  1:52666/140‑1 (MQ=255)                                                    ggTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTCTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGc                                                                          <  2:489711/140‑1 (MQ=255)                                                     gTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACg                                                                            <  1:215198/137‑1 (MQ=255)                                                     gTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTCACCAGTTCGGCAAACGGCACCTGAACGATACg                                                                            >  2:215198/1‑137 (MQ=255)                                                               ccGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATgg                                                               >  2:993545/1‑140 (MQ=255)                                                                cGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgc                                                              >  2:962673/1‑140 (MQ=255)                                                                  tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATAc                                                                             <  2:55720/123‑1 (MQ=255)                                                                  tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATAc                                                                             >  1:55720/1‑123 (MQ=255)                                                                  tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgcgc                                                            <  1:96886/140‑1 (MQ=255)                                                                  tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGgcgc                                                            >  2:719267/1‑140 (MQ=255)                                                                           gaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGc                                                   >  1:933204/1‑140 (MQ=255)                                                                            aCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGcc                                                  >  2:964995/1‑140 (MQ=255)                                                                              aTCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCa                                                 <  1:500133/139‑1 (MQ=255)                                                                              aTCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCa                                                 >  2:500133/1‑139 (MQ=255)                                                                              aTCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCAc                                                >  1:315891/1‑140 (MQ=255)                                                                                cACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCAACTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACgg                                              <  1:270385/140‑1 (MQ=255)                                                                                 aCTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGACGCGCTTTGTGCCACggg                                             >  1:964580/1‑140 (MQ=255)                                                                                     cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAg                                         >  2:1085704/1‑140 (MQ=255)                                                                                        gCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGa                                      >  1:468234/1‑140 (MQ=255)                                                                                            atGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGaaaaa                                  <  2:778051/140‑1 (MQ=255)                                                                                                 ccGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCt                             <  1:962673/140‑1 (MQ=255)                                                                                                    cATAAGGTTCCCAGTCGGCCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCa                          >  1:981082/1‑140 (MQ=255)                                                                                                      tAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCAcc                        <  2:468234/140‑1 (MQ=255)                                                                                                      tAAGGTTCCCAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCAcc                        <  2:167547/140‑1 (MQ=255)                                                                                                             cccAGTCGGTCGCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCa                 <  2:705115/140‑1 (MQ=255)                                                                                                                     gtcgCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATgcgc                  >  1:660554/1‑131 (MQ=255)                                                                                                                     gtcgCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATgcgc                  <  2:660554/131‑1 (MQ=255)                                                                                                                     gtcgCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTg         >  1:808292/1‑140 (MQ=255)                                                                                                                     gtcgCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTg         >  1:821740/1‑140 (MQ=255)                                                                                                                     gtcgCCATATGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTg         >  1:891212/1‑140 (MQ=255)                                                                                                                            tatGGAATACG‑CCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGaa  <  1:870340/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  TCGTCACCGGACGTGATGCCGGACGCGGTACGTAATCATTGCTCTCTGACAGGTTTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTTTATTATGGCGCGCTTTGTGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAA  >  NC_000913/3102201‑3102464

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