Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,198,000 T→C 40.5% I173I (ATT→ATC degS → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,198,0000TC40.5% 27.1 / 52.2 74I173I (ATT→ATCdegSserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base C (24/6);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 3.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 111 total aligned reads displayed.

ATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTC  >  minE/2197932‑2198065
                                                                    |                                                                 
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGTattatt                                                                   <  1:179/69‑1 (MQ=255)
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1966793/69‑1 (MQ=255)
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:2382658/69‑1 (MQ=255)
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1253792/69‑1 (MQ=255)
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:834280/69‑1 (MQ=255)
aTTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:174706/69‑1 (MQ=255)
   ggCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGt                                                                >  1:1535400/1‑69 (MQ=255)
   ggCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGt                                                                >  1:1110421/1‑69 (MQ=255)
   ggCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGt                                                                >  1:73120/1‑69 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:3242785/54‑1 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1473397/54‑1 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:3112975/54‑1 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGa                                                    >  1:1643532/1‑69 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGa                                                    >  1:673289/1‑69 (MQ=255)
               cTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGa                                                    >  1:2514891/1‑69 (MQ=255)
               cTGGCGATCGGTAACCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1427077/53‑1 (MQ=255)
                 ggCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1895266/52‑1 (MQ=255)
                 ggCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:2515452/52‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1841831/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1793516/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:3079592/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:2006425/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:25539/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:17404/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:379174/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:123319/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:2538957/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1628720/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1620686/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:3386204/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:2326805/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:735825/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:3083626/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:119330/46‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:272403/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:2970130/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:2131028/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:185854/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:1186531/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:1215016/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:1220811/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTc                                             >  1:1046533/1‑68 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGTCt                                            <  1:1096012/69‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGTCt                                            <  1:2074125/69‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGTCt                                            <  1:3275446/69‑1 (MQ=255)
                       cGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGTCt                                            <  1:2760663/69‑1 (MQ=255)
                        ggTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGattatt                                                                   <  1:1592260/45‑1 (MQ=255)
                          tAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGa                                          <  1:880829/68‑1 (MQ=255)
                          tAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGa                                          <  1:2079132/68‑1 (MQ=255)
                          tAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGa                                          <  1:2512971/68‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:2174969/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:1265875/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:729011/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:1552975/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:2171112/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:659826/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:553477/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:395004/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:1997663/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:1975004/67‑1 (MQ=255)
                                            gCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCaa                         <  1:385977/67‑1 (MQ=255)
                                             cAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAAc                      <  1:1910023/69‑1 (MQ=255)
                                             cAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAAc                      <  1:1512214/69‑1 (MQ=255)
                                              aGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAAc                      <  1:2355372/68‑1 (MQ=255)
                                                  cATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTtcc                  >  1:894686/1‑68 (MQ=255)
                                                  cATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTtcc                  >  1:3002931/1‑68 (MQ=255)
                                                  cATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTtcc                  >  1:2310947/1‑68 (MQ=255)
                                                  cATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTtcc                  >  1:1284323/1‑68 (MQ=255)
                                                  cATTACCAAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTtcc                  >  1:2737458/1‑68 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACTGGTCGAATCGGt                                              >  1:2966829/1‑38 (MQ=38)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATc                                                 >  1:178167/1‑35 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCgg                                               >  1:2336580/1‑37 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:2700503/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:938106/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1072050/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:855677/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:765869/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1109954/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1310667/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1314609/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1645826/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:25187/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:2783375/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:2926066/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:2651389/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:263246/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:304424/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:252372/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:3378017/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:1032042/1‑38 (MQ=255)
                                                    ttACCCAGGGGATTATCAGTGCCACGGGTCGAATCGGt                                              >  1:333300/1‑38 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:1671637/68‑1 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:2586842/68‑1 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:403631/68‑1 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:1103180/68‑1 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:1889442/68‑1 (MQ=255)
                                                        ccAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCaaa            <  1:1098990/68‑1 (MQ=255)
                                                                ttattaGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGCGGCAAAACTTCCTCCAAACCGATGCTTc  <  1:2138509/69‑1 (MQ=255)
                                                                 tattaGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTc  <  1:480920/69‑1 (MQ=255)
                                                                 tattaGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTc  <  1:2033377/69‑1 (MQ=255)
                                                                    |                                                                 
ATTGGCGACGTAGTACTGGCGATCGGTAACCCGTACAACCTCGGGCAGACCATTACCCAGGGGATTATTAGTGCCACGGGTCGAATCGGTCTGAACCCGACCGGGCGGCAAAACTTCCTCCAAACCGATGCTTC  >  minE/2197932‑2198065

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: