Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,408,526 A→G 46.9% P265P (CCA→CCG 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,5260AG46.9% 2.0 / 54.6 64P265P (CCA→CCGyejFfused predicted oligopeptide transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (8/26);  new base G (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 = 7.97e-01

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