Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,469 A→G 41.6% V350A (GTG→GCG)  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,4690AG41.6% 13.3 / 33.4 48V350A (GTG→GCG) caiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (7/21);  new base G (11/9);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 6.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

ACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGG  >  minE/34402‑34532
                                                                   |                                                               
aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:184136/68‑1 (MQ=255)
aCCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                                                                 <  1:839092/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:1352944/69‑1 (MQ=255)
   gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:3116157/69‑1 (MQ=255)
   gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCg                                                             <  1:279620/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:3103770/1‑69 (MQ=255)
                             tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:1673844/1‑69 (MQ=255)
                             tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:688466/1‑69 (MQ=255)
                             tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt                                   >  1:2005865/1‑69 (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)
                                 aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:2641760/69‑1 (MQ=255)
                                 aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:534624/69‑1 (MQ=255)
                                 aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:163544/69‑1 (MQ=255)
                                 aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc                               <  1:3303440/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:981021/1‑69 (MQ=255)
                                       gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:3307254/1‑69 (MQ=255)
                                       gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:39801/1‑69 (MQ=255)
                                       gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:617291/1‑69 (MQ=255)
                                       gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:2035523/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)
                                       gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg                         >  1:3080975/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:2233556/69‑1 (MQ=255)
                                                         gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2113037/69‑1 (MQ=255)
                                                         gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2704419/69‑1 (MQ=255)
                                                         gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:2087627/69‑1 (MQ=255)
                                                         gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:46654/69‑1 (MQ=255)
                                                         gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac       <  1:1216865/69‑1 (MQ=255)
                                                              aCCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATgg  >  1:1046631/1‑69 (MQ=255)
                                                                   |                                                               
ACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGG  >  minE/34402‑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: