Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,368,784 C→T 77.3% I381I (ATC→ATT fadB → fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,368,7840CT77.3% 42.2 / 7.4 35I381I (ATC→ATTfadBfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base T (10/17);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGAACAAGCAGCTTGAGCGCGGCAAGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTT  >  minE/2368718‑2368841
                                                                  |                                                         
gctgAACAAGCAGCTTGAGCGCGGCAAGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCa                                                         <  1:1104910/69‑1 (MQ=255)
                         aaGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCAGAATTCACCCAACGCTCGACTACGCCGGGttt                                <  1:2107447/69‑1 (MQ=255)
                         aaGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGttt                                <  1:1310956/69‑1 (MQ=255)
                         aaGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGttt                                <  1:1459754/69‑1 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGttt                                <  1:1040851/68‑1 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGAttt                                >  1:503975/1‑68 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                               >  1:1163726/1‑69 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                               >  1:2952528/1‑69 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                               >  1:2609505/1‑69 (MQ=255)
                          aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                               >  1:2547356/1‑69 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:468237/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:3000325/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:2840122/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:234460/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:2168927/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:511113/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:579788/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:1874574/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:1466281/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:934923/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:977662/69‑1 (MQ=255)
                           gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCGCGACTACGCCGGGTTTGa                              <  1:994982/69‑1 (MQ=255)
                            aTCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                              <  1:2859864/68‑1 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCa                                                     >  1:2512623/1‑40 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:523760/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:825187/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:2453271/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:2108030/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:1972981/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:1760199/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:1355434/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:1327877/1‑68 (MQ=255)
                                 tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                         >  1:1083495/1‑68 (MQ=255)
                                      tGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGata                   <  1:1880568/69‑1 (MQ=255)
                                                       aTCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGgtt  >  1:2482662/1‑69 (MQ=255)
                                                                  |                                                         
GCTGAACAAGCAGCTTGAGCGCGGCAAGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTT  >  minE/2368718‑2368841

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: