Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,694,489 G→A 25.5% P47S (CCT→TCT)  yrfF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,694,4890GA25.5% 89.7 / 24.8 74P47S (CCT→TCT) yrfFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (42/13);  new base A (9/10);  total (51/23)
Fisher's exact test for biased strand distribution p-value = 2.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GCCAGGCACCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAT  >  minE/2694424‑2694553
                                                                 |                                                                
gCCAGGCACCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGgtg                                                                <  1:2051567/68‑1 (MQ=255)
gCCAGGCACCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGgtg                                                                <  1:1600875/68‑1 (MQ=255)
       aCCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTa                                                        >  1:630736/1‑69 (MQ=255)
       aCCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTa                                                        >  1:2728564/1‑69 (MQ=255)
       aCCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTa                                                        >  1:1035381/1‑69 (MQ=255)
        ccTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAg                                                            >  1:82578/1‑64 (MQ=255)
        ccTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAg                                                            >  1:311526/1‑64 (MQ=255)
        ccTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAg                                                            >  1:2393884/1‑64 (MQ=255)
        ccTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAg                                                            >  1:1838279/1‑64 (MQ=255)
        ccTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAg                                                            >  1:973902/1‑64 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:868947/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:3180384/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:2730815/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:920551/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:971383/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGtgtg                                                              >  1:2196053/1‑56 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGttt                                                 >  1:3059977/1‑69 (MQ=255)
              aTACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGttt                                                 >  1:2500944/1‑69 (MQ=255)
                aCCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGtt                                                          >  1:1534630/1‑58 (MQ=255)
                aCCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGtt                                                          >  1:1219496/1‑58 (MQ=255)
                 ccTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGtt                                                          <  1:456701/57‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGGGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:1303165/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:3020480/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:324383/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:3052000/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:3039333/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:1928247/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:1927046/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:1635883/68‑1 (MQ=255)
                       tCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:679576/68‑1 (MQ=255)
                           gCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGa                                     >  1:1816737/1‑68 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:609199/65‑1 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:2825919/65‑1 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:2620725/65‑1 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:1488918/65‑1 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:779762/65‑1 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAATGCGtt                              >  1:633111/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:2756827/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:408485/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:982593/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:974277/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:3164653/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:488625/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:3435949/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:1504206/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:1477783/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:1721485/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:1225059/1‑69 (MQ=255)
                                 aaGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:1772487/1‑69 (MQ=255)
                                 aaGATAAATTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGtt                              >  1:2601403/1‑69 (MQ=255)
                                  aGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgg                                  <  1:3225217/64‑1 (MQ=255)
                                     tAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTgg                            >  1:2349629/1‑67 (MQ=255)
                                     tAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTgg                            >  1:924374/1‑67 (MQ=255)
                                         tttCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAg                       <  1:2847965/68‑1 (MQ=255)
                                         tttCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAg                       <  1:783875/68‑1 (MQ=255)
                                         tttCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAg                       <  1:909892/68‑1 (MQ=255)
                                         tttCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAg                       <  1:1262930/68‑1 (MQ=255)
                                         tttCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAg                       <  1:1162279/68‑1 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGc                 >  1:218484/1‑68 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:1345787/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:1205589/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:3028005/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:1102189/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:617175/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:2044181/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:25159/1‑69 (MQ=255)
                                               cGGCGCTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATCCACGAAGGCGTTGGTCCAGGGCAGCt                >  1:2781338/1‑69 (MQ=255)
                                                    cTACGTTCTTCAGATGTGAGTTTACGCGTTTGCGCATcc                                         <  1:803264/39‑1 (MQ=255)
                                                      aCGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGccgccg         >  1:941767/1‑69 (MQ=255)
                                                       cGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAgt                                  <  1:2813962/43‑2 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAt  >  1:1482280/1‑68 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAt  >  1:1563939/1‑68 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAt  >  1:1759979/1‑68 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAt  >  1:189358/1‑68 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAt  >  1:2203197/1‑68 (MQ=255)
                                                              cAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTATGGAt  >  1:822611/1‑68 (MQ=255)
                                                                 |                                                                
GCCAGGCACCTGTAATACCTGCGTCAGGCTCTCAAGATAATTTTCAACGGCGCTACGTTCTTCAGGTGTGAGTTTACGCGTTTGCGCATCCGCGAAGGCGTTGGTCCAGGGCAGCTGACGCCGTCTGGAT  >  minE/2694424‑2694553

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