Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,516,656 C→T A102T (GCT→ACT)  trmH ← tRNA (Guanosine‑2'‑O‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,516,6560CT100.0% 137.1 / NA 42A102T (GCT→ACT) trmHtRNA (Guanosine‑2'‑O‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/30);  total (12/30)

AAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCA  >  minE/2516610‑2516717
                                              |                                                             
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAGATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:1530038/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:87756/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:851593/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:638459/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:502249/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:496753/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:1622013/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:3185965/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:288423/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:1971525/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:243188/1‑68 (MQ=255)
aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGc                                          >  1:2225114/1‑68 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:1165058/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:963953/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:3214175/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:2576710/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:1151505/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:2058534/71‑1 (MQ=255)
      cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:1171011/71‑1 (MQ=255)
       aGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:1942863/70‑1 (MQ=255)
       aGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:3003286/70‑1 (MQ=255)
             ggCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:1470193/64‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:1873363/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:913375/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:813159/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:658885/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:1427599/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:630046/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:506682/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:1560866/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:2185346/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:192392/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:3067464/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:2717260/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:2508186/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:2351812/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:2086784/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCGGGCCTTTGAGATGAGCGACGGc   <  1:2961677/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGGGTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:246058/71‑1 (MQ=255)
                                    aaaTCGACAGTGTTATCAGAAAGATGGCTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc   <  1:1344893/71‑1 (MQ=255)
                                     aaTCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCa  <  1:1994234/71‑1 (MQ=255)
                                        cGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                 <  1:689327/37‑1 (MQ=255)
                                              |                                                             
AAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCA  >  minE/2516610‑2516717

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: