Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 34,045 T→C 42.8% Q491Q (CAA→CAG 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,0450TC42.8% 15.6 / 47.6 63Q491Q (CAA→CAGcaiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (15/12);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCT  >  minE/33987‑34111
                                                          |                                                                  
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAg                                                           >  1:94082/1‑68 (MQ=255)
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGa                                                          >  1:2855035/1‑69 (MQ=255)
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGa                                                          >  1:2273947/1‑69 (MQ=255)
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGa                                                          >  1:1341566/1‑69 (MQ=255)
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGa                                                          >  1:1401290/1‑69 (MQ=255)
tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGa                                                          >  1:1643358/1‑69 (MQ=255)
 aCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                         >  1:1122937/1‑69 (MQ=255)
 aCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                         >  1:3025273/1‑69 (MQ=255)
 aCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                         >  1:2127830/1‑69 (MQ=255)
 aCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                         >  1:2102683/1‑69 (MQ=255)
 aCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                         >  1:3121161/1‑69 (MQ=255)
            cTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAAtt                                              >  1:1035703/1‑69 (MQ=255)
            cTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAAtt                                              >  1:2199895/1‑69 (MQ=255)
                 cTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2549786/65‑1 (MQ=255)
                 cTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2563313/65‑1 (MQ=255)
                  tGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                             <  1:2579541/64‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:3086953/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2537328/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:302927/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:213673/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:130284/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2016752/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2586/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:1664055/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2656317/65‑1 (MQ=255)
                           gATAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGAAGAAGAATTCCTCTTCGCTc                                   <  1:1963413/65‑1 (MQ=255)
                            aTAAGAGGGCACTTTAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:251102/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2086864/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:1273347/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2281558/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2929916/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:1793624/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2653237/64‑1 (MQ=255)
                              aaGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                 <  1:2717591/64‑1 (MQ=255)
                                 agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:2828453/1‑69 (MQ=255)
                                 agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:826738/1‑69 (MQ=255)
                                 agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:851926/1‑69 (MQ=255)
                                 agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                         >  1:847984/1‑69 (MQ=255)
                                    gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:96204/69‑1 (MQ=255)
                                    gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:218100/69‑1 (MQ=255)
                                    gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:1345692/69‑1 (MQ=255)
                                    gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:2631179/69‑1 (MQ=255)
                                     cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      <  1:262235/68‑1 (MQ=255)
                                        tttAAATTTCGCCATATTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTc                                   <  1:2591434/52‑1 (MQ=255)
                                                   ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:2166235/68‑1 (MQ=255)
                                                   ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:2684036/68‑1 (MQ=255)
                                                   ccATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACcaca        <  1:1596978/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:2960105/1‑69 (MQ=255)
                                                       atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1036620/1‑69 (MQ=255)
                                                       atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:2766186/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:145093/1‑69 (MQ=255)
                                                       atTCTGTTCACAGAAGCAGAAGAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGc   >  1:1541815/1‑69 (MQ=255)
                                                        ttttttTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCt  >  1:1445782/1‑69 (MQ=255)
                                                        ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:2663208/1‑49 (MQ=255)
                                                        ttttGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCtt                      >  1:990473/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)
                                                          |                                                                  
TACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCT  >  minE/33987‑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: