Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,368,808 A→G 63.9% G389G (GGA→GGG 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,8080AG63.9% 36.3 / 32.3 58G389G (GGA→GGGfadBfused 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 A (12/9);  new base G (9/28);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 2.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

AAGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGAAAAAAGC  >  minE/2368743‑2368867
                                                                 |                                                           
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:2547356/1‑69 (MQ=255)
 aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                                                         >  1:2609505/1‑69 (MQ=255)
 aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                                                         >  1:1163726/1‑69 (MQ=255)
 aGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTg                                                         >  1:2952528/1‑69 (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:234460/69‑1 (MQ=255)
  gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                                                        <  1:1874574/69‑1 (MQ=255)
  gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                                                        <  1:468237/69‑1 (MQ=255)
  gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                                                        <  1:1466281/69‑1 (MQ=255)
  gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                                                        <  1:3000325/69‑1 (MQ=255)
  gATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGa                                                        <  1:2840122/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)
        tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                                                   >  1:523760/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:2453271/1‑68 (MQ=255)
        tGGTCTGAAACTGGCTGGCGTGATCTCCACAATTCACCCAACGCTCGACTACGCCGGGTTTGACcgcg                                                   >  1:825187/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)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:2608145/69‑1 (MQ=255)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:365776/69‑1 (MQ=255)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:3378632/69‑1 (MQ=255)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:1110861/69‑1 (MQ=255)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:305797/69‑1 (MQ=255)
                                            cccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:2994369/69‑1 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATtgt                                         >  1:1097065/1‑41 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaaa             >  1:2444126/1‑69 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaaa             >  1:3269832/1‑69 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaaa             >  1:2838055/1‑69 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaaa             >  1:2420718/1‑69 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaa              <  1:1093109/68‑1 (MQ=255)
                                             ccAACGCTCGACTACGCCGGATTTGACAGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGaaa             >  1:100157/1‑69 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTTTTGAAAACCCGAAAGTGaaaa      <  1:2238829/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAATTGaaaa      <  1:3649/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:1027941/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:785278/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:2004492/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:2108912/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:2239766/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:3191383/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:245142/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:2514123/69‑1 (MQ=255)
                                                    tCGACTACGCCGGGTTTGACCGCGTGGATGTTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGaaaa      <  1:2536432/69‑1 (MQ=255)
                                                        cTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGAAAAAAGc  <  1:39911/69‑1 (MQ=255)
                                                                 |                                                           
AAGATCGATGGTCTGAAACTGGCTGGCGTGATCTCCACAATCCACCCAACGCTCGACTACGCCGGATTTGACCGCGTGGATATTGTGGTAGAAGCGGTTGTTGAAAACCCGAAAGTGAAAAAAGC  >  minE/2368743‑2368867

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: