Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,545,265 A→G 5.7% F78L (TTT→CTT)  rfaC ← ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,545,2650AG 5.7% 372.5 / 9.1 122F78L (TTT→CTT) rfaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base A (35/80);  new base G (2/5);  total (37/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 127 total aligned reads displayed.

CCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAT  >  minE/2545196‑2545334
                                                                     |                                                                     
ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                                      <  1:1621691/71‑1 (MQ=255)
ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                                      <  1:4344104/71‑1 (MQ=255)
ccgcGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                                      <  1:3979380/71‑1 (MQ=255)
ccgcGCTTTTTACAAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAAc                                                                      >  1:5088682/1‑71 (MQ=255)
 cgcgCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACg                                                                     <  1:2466778/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:2783851/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:1150426/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:1198068/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:498238/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:2125700/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:2036999/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:3759455/71‑1 (MQ=255)
     cTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:3570713/71‑1 (MQ=255)
      tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                                   <  1:1000565/68‑1 (MQ=255)
      tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                                   <  1:4257384/68‑1 (MQ=255)
      tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                                   <  1:4336245/68‑1 (MQ=255)
      tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                                   <  1:1922041/68‑1 (MQ=255)
      tttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCt                                                                   <  1:4965671/68‑1 (MQ=255)
       ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                               >  1:3616911/1‑71 (MQ=255)
       ttttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGc                                                               >  1:4708153/1‑71 (MQ=255)
        tttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGTTTCGTGAAGCGCtac                                                                 >  1:1482223/1‑66 (MQ=37)
        tttACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCg                                                              >  1:4605049/1‑71 (MQ=255)
            ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                                           <  1:550479/70‑1 (MQ=255)
            ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                                           <  1:2505468/70‑1 (MQ=255)
            ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                                           <  1:4963426/70‑1 (MQ=255)
            ccAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTc                                                           <  1:3365253/70‑1 (MQ=255)
                  ccTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtttt                                                    <  1:4574140/71‑1 (MQ=255)
                  ccTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtttt                                                    <  1:3628929/71‑1 (MQ=255)
                    tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:4996978/56‑1 (MQ=255)
                    tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:2176097/56‑1 (MQ=255)
                    tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:4081037/56‑1 (MQ=255)
                    tGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCttt                                                                 <  1:1905626/56‑1 (MQ=255)
                     gAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                                 >  1:2790271/1‑71 (MQ=255)
                     gAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                                 >  1:3432281/1‑71 (MQ=255)
                      aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                                  <  1:3329540/69‑1 (MQ=255)
                      aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                                  <  1:2647311/69‑1 (MQ=255)
                      aGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                                 <  1:542549/70‑1 (MQ=255)
                            gATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGcc                                           >  1:2596579/1‑70 (MQ=255)
                              tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        <  1:2421958/71‑1 (MQ=255)
                              tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        <  1:3017815/71‑1 (MQ=255)
                              tAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        <  1:4275101/71‑1 (MQ=255)
                               aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTAt                                                  >  1:3659309/1‑60 (MQ=255)
                               aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATg                                                 >  1:3548593/1‑61 (MQ=255)
                               aaCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTAGCGTTCCGCTTTTAt                                                  >  1:4556396/1‑60 (MQ=255)
                                aCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGTTTTTATGGGGGCCGAGaa                                      <  1:3869313/71‑1 (MQ=255)
                                aCTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGaa                                      <  1:1115527/71‑1 (MQ=255)
                                      tCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                                     <  1:3859431/50‑1 (MQ=255)
                                        aTAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtt                                                      >  1:922801/1‑47 (MQ=255)
                                        aTAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCtt                                                      >  1:1462039/1‑47 (MQ=255)
                                         tAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              >  1:2320381/1‑70 (MQ=255)
                                          aGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            >  1:992030/1‑71 (MQ=255)
                                          aGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            >  1:1843934/1‑71 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                                     <  1:4479432/44‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                                     <  1:2736096/44‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCttt                                                     <  1:3931862/44‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                                   <  1:5152062/46‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                                   <  1:2908587/46‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTa                                                   <  1:1630178/46‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCATGCTTTa                            <  1:8017/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:129987/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:1980312/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:2244955/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:4216265/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:4158796/69‑1 (MQ=255)
                                            ttCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            <  1:4720893/69‑1 (MQ=255)
                                              ctctGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGc                         >  1:1103572/1‑70 (MQ=255)
                                                  gCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGc                     >  1:863079/1‑70 (MQ=255)
                                                   cTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCg                    <  1:4155942/70‑1 (MQ=255)
                                                   cTTGTAGCGCTTCACGAAACGATTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCg                    <  1:4631338/70‑1 (MQ=255)
                                                      gTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg                 >  1:2249930/1‑70 (MQ=255)
                                                      gTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg                 <  1:889786/70‑1 (MQ=255)
                                                       tAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt                <  1:1387420/70‑1 (MQ=255)
                                                       tAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt                <  1:4423790/70‑1 (MQ=255)
                                                        aGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGt                <  1:1412717/69‑1 (MQ=255)
                                                          cgcTTCACGAAACGCTTTGCGTTCCTCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCTGCGACg                 <  1:4195220/66‑1 (MQ=255)
                                                            cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:4403758/51‑1 (MQ=255)
                                                            cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:2479647/51‑1 (MQ=255)
                                                            cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:4659195/51‑1 (MQ=255)
                                                            cTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGCCCGAGAACCAGGCtt                              <  1:4432958/51‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:5065098/50‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:4027466/50‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:25835/50‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:865078/50‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg                 <  1:615741/63‑1 (MQ=255)
                                                             ttCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACg                 <  1:3004122/63‑1 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGcc                                           >  1:1089542/1‑36 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        >  1:1082024/1‑39 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        >  1:1209617/1‑39 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        >  1:2086242/1‑39 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCgag                                        >  1:2954640/1‑39 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAgg                                 >  1:1387069/1‑46 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            >  1:3884007/1‑51 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            >  1:2674269/1‑51 (MQ=255)
                                                              tCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTa                            >  1:3322721/1‑51 (MQ=255)
                                                                 cGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCtt                              <  1:1149178/46‑1 (MQ=255)
                                                                    aaCGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAt  <  1:1605771/71‑1 (MQ=255)
                                                                    aaCGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAt  <  1:583885/71‑1 (MQ=255)
                                                                     aCGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAt  >  1:1745758/1‑70 (MQ=255)
                                                                     aCGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAt  >  1:653097/1‑70 (MQ=255)
                                                                     aCGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAt  >  1:3218194/1‑70 (MQ=255)
                                                                     |                                                                     
CCGCGCTTTTTACCAGCCCCTGAGCGTCGATAACTGCGTCATAGTTCTCTGCTTGTAGCGCTTCACGAAACGCTTTGCGTTCCGCTTTTATGGGGGCCGAGAACCAGGCTTTACGCCAGCGACGTATTGCCACAGGAAT  >  minE/2545196‑2545334

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: