Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 816,410 C→T 60.6% S161S (AGC→AGT hemA → glutamyl tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE816,4100CT60.6% 12.6 / 23.3 38S161S (AGC→AGThemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (4/11);  new base T (8/15);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTCAC  >  minE/816353‑816478
                                                         |                                                                    
ccAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctgt                                                           >  1:1769583/1‑69 (MQ=255)
ccAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctgt                                                           >  1:2698738/1‑69 (MQ=255)
ccAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctgt                                                           >  1:2994250/1‑69 (MQ=255)
ccAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctg                                                            >  1:3031373/1‑68 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:935036/69‑1 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:762103/69‑1 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:334684/69‑1 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:2032470/69‑1 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:1049301/69‑1 (MQ=255)
      aTCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:974425/69‑1 (MQ=255)
       tCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:967409/68‑1 (MQ=255)
       tCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:1088819/68‑1 (MQ=255)
       tCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:1656302/68‑1 (MQ=255)
       tCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:3204952/68‑1 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:513857/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:1027471/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:2926061/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:2720268/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:2590325/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:1012102/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:2031841/1‑69 (MQ=255)
           tctcTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGTGCTGTGTCTGTCGCTTTTGCgg                                                >  1:2586350/1‑69 (MQ=255)
                        aCGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCttt                                                     <  1:1204009/51‑1 (MQ=255)
                                                 ggTGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGctctct          <  1:3393589/69‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2662888/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2734877/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:1049542/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:767955/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:1767550/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:388771/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:3284860/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2044872/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:240288/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2944757/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2548076/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:2871988/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGTCGCTTTTGCGGCTTGTACGATGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:347408/68‑1 (MQ=255)
                                                   tGCCAGTGCTGTGTCTGGCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTa         <  1:930722/68‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:329934/69‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:2032803/69‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:3185206/69‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:1467564/69‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:765769/69‑1 (MQ=255)
                                                         cgcTGTGTCTGTCGCTATTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTcac  <  1:2761382/69‑1 (MQ=255)
                                                         |                                                                    
CCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATCTTTGAATCGCTCTCTACGGTCAC  >  minE/816353‑816478

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