Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,676,470 T→A 83.4% P73P (CCT→CCA bioH → carboxylesterase of pimeloyl‑CoA synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,676,4700TA83.4% 83.1 / 9.4 42P73P (CCT→CCAbioHcarboxylesterase of pimeloyl‑CoA synthesis
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base A (15/20);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

GGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTT  >  minE/2676413‑2676529
                                                         |                                                           
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:2797456/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:944463/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:756307/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:646395/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:494817/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:380540/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:365773/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:218651/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:2038306/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:12236/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:1513076/1‑68 (MQ=255)
gggATTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCa                                                   >  1:206680/1‑68 (MQ=255)
   aTTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCATTTg                                               >  1:2092178/1‑69 (MQ=255)
   aTTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCATTTg                                               >  1:2949905/1‑69 (MQ=255)
   aTTTGGTGCGCTATCACTTGCTGATATGGCCGAAGTCGTGCTGCGACAGGCACCAGATAAAGCCATTTg                                               >  1:248100/1‑69 (MQ=255)
          gcgcTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGc                                        <  1:3217080/69‑1 (MQ=255)
          gcgcTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGc                                        <  1:808959/69‑1 (MQ=255)
          gcgcTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGc                                        <  1:2145700/69‑1 (MQ=255)
          gcgcTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGc                                        <  1:2362119/69‑1 (MQ=255)
                               gCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCtggt                    >  1:1237504/1‑68 (MQ=255)
                               gCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCtggt                    >  1:627401/1‑68 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGTAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:748157/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2200427/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:1162395/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:860025/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:840830/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2193695/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:700041/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:559603/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2432183/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:183190/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:322476/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2406677/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2879801/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:1045963/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2640189/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2614997/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:2600790/69‑1 (MQ=255)
                                                aCAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCCGGCTGGTGGCAAGCCAGATTGCGtt  <  1:3068496/69‑1 (MQ=255)
                                                 cAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCATATTGCGtt  <  1:1721778/68‑1 (MQ=255)
                                                 cAGGCACCAGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGtt  <  1:836879/68‑1 (MQ=255)
                                                  aGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAg          <  1:669730/59‑1 (MQ=255)
                                                         |                                                           
GGGATTTGGTGCGCTGTCACTTGCTGATATGGCCGAAGCCGTGCTGCAACAGGCACCTGATAAAGCCATTTGGTTAGGCTGGAGTCTGGGCGGGCTGGTGGCAAGCCAGATTGCGTT  >  minE/2676413‑2676529

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: