Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,063 A→G 41.5% F485F (TTT→TTC 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,0630AG41.5% 19.3 / 46.5 65F485F (TTT→TTCcaiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base G (15/12);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCT  >  minE/33999‑34111
                                                                |                                                
cTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAAtt                                              >  1:1035703/1‑69 (MQ=255)
cTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAAtt                                              >  1:2199895/1‑69 (MQ=255)
     cTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2563313/65‑1 (MQ=255)
     cTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2549786/65‑1 (MQ=255)
      tGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2579541/64‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:213673/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2656317/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2586/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2537328/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2016752/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:302927/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:1664055/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:130284/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:3086953/65‑1 (MQ=255)
               gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGAAGAAGAATTCCTCTTCGCTc                                   <  1:1963413/65‑1 (MQ=255)
                aTAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:251102/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2929916/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2281558/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2086864/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2653237/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:1273347/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:1793624/64‑1 (MQ=255)
                  aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2717591/64‑1 (MQ=255)
                     agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:826738/1‑69 (MQ=255)
                     agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:851926/1‑69 (MQ=255)
                     agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:2828453/1‑69 (MQ=255)
                     agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:847984/1‑69 (MQ=255)
                        gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:2631179/69‑1 (MQ=255)
                        gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:1345692/69‑1 (MQ=255)
                        gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:218100/69‑1 (MQ=255)
                        gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:96204/69‑1 (MQ=255)
                         cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:262235/68‑1 (MQ=255)
                            tttAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2591434/52‑1 (MQ=255)
                                       ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:2684036/68‑1 (MQ=255)
                                       ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:1596978/68‑1 (MQ=255)
                                       ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:2166235/68‑1 (MQ=255)
                                       ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:372740/68‑1 (MQ=255)
                                           atTCTGTTCACATAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1664905/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGGTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2860974/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGt                              >  1:174789/1‑42 (MQ=25)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2766186/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1036620/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2960105/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2547904/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2477517/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2215805/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2144981/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1928805/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1278273/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1588034/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1541815/1‑69 (MQ=255)
                                           atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:145093/1‑69 (MQ=255)
                                            ttttttTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:1445782/1‑69 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:990473/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:2663208/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:1486895/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:557297/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:420258/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:2895117/1‑49 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:547/1‑69 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:446581/1‑69 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:17991/1‑69 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:984811/1‑69 (MQ=255)
                                            ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACATTCATTCAGCACCACAAATGCt  >  1:748321/1‑69 (MQ=255)
                                              ttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        >  1:730160/1‑61 (MQ=255)
                                                                |                                                
CTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCT  >  minE/33999‑34111

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