Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,427 C→T 25.0% R364H (CGC→CAC)  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,4270CT25.0% 60.2 / 18.4 36R364H (CGC→CAC) caiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base T (2/7);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT  >  minE/34362‑34496
                                                                 |                                                                     
tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGAt                                                                    <  1:3302761/69‑1 (MQ=255)
tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGa                                                                     <  1:2037027/68‑1 (MQ=255)
  aaGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtgtg                                                                <  1:2150807/71‑1 (MQ=255)
  aaGATGGTGTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtgtg                                                                <  1:2769180/71‑1 (MQ=255)
         ttttCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGa                                                               <  1:1298080/65‑1 (MQ=255)
            tCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTcgc                                                        >  1:3873042/1‑69 (MQ=255)
              ccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGa                                                     <  1:1055756/70‑1 (MQ=255)
                  gTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGa                                                                     <  1:1147074/50‑1 (MQ=255)
                   tATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCg                                               <  1:3799434/71‑1 (MQ=255)
                     tGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCt                                             <  1:886354/71‑1 (MQ=255)
                       ccTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTc                                           <  1:3586117/71‑1 (MQ=255)
                        cTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg                                          >  1:2089950/1‑71 (MQ=255)
                          ttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTcc                                                >  1:2658201/1‑63 (MQ=255)
                          ttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg                                          >  1:3727501/1‑69 (MQ=255)
                              tGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa                                    <  1:158468/71‑1 (MQ=255)
                               gCAGATTTCACCGATCTCCCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa                                    >  1:2519795/1‑70 (MQ=255)
                                cAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaaaa                                  <  1:1856596/71‑1 (MQ=255)
                                cAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaaa                                   <  1:1205126/70‑1 (MQ=255)
                                   aTTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAc                                 >  1:2196250/1‑69 (MQ=255)
                                     ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCg                                               >  1:2720048/1‑53 (MQ=255)
                                     ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa                                    <  1:2254911/64‑1 (MQ=255)
                                       cACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCAcc                           >  1:2761246/1‑71 (MQ=255)
                                          cGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGa                        >  1:3130442/1‑71 (MQ=255)
                                          cGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGa                        >  1:3396342/1‑71 (MQ=255)
                                          cGATCTCACCAGCCGGGAGCGCGCGATTGTGATCGTCGCGGATCTCCGCt                                             >  1:3198314/1‑50 (MQ=255)
                                              ctcACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGAtct                                                  <  1:3164973/41‑1 (MQ=255)
                                              ctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCAGCTTCGTAGCAAAACCCCAccc                          <  1:750135/65‑1 (MQ=255)
                                                cACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATc                  >  1:212119/1‑71 (MQ=255)
                                                 aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCaa                    >  1:416049/1‑68 (MQ=255)
                                                 aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCa                     <  1:2889345/67‑1 (MQ=255)
                                                  ccAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGa                <  1:3683040/71‑1 (MQ=255)
                                                   cAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc               <  1:3800268/71‑1 (MQ=255)
                                                   cAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc               >  1:3153023/1‑71 (MQ=255)
                                                   cAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGAc               <  1:2025730/71‑1 (MQ=255)
                                                       cGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGc            >  1:2355506/1‑70 (MQ=255)
                                                                gCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGtt  >  1:937090/1‑71 (MQ=255)
                                                                 |                                                                     
TGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT  >  minE/34362‑34496

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: