Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,217,619 A→G 54.3% V332V (GTT→GTC yhdA ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,217,6190AG54.3% 8.2 / 46.8 59V332V (GTT→GTCyhdAconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (7/20);  new base G (14/18);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAACAGCGCAG  >  minE/2217551‑2217686
                                                                    |                                                                   
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:2197534/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:304308/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:1110916/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:1230950/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:2955774/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:2732849/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:2201457/1‑69 (MQ=255)
cGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCa                                                                     >  1:2322249/1‑69 (MQ=255)
       cgcAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCt                                                               >  1:333303/1‑68 (MQ=255)
        gaaGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCttt                                                             <  1:2420649/67‑1 (MQ=255)
            aTACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGAtt                                                         <  1:980103/69‑1 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:2465533/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:654389/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:613287/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:2618708/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:69644/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:503314/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:3258640/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:699252/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:94570/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:846156/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:1256353/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:936984/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:900189/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTaa                                                       >  1:1410354/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTaa                                                       >  1:3331023/1‑69 (MQ=255)
              aCCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTa                                                        >  1:349857/1‑68 (MQ=255)
                 aaTGTGTATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:1023593/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGATCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:690668/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:712431/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:124343/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:428778/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2630656/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:922758/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:1542725/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:250481/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:1587354/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2423475/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:1614270/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2420251/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2236978/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:716448/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGAtt                                                         <  1:7380/64‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACt                                                     <  1:2769940/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACt                                                     <  1:1043311/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACt                                                     <  1:1665448/68‑1 (MQ=255)
                 aaTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTAGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2022251/68‑1 (MQ=255)
                  aTGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:2440388/67‑1 (MQ=255)
                       gATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGc                                              <  1:902589/69‑1 (MQ=255)
                       gATCATATCGTCGCGATAGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGc                                              <  1:2577109/69‑1 (MQ=255)
                        aTCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCGACGGCTTTGATTAACt                                                     <  1:938957/61‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:2629802/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:728442/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:1784126/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:1972442/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:2269546/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:1635988/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:1361947/69‑1 (MQ=255)
                             atCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:2856701/69‑1 (MQ=255)
                              tcgtcgCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGct                                        <  1:366788/68‑1 (MQ=255)
                                  cgcgATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTg                                   >  1:1497102/1‑69 (MQ=255)
                                     gATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTa                                                        <  1:489997/45‑1 (MQ=255)
                                     gATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTa                                                        <  1:3246064/45‑1 (MQ=255)
                                     gATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTa                                                        <  1:306763/45‑1 (MQ=255)
                                                ttATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTc                     >  1:2950732/1‑69 (MQ=255)
                                                ttATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTc                     >  1:3435812/1‑69 (MQ=255)
                                                ttATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTc                     >  1:490641/1‑69 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAAc         >  1:1564192/1‑61 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAAc         >  1:478418/1‑61 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAAc         >  1:467218/1‑61 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAAc         >  1:2701301/1‑61 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAACAGCgcag  >  1:2455602/1‑68 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAACAGCgcag  >  1:1055134/1‑68 (MQ=255)
                                                                    aaCGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAACAGCgcag  >  1:3196702/1‑68 (MQ=255)
                                                                    |                                                                   
CGCCAGGCGCAGATACCAATGTGGATCATATCGTCGCGATCGAGCATTTTATTGTTCGGCAAGGTATCAACGGCTTTGATTAACTGACCGGCGATGCTCTCTGCCTCTTTTAACGTCCGGTGCGGTAACAGCGCAG  >  minE/2217551‑2217686

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: