Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,309,857 A→G 23.7% A80A (GCA→GCG hslV → peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,309,8570AG23.7% 133.1 / 35.3 105A80A (GCA→GCGhslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base G (12/13);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

ACTGTTTGAACGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGAA  >  minE/2309794‑2309920
                                                               |                                                               
aCTGTTTGAACGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAg                                                             <  1:2339015/68‑1 (MQ=255)
aCTGTTTGAACGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAg                                                             <  1:632788/68‑1 (MQ=255)
aCTGTTTGAACGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAg                                                             <  1:2969601/68‑1 (MQ=255)
         aCGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTa                                                    <  1:440359/68‑1 (MQ=255)
         aCGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTa                                                    <  1:777044/68‑1 (MQ=255)
         aCGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTa                                                    <  1:19223/68‑1 (MQ=255)
           gTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCg                                                 >  1:1386432/1‑69 (MQ=255)
                cTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCg                                             <  1:307670/68‑1 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:2276169/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:458276/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:2443452/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:1326689/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:2361269/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgctg                                      >  1:663286/1‑68 (MQ=255)
                       tGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGAATGGCGTACCGATCGCAtgctg                                      >  1:1247131/1‑68 (MQ=255)
                       tGCATAAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCAtgct                                       >  1:3347140/1‑67 (MQ=255)
                         cATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCg                                                      >  1:2384075/1‑50 (MQ=255)
                         cATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGt                                                     >  1:2002093/1‑51 (MQ=255)
                         cATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGt                                                     >  1:246686/1‑51 (MQ=255)
                         cATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGt                                                     >  1:1761509/1‑51 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCaaa                                >  1:518945/1‑68 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAAc                               >  1:465672/1‑69 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAAc                               >  1:1668713/1‑69 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAAc                               >  1:1708502/1‑69 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAAc                               >  1:1348056/1‑69 (MQ=255)
                             aGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGAGCAAAc                               >  1:3423623/1‑69 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:1818987/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:433552/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:2253636/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:1305084/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:3262712/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:503419/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:276021/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:2763382/69‑1 (MQ=255)
                                 ccATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:1610642/69‑1 (MQ=255)
                                    tCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGc                        <  1:1355972/69‑1 (MQ=255)
                                    tCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGc                        <  1:1248861/69‑1 (MQ=255)
                                      tGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCTTACCGATCGCATGCTGCGCaaa                                <  1:2656221/59‑1 (MQ=255)
                                      tGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCTATCGCATGCTGCGCaaa                                <  1:289564/59‑1 (MQ=255)
                                      tGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCaaa                                <  1:178475/59‑1 (MQ=255)
                                      tGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCaaa                                <  1:1943905/59‑1 (MQ=255)
                                      tGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCaaa                                <  1:2941024/59‑1 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:3204488/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:2373223/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:2301393/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:1014644/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:2292628/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:1164234/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:1187888/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:2016712/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgc                   >  1:1515894/1‑69 (MQ=255)
                                         tCAAAGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActg                    >  1:2340317/1‑68 (MQ=255)
                                           aaaGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgct                  >  1:3427317/1‑68 (MQ=255)
                                           aaaGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgct                  >  1:3372591/1‑68 (MQ=255)
                                           aaaGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgct                  >  1:2956915/1‑68 (MQ=255)
                                           aaaGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCActgct                  >  1:1434887/1‑68 (MQ=255)
                                             aGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTgg                <  1:3290839/68‑1 (MQ=255)
                                             aGCCGCCGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTgg                <  1:1519425/68‑1 (MQ=255)
                                              gccgccGTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTgg                <  1:1645389/67‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:765918/49‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:156751/49‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:1574748/49‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:1647069/49‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:637285/49‑1 (MQ=255)
                                                  ccgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:180451/49‑1 (MQ=255)
                                                   cgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:2738659/48‑1 (MQ=255)
                                                   cgTTGAGCTGGCGAAAGACTGGCGTACCGATCGCATGCTGCGCAAACt                              <  1:458768/48‑1 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        <  1:2312484/69‑1 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:1699299/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:2350309/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:403050/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:355027/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:3432808/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        <  1:1189843/69‑1 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:3324530/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:2436564/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:1459720/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:1845597/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:3071065/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:2556755/1‑69 (MQ=255)
                                                    gTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        >  1:692336/1‑69 (MQ=255)
                                                     ttGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTcgcg        <  1:1388575/68‑1 (MQ=255)
                                                       gAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGc                        >  1:2367825/1‑50 (MQ=255)
                                                        aGCTGGCGAAAGACTGGCGTACCGAGCGCATGCTGCGCAAACTTGAAGCAGTGCTgg                <  1:2367276/57‑1 (MQ=255)
                                                        aGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGa                           <  1:2013339/46‑1 (MQ=255)
                                                          cTGTCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:2086150/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGaa  >  1:3371829/1‑69 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGaa  >  1:1171368/1‑69 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:843861/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:928661/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:3094812/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:2875197/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:2771109/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:2176979/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   >  1:1549738/1‑68 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:1510223/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:1242940/68‑1 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCAATGCTGGCAGTCGCGGATGaa  >  1:1071784/1‑69 (MQ=255)
                                                          cTGGCAAAAGACTGGCGTACCGATCGCATGATGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGaa  >  1:2514296/1‑69 (MQ=255)
                                                           tGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGa   <  1:1733278/67‑1 (MQ=255)
                                                               |                                                               
ACTGTTTGAACGTAAACTGGAAATGCATCAGGGCCATCTGGTCAAAGCCGCCGTTGAGCTGGCAAAAGACTGGCGTACCGATCGCATGCTGCGCAAACTTGAAGCACTGCTGGCAGTCGCGGATGAA  >  minE/2309794‑2309920

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: