Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,532 G→A 46.8% P267P (CCG→CCA yejF → fused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,408,5320GA46.8% 6.4 / 52.2 64P267P (CCG→CCAyejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (8/26);  new base A (14/16);  total (22/42)
Fisher's exact test for biased strand distribution p-value = 6.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CACCCACTCATCCTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCATT  >  minE/1408471‑1408592
                                                             |                                                            
cACCCACTCATCCTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAg                                                       <  1:3286970/69‑1 (MQ=255)
   ccACTCATCCTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGaa                                                     <  1:805914/68‑1 (MQ=255)
   ccACTCATCCTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGaa                                                     <  1:749067/68‑1 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:686253/1‑69 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:2867504/1‑69 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:682365/1‑69 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:170363/1‑69 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:3338646/1‑69 (MQ=255)
            cTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCaa                                           >  1:1769735/1‑69 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTTAACGTTGCTGGATGTTg                            <  1:3023576/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:2791023/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:2933077/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:3130351/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:3207870/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:2110158/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:1893609/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:1739016/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:1720089/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:519928/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:152035/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:1308846/69‑1 (MQ=255)
                           tACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTg                            <  1:1249627/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:3405314/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:335550/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:395955/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:3350410/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:3348959/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:3293252/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:460882/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:3140664/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:918743/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:2731508/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:1378990/69‑1 (MQ=255)
                            aCTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:2491926/69‑1 (MQ=255)
                             cTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:2800901/68‑1 (MQ=255)
                             cTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:2955465/68‑1 (MQ=255)
                             cTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:494897/68‑1 (MQ=255)
                             cTCAACAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGa                           <  1:1872699/68‑1 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATTTTGAACAGc                      >  1:2548203/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:2759635/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:833307/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:1088672/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:535454/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:446865/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:2226866/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:2906270/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:3148283/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:2305741/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:2577156/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:3126153/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCTAGATGTTGAACAGc                      >  1:3034219/1‑68 (MQ=255)
                                  cAGTGAACCATCTGGCGACCCGGTGCCATTGCCAGAACCCGCCTCAACGTTGCCAGATGTTGAACAGc                      >  1:1125806/1‑68 (MQ=255)
                                   aGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGc                      >  1:1671635/1‑67 (MQ=255)
                                      gAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAg                 >  1:3171000/1‑69 (MQ=255)
                                                     ccAGTGCCGTTGGCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:2125586/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:1729434/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:1002156/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:2630966/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:1208985/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:1915034/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:1992007/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:2372658/69‑1 (MQ=255)
                                                     ccAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:2290675/69‑1 (MQ=255)
                                                     ccAGTGCCGGTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCAtt  <  1:443479/69‑1 (MQ=255)
                                                             |                                                            
CACCCACTCATCCTTACACACAAAAGCTACTCAACAGTGAACCGTCAGGCGATCCAGTGCCGTTGCCAGAACCTGCCTCAACGTTGCTGGATGTTGAACAGCTTCAGGTTGCCTTCCCCATT  >  minE/1408471‑1408592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: