Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 281,478 T→C 25.8% L179L (TTA→TTG ampG ← muropeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE281,4780TC25.8% 58.9 / 18.5 54L179L (TTA→TTGampGmuropeptide transporter
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (12/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGGCCGCC  >  minE/281410‑281537
                                                                    |                                                           
ggAATGGTGTCGGTTGGTTCTGGTGCGAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2327439/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACGGTGCCGCCAtt                                                             >  1:3076646/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:3193327/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2966599/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2881519/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:478768/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:3311158/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2600785/1‑69 (MQ=255)
ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2702835/1‑69 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2052063/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2024702/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1958790/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2187735/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2551033/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2769641/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:881674/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1356333/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1224604/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:3041941/1‑67 (MQ=255)
   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:3106592/1‑67 (MQ=255)
      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:3254765/69‑1 (MQ=255)
      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:2286368/69‑1 (MQ=255)
      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:1792466/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1875300/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:882946/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:115679/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:3050349/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1506766/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:412492/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1687028/69‑1 (MQ=255)
        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:186186/69‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:181923/61‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:582603/61‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:3409504/61‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:1098489/61‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:1266137/61‑1 (MQ=255)
           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACAt                                                  >  1:2261835/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:249397/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:814367/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1084346/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:3137295/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:3105012/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1170801/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1415965/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:335578/1‑69 (MQ=255)
                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1515012/1‑69 (MQ=255)
                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:753285/1‑68 (MQ=255)
                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:2824454/1‑68 (MQ=255)
                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:1654413/1‑68 (MQ=255)
                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:1045652/1‑69 (MQ=255)
                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:361394/1‑69 (MQ=255)
                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:2491317/1‑69 (MQ=255)
                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:2108664/1‑69 (MQ=255)
                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:3219127/69‑1 (MQ=255)
                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:1488934/69‑1 (MQ=255)
                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:1013808/69‑1 (MQ=255)
                                          ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   <  1:2824793/69‑1 (MQ=255)
                                          ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   >  1:2615649/1‑69 (MQ=255)
                                            aCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   >  1:2838540/1‑67 (MQ=255)
                                             cAGGGGATCAACAGTGTCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCcac                <  1:2651231/69‑1 (MQ=255)
                                                     cAACAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa        <  1:2616209/69‑2 (MQ=255)
                                                      aaCAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa        <  1:2592896/68‑2 (MQ=255)
                                                            gccgccATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGgccgcc  >  1:229712/1‑68 (MQ=255)
                                                                    |                                                           
GGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGGCCGCC  >  minE/281410‑281537

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: