Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,465 C→G 44.2% G351G (GGG→GGC 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,4650CG44.2% 17.3 / 64.6 61G351G (GGG→GGCcaiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (7/27);  new base G (11/16);  total (18/43)
Fisher's exact test for biased strand distribution p-value = 9.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGG  >  minE/34399‑34532
                                                                  |                                                                   
ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAccc                                                                   >  1:2010824/1‑69 (MQ=255)
ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAccc                                                                   >  1:1829472/1‑69 (MQ=255)
ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAccc                                                                   >  1:3379692/1‑69 (MQ=255)
ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAccc                                                                   >  1:3246369/1‑69 (MQ=255)
   aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:839092/68‑1 (MQ=255)
   aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:184136/68‑1 (MQ=255)
   aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:1692420/68‑1 (MQ=255)
   aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCAGCTTCGTAGCAAAACCCCa                                                                 <  1:2054549/68‑1 (MQ=255)
   aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGACCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:863281/68‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:316767/66‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:1883945/66‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:1413923/66‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:1146044/66‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:1133455/66‑1 (MQ=255)
      gagctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTACGTAGCAAAAGcccgcccg                                                             <  1:759879/66‑1 (MQ=255)
      gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:3116157/69‑1 (MQ=255)
      gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:279620/69‑1 (MQ=255)
      gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:1352944/69‑1 (MQ=255)
      gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:1266188/69‑1 (MQ=255)
       agctcACCAGCCGGGAGGGGACGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAGcccgcccg                                                             <  1:2985139/66‑1 (MQ=255)
                  cGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:1750914/53‑1 (MQ=255)
                   gggAGCGGGCGATTGTGATCGTCGCGGATCTCAGCTTCGTAGCAAAACCCCACCCg                                                             <  1:1707255/56‑1 (MQ=255)
                    ggAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGcc                                               <  1:2438480/69‑1 (MQ=255)
                     gAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGcc                                               <  1:1331444/68‑1 (MQ=255)
                       gCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAACCCCCACCCGACCAATCGACGGCCAGc                                            >  1:2631419/1‑69 (MQ=255)
                       gCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGc                                            >  1:2772182/1‑69 (MQ=255)
                                tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:1673844/1‑69 (MQ=255)
                                tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:2005865/1‑69 (MQ=255)
                                tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:3103770/1‑69 (MQ=255)
                                tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:688466/1‑69 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:163544/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:3303440/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:534624/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:2641760/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:698033/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:2156192/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:1238546/69‑1 (MQ=255)
                                    aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:913686/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:2035523/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:981021/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:3080975/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:2054424/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:2724977/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:1742977/1‑69 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                         <  1:754091/69‑1 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                         <  1:2150487/69‑1 (MQ=255)
                                          gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg                         <  1:858926/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:46654/69‑1 (MQ=255)
                                                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2704419/69‑1 (MQ=255)
                                                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2233556/69‑1 (MQ=255)
                                                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2113037/69‑1 (MQ=255)
                                                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:1216865/69‑1 (MQ=255)
                                                            gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2087627/69‑1 (MQ=255)
                                                                 aCCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATgg  >  1:1046631/1‑69 (MQ=255)
                                                                  |                                                                   
TTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGG  >  minE/34399‑34532

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: