Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 83,853 A→C 48.2% I237I (ATA→ATC ftsW → integral membrane protein involved in stabilising FstZ ring during cell division

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE83,8530AC48.2% ‑3.6 / 99.0 81I237I (ATA→ATCftsWintegral membrane protein involved in stabilising FstZ ring during cell division
Reads supporting (aligned to +/- strand):  ref base A (17/25);  new base C (18/21);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATCC  >  minE/83792‑83912
                                                             |                                                           
cGAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCg                                                      <  1:2237772/69‑1 (MQ=255)
 gAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGa                                                     <  1:701378/69‑1 (MQ=255)
 gAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGa                                                     <  1:1456784/69‑1 (MQ=255)
 gAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGa                                                     <  1:1236373/69‑1 (MQ=255)
 gAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGa                                                     <  1:769251/69‑1 (MQ=255)
  aaaTTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGa                                                     <  1:2965113/68‑1 (MQ=255)
                 aTTGCCATTATCGGTATTGGCATTTCAGCGGTTGTGTTGCTGATCCTCGc                                                        >  1:1953765/1‑50 (MQ=255)
                 aTTGCCATTATCGGTATGGGGATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:97300/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:3311185/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:256709/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:709608/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:1273771/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:2157824/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:3353944/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:2135541/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:2621414/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:1694388/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:1030667/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:2099973/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:17804/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:366192/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccgc                                     >  1:2831272/1‑69 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccg                                      >  1:1088743/1‑68 (MQ=255)
                 aTTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATCCTCGCCGAACCGTACCGTATccg                                      >  1:836036/1‑68 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:2650890/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:2511996/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:2298511/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:3231679/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1757532/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1739896/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:687742/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1450448/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1236789/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1098594/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCgtgt                                <  1:1038032/65‑1 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTAc                             >  1:1826778/1‑68 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTAc                             >  1:1825767/1‑68 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTAc                             >  1:2113486/1‑68 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTAc                             >  1:2187931/1‑68 (MQ=255)
                          aTCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTAc                             >  1:978764/1‑68 (MQ=255)
                             ggTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCa                         >  1:49663/1‑69 (MQ=255)
                             ggTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCa                         >  1:2265279/1‑69 (MQ=255)
                                aTGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTc                      >  1:2775796/1‑69 (MQ=255)
                                aTGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTc                      >  1:1520839/1‑69 (MQ=255)
                                aTGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTc                      >  1:3109122/1‑69 (MQ=255)
                                aTGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTc                      >  1:1561848/1‑69 (MQ=255)
                                aTGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTc                      >  1:3211888/1‑69 (MQ=255)
                                   ggCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAaccgtaccg                                            >  1:2572286/1‑44 (MQ=255)
                                   ggCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAaccgtaccg                                            >  1:1076254/1‑44 (MQ=255)
                                   ggCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAaccgtaccg                                            >  1:1817003/1‑44 (MQ=255)
                                   ggCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAaccgtaccg                                            >  1:1650013/1‑44 (MQ=255)
                                   ggCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAaccgtaccg                                            >  1:621506/1‑44 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCTAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTTGGaa       <  1:574519/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTTGGaa       <  1:3069546/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:1910905/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:333017/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:355019/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:69186/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:1229667/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:828284/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:1208498/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:969727/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:1031222/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2197236/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2264098/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2119027/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:3300585/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2110491/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2898302/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:3210991/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:3097855/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:3330393/67‑1 (MQ=255)
                                                 tgtgTTGCTGATCCTCGCAGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGaa       <  1:2992163/67‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:1264658/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:202576/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:209914/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:92473/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:3267134/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:1918078/68‑1 (MQ=255)
                                                     ttGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:338128/68‑1 (MQ=255)
                                                      tGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATcc  <  1:2022544/67‑1 (MQ=255)
                                                             |                                                           
CGAAATTGTGGCAGTTCATTGCCATTATCGGTATGGGCATTTCAGCGGTTGTGTTGCTGATACTCGCCGAACCGTACCGTATCCGCCGTGTTACCGCATTCTGGAACCCGTGGGAAGATCC  >  minE/83792‑83912

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: