Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,196,732 C→A 32.4% G236G (GGC→GGA degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,196,7320CA32.4% 88.0 / 46.4 74G236G (GGC→GGAdegQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base A (13/11);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTGG  >  minE/2196674‑2196787
                                                          |                                                       
cACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGc                                                >  1:2567509/1‑68 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:2386400/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:250226/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:781797/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:548455/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:2064589/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:2821365/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:3202322/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:2022614/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:2873612/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:1620094/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:3091647/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:1711737/1‑69 (MQ=255)
     ttAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCggg                                          >  1:1744775/1‑69 (MQ=255)
          ccTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGggattgg                                     >  1:2288996/1‑69 (MQ=255)
             tAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGggattggatt                                  >  1:2908533/1‑69 (MQ=255)
             tAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGggattggat                                   >  1:1371695/1‑68 (MQ=255)
             tAACGGTGAGTTAATTGGCAGCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGggattggatt                                  >  1:3052392/1‑69 (MQ=255)
                cGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGAttt                                 <  1:501677/67‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGAttt                                 <  1:1353307/66‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:3128764/69‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:1395610/69‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:1358447/69‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:1236261/69‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:2353887/69‑1 (MQ=255)
                 ggTGAGTTAATTGGCATCAACACTGCAATCCTTGAGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGcc                              <  1:238179/69‑1 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATc                           <  1:857958/69‑1 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATc                           <  1:2582798/69‑1 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2233128/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2173058/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:1863498/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2448924/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2561375/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:1361101/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:332380/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2921845/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:2791998/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:202096/69‑2 (MQ=255)
                    gAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGAGGCGGTAGCGTCGGGATTGGATTTGCCAtt                           <  1:3297913/69‑2 (MQ=255)
                            tGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                   <  1:3195417/69‑1 (MQ=255)
                            tGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                   <  1:653968/69‑1 (MQ=255)
                            tGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                   <  1:2588073/69‑1 (MQ=255)
                            tGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                   <  1:203164/69‑1 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGc                               >  1:1711120/1‑50 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGc                               >  1:289846/1‑50 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGc                               >  1:2770656/1‑50 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGc                               >  1:802172/1‑50 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGc                               >  1:80237/1‑50 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATg               >  1:3416800/1‑66 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATg               >  1:3186949/1‑66 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:2207167/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:2955960/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:1259986/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:3359550/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:2140669/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:376986/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:384932/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:459812/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:921555/1‑68 (MQ=255)
                                   aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgc             >  1:1329867/1‑68 (MQ=255)
                                       cTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGaa        <  1:1094063/69‑1 (MQ=255)
                                            aTCCTTTCGCCTGGCGGCGGGAGCGTCGGGATTGGAttt                                 <  1:2152894/39‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGTGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg          <  1:180876/59‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGCGCGTCGGGATTGGAttt                                 >  1:234864/1‑36 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTg                                >  1:2430140/1‑37 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg          <  1:744283/59‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg          <  1:3296262/59‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg          <  1:818483/59‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg          <  1:854298/59‑1 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGGACACTgg  >  1:3120671/1‑67 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTgg  >  1:9778/1‑67 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTgg  >  1:1282674/1‑67 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTgg  >  1:1347769/1‑67 (MQ=255)
                                               cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTgg  >  1:1870436/1‑67 (MQ=255)
                                                          |                                                       
CACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTGG  >  minE/2196674‑2196787

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: