Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 161,135 C→T 21.1% intergenic (+2/‑153) degP → / → cdaR serine endoprotease (protease Do), membrane‑associated/DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE161,1350CT21.1% 200.3 / 37.9 136intergenic (+2/‑153)degP/cdaRserine endoprotease (protease Do), membrane‑associated/DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (76/31);  new base T (6/23);  total (82/54)
Fisher's exact test for biased strand distribution p-value = 2.37e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 148 total aligned reads displayed.

GCAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCA  >  minE/161066‑161203
                                                                     |                                                                    
gCAAACCGTCTGTCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAAtc                                                                      <  1:1194525/69‑1 (MQ=255)
 cAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAAtc                                                                      <  1:2574641/69‑1 (MQ=255)
 cAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAAtc                                                                      <  1:1514779/69‑1 (MQ=255)
      cGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTccc                                                                  >  1:2587637/1‑68 (MQ=255)
      cGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTccc                                                                  >  1:1332033/1‑68 (MQ=255)
                  cACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTcc                                                      >  1:1051376/1‑68 (MQ=255)
                  cACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTcc                                                      >  1:1122356/1‑68 (MQ=255)
                  cACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTcc                                                      >  1:124439/1‑68 (MQ=255)
                  cACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTcc                                                      >  1:892915/1‑68 (MQ=255)
                  cACTCAACATTCAGCGCGGCGACAGCACCATATACCTGTTAATGCAGTAATCTCCATCAACCCCTTcc                                                      >  1:2264615/1‑68 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTc                                                                >  1:1931656/1‑54 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTc                                                                >  1:1695419/1‑54 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTc                                                                >  1:1532941/1‑54 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTc                                                                >  1:786462/1‑54 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTc                                                                >  1:2848433/1‑54 (MQ=255)
                      cAACATTCAGCGCGGCGACAGCACCATCTACCTGATAATGCAGTAATCTCCCTc                                                                >  1:1028147/1‑54 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:932696/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:220364/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:1593242/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:3208150/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:20345/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:1761728/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:3419909/69‑1 (MQ=255)
                           ttCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            <  1:3368176/69‑1 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAAc                                               >  1:2375318/1‑64 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:2387405/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:2887928/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:3269256/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:559469/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:71805/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:136360/1‑67 (MQ=255)
                             cAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:1253156/1‑67 (MQ=255)
                             cAGCGCGGCGAAAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACggg                                            >  1:1402505/1‑67 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:3273795/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:3407605/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:375460/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:993631/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:56385/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:2636547/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:829481/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:944893/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:980965/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:2468646/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:1282747/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:1290547/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:1845454/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:1203353/68‑1 (MQ=255)
                               gcgcGGCGACAGCACCATCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:189843/68‑1 (MQ=255)
                                cgcgGCGACAGCACCATCTTCCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAg                                         <  1:2875712/67‑1 (MQ=255)
                                            aCCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTAc                            >  1:2362864/1‑68 (MQ=255)
                                            aCCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTAc                            >  1:1129993/1‑68 (MQ=255)
                                            aCCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACa                           >  1:2539572/1‑69 (MQ=255)
                                            aCCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACa                           >  1:3120628/1‑69 (MQ=255)
                                            aCCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACa                           >  1:1721905/1‑69 (MQ=255)
                                               aTCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATc                        <  1:240230/69‑1 (MQ=255)
                                               aTCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATc                        <  1:1879875/69‑1 (MQ=255)
                                               aTCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATc                        <  1:911731/69‑1 (MQ=255)
                                               aTCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATc                        <  1:3233612/69‑1 (MQ=255)
                                                tCTACCTGTTAATGCAGTAATTTCCGTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAAt                         >  1:1590198/1‑67 (MQ=255)
                                                tCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCt                       >  1:2472658/1‑69 (MQ=255)
                                                tCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCt                       >  1:2923317/1‑69 (MQ=255)
                                                tCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCt                       >  1:2215290/1‑69 (MQ=255)
                                                tCTACCTGTTAATGCAGTAATTTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCt                       >  1:2784830/1‑69 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGt                                     >  1:2528645/1‑53 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:2089241/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:1877341/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:3031763/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:3052227/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:231445/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:664328/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:3161257/1‑55 (MQ=255)
                                                  tACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTtc                                   >  1:1537727/1‑55 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgtg                    >  1:550318/1‑66 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgtg                    >  1:3200839/1‑66 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgtg                    >  1:310471/1‑66 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgtg                    >  1:2946964/1‑66 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgtg                    >  1:2761753/1‑66 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCtgt                     >  1:3422028/1‑65 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGa                   <  1:1470108/67‑1 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGa                   <  1:829207/67‑1 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGa                   <  1:2758163/67‑1 (MQ=255)
                                                      tGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGa                   <  1:1633189/67‑1 (MQ=255)
                                                              cAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTcaccac         <  1:3097566/69‑1 (MQ=255)
                                                              cAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTcaccac         <  1:1306583/69‑1 (MQ=255)
                                                              cAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACATcaccac         <  1:2768675/69‑1 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:107433/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:515253/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:3430700/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:1994568/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:2974341/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:2920013/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:282444/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:263783/1‑69 (MQ=255)
                                                                   aTCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGGTCTCCTTACAATCTGTGAACTTCACCACAACTc    >  1:2474168/1‑69 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATc                        <  1:2467849/47‑1 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACTACTCCa  <  1:437123/69‑1 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCa  <  1:2082406/69‑1 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCa  <  1:3182586/69‑1 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCa  <  1:3174978/69‑1 (MQ=255)
                                                                     ctcCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCa  <  1:26727/69‑1 (MQ=255)
                                                                     |                                                                    
GCAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGCAGTAATCTCCCTCAACCCCTTCCTGAAAACGGGAAGGGGTTCTCCTTACAATCTGTGAACTTCACCACAACTCCA  >  minE/161066‑161203

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: