Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,721,260 C→T 51.5% E107E (GAG→GAA yhfK ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,721,2600CT51.5% ‑4.4 / 59.0 66E107E (GAG→GAAyhfKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base T (9/25);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GAGGGTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAG  >  minE/2721200‑2721328
                                                            |                                                                    
gAGGGTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGt                                                              <  1:757908/69‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:3195695/68‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:1777712/68‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:2807031/68‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:370637/68‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:2156568/68‑1 (MQ=255)
   ggTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGa                                                            <  1:3052090/68‑1 (MQ=255)
       aaaaTGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCTGTGACGCCa                                                       <  1:1091488/69‑1 (MQ=255)
       aaaaTGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCa                                                       <  1:2498572/69‑1 (MQ=255)
       aaaaTGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCa                                                       <  1:2589290/69‑1 (MQ=255)
       aaaaTGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCa                                                       <  1:2636010/69‑1 (MQ=255)
                         gATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCg                                     <  1:406894/69‑1 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATc                                       >  1:3415539/1‑58 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGc                               >  1:530197/1‑66 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaa                            <  1:3047387/69‑1 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:420581/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:24396/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:2594665/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:2714398/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:2849404/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:928835/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:1301223/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:1873526/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:1489261/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:229319/1‑67 (MQ=255)
                                  aGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCa                              >  1:3145444/1‑67 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCTTTAATCCGGTCAGCaaaaa                          <  1:1485850/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGTTCAGCaaaaa                          <  1:505749/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:355829/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2522331/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2510253/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2203342/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:752723/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:950795/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1275453/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1278758/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1384001/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:924278/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1721888/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1724759/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:1865487/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2403511/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:746332/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2135764/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2150465/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:594070/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2213610/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2263909/69‑1 (MQ=255)
                                    cAATTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2352229/69‑1 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGTTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCaaaaa                          <  1:2458052/68‑1 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCg                                     <  1:1312164/57‑1 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:2877479/1‑69 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:2545001/1‑69 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:2311847/1‑69 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:2194076/1‑69 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:2155670/1‑69 (MQ=255)
                                     aaTTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAg                         >  1:187678/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTGCCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:929698/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:2159442/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:1083356/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:1543676/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:179223/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:315212/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:1019926/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:2572990/1‑69 (MQ=255)
                                                     gCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCtt         >  1:325127/1‑69 (MQ=255)
                                                            cTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCag  <  1:208910/69‑1 (MQ=255)
                                                            |                                                                    
GAGGGTAAAAATGGCGGCGAGCAGCGATGCAGGAAGCAATTTTGCGTGCAATGGCCCCAGCTCAGCAGTGACGCCAAGTACCAGCGTTAATCCGGTCAGCAAAAAGGGCAGGGGAACATCTTTTGCCAG  >  minE/2721200‑2721328

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: