Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,662,890	A→G	L145P (CTG→CCG)	trmJ ←	tRNA Cm32/Um32 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,662,890	0	A	G	100.0%	117.8 / NA	34	L145P (CTG→CCG)	trmJ	tRNA Cm32/Um32 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/16);  total (18/16)

CATAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACG  >  NC_000913/2662762‑2663024                                                                                                                                 |                                                                                                                                       cATAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAg                                                                                                                             >  1:586631/1‑140 (MQ=255) cATAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAg                                                                                                                             >  1:655961/1‑140 (MQ=255)  aTAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGc                                                                                                                            <  1:415672/140‑1 (MQ=255)   tAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCt                                                                                                                           <  1:848257/140‑1 (MQ=255)     aaaaCGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGt                                                                                                                         >  2:441715/1‑140 (MQ=255)                      cGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCa                                                                                                        >  1:779277/1‑140 (MQ=255)                        aCCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAAt                                                                                                      <  2:816201/140‑1 (MQ=255)                          cAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATcg                                                                                                    <  1:220344/140‑1 (MQ=255)                          cAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATcg                                                                                                    <  2:586631/140‑1 (MQ=255)                                   cGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGAt                                                                                           <  2:760460/140‑1 (MQ=255)                                          tctTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATT‑CGGGTTAGCCGCAATCGCGACa                                                                                               >  1:164051/1‑128 (MQ=255)                                          tctTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATT‑CGGGTTAGCCGCAATCGCGACa                                                                                               <  2:164051/128‑1 (MQ=255)                                                tGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTg                                                                              >  1:654970/1‑140 (MQ=255)                                                   tCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCaa                                                                           <  2:779277/140‑1 (MQ=255)                                                     gACCTGCCCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAct                                                                         <  1:751946/140‑1 (MQ=255)                                                      aCCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAActc                                                                        >  1:950837/1‑140 (MQ=255)                                                                              gccagccaGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTc                                                >  1:721474/1‑140 (MQ=255)                                                                               ccagccagGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcg                                               <  1:441715/140‑1 (MQ=255)                                                                                 agccagGCCATGTGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTcgcg                                             >  2:481786/1‑140 (MQ=255)                                                                                  gccagGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGa                                            <  2:655961/140‑1 (MQ=255)                                                                                    cagGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAcc                                          >  2:1008995/1‑140 (MQ=255)                                                                                      ggCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAcg                                                     <  2:27614/127‑1 (MQ=255)                                                                                      ggCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAcg                                                     >  1:27614/1‑127 (MQ=255)                                                                                                  ttCATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCa                            <  1:481786/140‑1 (MQ=255)                                                                                                    cATAGGCAATAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCACGACCAAACACCAGCGCCAcc                          <  2:721474/140‑1 (MQ=255)                                                                                                            aTAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAtt                  >  2:673143/1‑140 (MQ=255)                                                                                                            aTAACCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAtt                  >  2:757894/1‑140 (MQ=255)                                                                                                              aaCCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                              >  2:388688/1‑126 (MQ=255)                                                                                                              aaCCTGAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgcgc                              <  1:388688/126‑1 (MQ=255)                                                                                                                   gAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATTATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCt           >  1:128651/1‑140 (MQ=255)                                                                                                                   gAACCGCCATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCt           >  2:669692/1‑140 (MQ=255)                                                                                                                         ccATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCg     >  1:195370/1‑140 (MQ=255)                                                                                                                          cATCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGa    >  2:711534/1‑140 (MQ=255)                                                                                                                            tCGCCGGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACg  <  1:757894/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       CATAAAAACGCTCCAGATCATCGACCAGCGGATACGGCGTCTCTTCATGCTCGACCTGCTCGCCGTTTTCCTGAGTTGCCAGCCAGGCCATGCGCACTTCATAGGCAATAACCTGAACCGCCATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCAGCGACG  >  NC_000913/2662762‑2663024

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