Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 923,755 G→A 34.5% G197S (GGT→AGT)  mppA → murein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE923,7550GA34.5% 59.7 / 44.1 81G197S (GGT→AGT) mppAmurein tripeptide (L‑ala‑gamma‑D‑glutamyl‑meso‑DAP) transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (16/37);  new base A (20/8);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 8.23e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

GCCGTGGTTTGTGAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTA  >  minE/923688‑923809
                                                                   |                                                      
gCCGTGGTTTGTGAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCg                                                        >  1:169805/1‑68 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:2975304/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:1275425/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:414254/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:2395983/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:724243/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:2757523/1‑69 (MQ=255)
         tgtgAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTg                                              >  1:1631709/1‑69 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:2202344/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:2691679/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:3417280/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:3159085/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:1804708/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:2413331/65‑1 (MQ=255)
             aaTTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:336762/65‑1 (MQ=255)
             aaTTTAACCCCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTg                                              <  1:3119369/65‑1 (MQ=255)
               tttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaa                                        >  1:2713360/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTgg                                             >  1:2070566/1‑63 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:2190668/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:1570688/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:249731/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:1815165/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:947265/1‑69 (MQ=255)
                ttAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGaaa                                       >  1:1292667/1‑69 (MQ=255)
                            tttGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTg                           <  1:1925555/69‑1 (MQ=255)
                              tGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGAt                         >  1:2615825/1‑69 (MQ=255)
                                ccTTCTTCCCGGTGCAAAAAGGCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATc                        <  1:1861062/68‑1 (MQ=255)
                                ccTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATc                        <  1:1580740/68‑1 (MQ=255)
                                ccTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATc                        <  1:1512129/68‑1 (MQ=255)
                                  ttcttcCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:2340820/69‑1 (MQ=255)
                                  ttcttcCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:2486402/69‑1 (MQ=255)
                                  ttcttcCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:157138/69‑1 (MQ=255)
                                  ttcttcCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:3094948/69‑1 (MQ=255)
                                  ttcttcCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:3055490/69‑1 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGAc                                           >  1:2623293/1‑40 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:3134903/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2745618/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2805665/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:3032285/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2166842/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:3310546/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2359369/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2296435/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:2241590/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:364083/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:408854/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:45974/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1681731/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1630464/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:868058/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1422249/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1075236/1‑53 (MQ=255)
                                         cGGTGCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1007577/1‑53 (MQ=255)
                                         cGGTCCAAAAAGCCAACGTAGAAAGCAGTAAAGAGTGGACGAAACCCGGAAAt                              >  1:1557602/1‑53 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:3374037/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:99436/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:3005833/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:1520255/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:1555740/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:2118961/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:2399991/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:2723539/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:2767507/69‑1 (MQ=255)
                                           gTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCt            <  1:2805731/69‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:2447974/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:2289614/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:3420883/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:3422474/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:349924/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:1938417/57‑1 (MQ=255)
                                              cAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:872126/57‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:948534/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:1124009/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:919997/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:1015291/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:1883956/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:2146159/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:2512600/69‑1 (MQ=255)
                                                     ccAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:2810477/69‑1 (MQ=255)
                                                      cAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTa  <  1:1983289/68‑1 (MQ=255)
                                                          gTAGAAAGCGGTAAAGAGTGGACGAAACCCGTAAATCTGATCGGc                     <  1:918443/45‑1 (MQ=255)
                                                              aaaGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGc                     <  1:825151/41‑1 (MQ=255)
                                                                   |                                                      
GCCGTGGTTTGTGAATTTAACCGCTAACTTTGCCTTCTTCCCGGTGCAAAAAGCCAACGTAGAAAGCGGTAAAGAGTGGACGAAACCCGGAAATCTGATCGGCAATGGCGCTTATGTTCTTA  >  minE/923688‑923809

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: