Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,492 T→C 48.3% R342R (CGA→CGG caiC ← predicted crotonobetaine CoA ligase:carnitine CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE34,4920TC48.3% ‑5.5 / 49.6 56R342R (CGA→CGGcaiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (9/20);  new base C (18/9);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-02

TGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAA  >  minE/34431‑34560
                                                             |                                                                    
tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                                               >  1:3103770/1‑69 (MQ=255)
tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                                               >  1:2005865/1‑69 (MQ=255)
tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                                               >  1:688466/1‑69 (MQ=255)
tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                                               >  1:1673844/1‑69 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:913686/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:2156192/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:2641760/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:3303440/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:534624/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:163544/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:1238546/69‑1 (MQ=255)
    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:698033/69‑1 (MQ=255)
     tcgtcgCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                                                           <  1:738024/68‑1 (MQ=255)
      cgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg                                                         <  1:1791720/69‑1 (MQ=255)
      cgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg                                                         <  1:2942670/69‑1 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCa                                                          >  1:1441268/1‑64 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAgggct                                                     >  1:2083482/1‑66 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:617291/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:39801/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:3307254/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:3080975/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:2724977/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:981021/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:1742977/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:2035523/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                                                     >  1:2054424/1‑69 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                                                     <  1:2150487/69‑1 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                                                     <  1:858926/69‑1 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                                                     <  1:754091/69‑1 (MQ=255)
          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                                                     <  1:1685169/69‑1 (MQ=255)
                     cTTCGTAGCAAAACCCCACCCGACCAATCGACGTCCAGCGTCGTTTATCGCCAGGACGATCGCCGAt                                            <  1:2875741/67‑1 (MQ=255)
                     cTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGAt                                            <  1:2991562/67‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:1216865/69‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:46654/69‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:2087627/69‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:2113037/69‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:2704419/69‑1 (MQ=255)
                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac                                   <  1:2233556/69‑1 (MQ=255)
                                 aCCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATgg                              >  1:1046631/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:1763214/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:525993/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:74124/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:3120966/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:2909198/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:2127314/1‑69 (MQ=255)
                                        ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt                       >  1:2492063/1‑69 (MQ=255)
                                             cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt                                >  1:3317422/1‑55 (MQ=255)
                                             cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt                                >  1:3228199/1‑55 (MQ=255)
                                             cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt                                >  1:2322124/1‑55 (MQ=255)
                                             cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt                                >  1:2145839/1‑55 (MQ=255)
                                              aaTCGACGGCCAGCGTCGTTTATCGCCAGGCCGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                 <  1:1864592/69‑1 (MQ=255)
                                              aaTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                 <  1:892129/69‑1 (MQ=255)
                                              aaTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                 <  1:940830/69‑1 (MQ=255)
                                               aTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc                 <  1:1141228/68‑1 (MQ=255)
                                                  gACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATaa             <  1:1710907/69‑1 (MQ=255)
                                                  gACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATaa             <  1:1154894/69‑1 (MQ=255)
                                                             tcgtTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCaa  >  1:1481386/1‑69 (MQ=255)
                                                             |                                                                    
TGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAA  >  minE/34431‑34560

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: