Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,193,222 T→A 47.2% V147V (GTA→GTT fadD ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,193,2220TA47.2% 1.9 / 82.2 70V147V (GTA→GTTfadDacyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (21/12);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TTTAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACA  >  minE/1193156‑1193288
                                                                  |                                                                  
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTGAGAATTAc                                                                  >  1:1930767/1‑69 (MQ=255)
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:381919/1‑69 (MQ=255)
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:857753/1‑69 (MQ=255)
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:1294466/1‑69 (MQ=255)
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:581883/1‑69 (MQ=255)
tttAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTAc                                                                  >  1:2067728/1‑69 (MQ=255)
      aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:1013237/69‑1 (MQ=255)
      aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:2703162/69‑1 (MQ=255)
      aacGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTg                                                            <  1:678085/69‑1 (MQ=255)
                  tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:61011/69‑1 (MQ=255)
                  tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:2068553/69‑1 (MQ=255)
                  tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:199135/69‑1 (MQ=255)
                  tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:1528251/69‑1 (MQ=255)
                  tACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGcttt                                                <  1:1399077/69‑1 (MQ=255)
                   aCCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGttt                                                 <  1:1060670/67‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2479076/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3387999/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3200109/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2904346/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:506101/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1123599/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2526621/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:826259/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1755477/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1873045/68‑1 (MQ=255)
                                  gTAGATAGCTGATCGCCCATACGGGTCAGAATAACGTGCTGAACGGCGGGTTTATCAACCACTTTTTc                                 <  1:63189/68‑1 (MQ=255)
                                     gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTTTCAACCACTTTTTCCAgt                             >  1:131118/1‑69 (MQ=255)
                                     gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:1465220/1‑69 (MQ=255)
                                     gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:2117789/1‑69 (MQ=255)
                                     gATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgt                             >  1:1846808/1‑69 (MQ=255)
                                        aGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgtgt                           >  1:1376675/1‑68 (MQ=255)
                                        aGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAgtgt                           >  1:388671/1‑68 (MQ=255)
                                                  ccATACGGGTCAGAATTACGTGCTGAACGGCGGtttt                                                <  1:305981/37‑1 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3034641/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2990387/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:1620330/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3102201/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:805130/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:3280145/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2896607/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:285904/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2778637/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2716363/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2712581/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:266609/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2424192/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:508315/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2413383/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:2285532/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:614062/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:1325561/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:85307/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAAACACTTTTTCCAGAGTGTGGGCAAAGTTa              >  1:409850/1‑68 (MQ=255)
                                                     tACGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAAACACTTTTTCCAGAGTGCGGGCAAAGTTa              >  1:2500694/1‑68 (MQ=255)
                                                      aCGGGTCAGAATAACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:992221/48‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3164954/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:1005772/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:988995/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:683618/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:602527/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:2322562/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:450390/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTc                                 <  1:3306108/47‑1 (MQ=255)
                                                       cGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGac           >  1:2180596/1‑69 (MQ=255)
                                                            cAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:662781/1‑69 (MQ=255)
                                                            cAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2976133/1‑69 (MQ=255)
                                                             aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCTAAGTTAGACACGAt      >  1:935892/1‑68 (MQ=255)
                                                             aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2423991/1‑68 (MQ=255)
                                                             aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:2130187/1‑68 (MQ=255)
                                                             aGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGAt      >  1:1833096/1‑68 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:2711638/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1839383/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1736746/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1629760/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1358054/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:1185488/68‑1 (MQ=255)
                                                                 ttACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACa  <  1:3135452/68‑1 (MQ=255)
                                                                  |                                                                  
TTTAACAACGAAATTGACTACCGTGCCTTTTGCCGTAGATAGCTGATCGCCCATACGGGTCAGAATTACGTGCTGAACGGCGGTTTTATCAACCACTTTTTCCAGTGTGTGAGCAAAGTTAGACACGATAACA  >  minE/1193156‑1193288

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: