Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 288,179 T→C 49.1% L275L (TTG→CTG)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE288,1790TC49.1% ‑5.4 / 51.0 58L275L (TTG→CTG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base C (19/9);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGG  >  minE/288113‑288247
                                                                  |                                                                    
gCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:1920474/1‑67 (MQ=255)
gCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:1822801/1‑67 (MQ=255)
gCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:2668172/1‑67 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:2846485/1‑66 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:939897/1‑66 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:88291/1‑66 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:3120060/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:3273467/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:3100277/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:80277/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:1421732/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:3088472/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:1165099/69‑1 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCAGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:3120165/1‑66 (MQ=255)
 cccTGAAGCGCAAAATCGACGCGGCGAAGATGCAGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:3123873/1‑66 (MQ=255)
 cacTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                                                                      >  1:144695/3‑66 (MQ=255)
  ccTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:1833257/68‑1 (MQ=255)
               aTCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGc                                                                   <  1:1768383/55‑1 (MQ=255)
                tCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga                                                    >  1:2819035/1‑69 (MQ=255)
                tCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga                                                    >  1:2746955/1‑69 (MQ=255)
                tCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga                                                    >  1:1724562/1‑69 (MQ=255)
                tCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga                                                    >  1:2544084/1‑69 (MQ=255)
                tCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga                                                    >  1:3405265/1‑69 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:939216/69‑1 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:1140653/69‑1 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:3377502/69‑1 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:457602/69‑1 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:2622142/69‑1 (MQ=255)
                 cGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:2711336/69‑1 (MQ=255)
                   aCGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:426478/67‑1 (MQ=255)
                     gcggcgAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct                                               <  1:967045/69‑1 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:3241164/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:3104029/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:3429103/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:50052/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:51355/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:571457/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:601680/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:649824/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:958979/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1134111/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:3053286/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:2963053/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1657013/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:25864/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:2567893/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1817595/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1822/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1473870/1‑59 (MQ=255)
                                  cGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGCTGCAGAAGCTGAAAATGATGTCTCCg                                            >  1:1937704/1‑59 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:565670/48‑1 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:2933040/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:1702604/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:3409728/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:2121303/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:1210031/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:2259604/1‑69 (MQ=255)
                                      agagGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGc                              >  1:119906/1‑69 (MQ=255)
                                          gCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGGTGAAAATGATGTCTCCGATGTCGGCAGAAGCGAc                           <  1:738161/68‑1 (MQ=255)
                                            aaaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGt                        <  1:1377719/69‑1 (MQ=255)
                                            aaaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGt                        <  1:289283/69‑1 (MQ=255)
                                            aaaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCACAAGCGACCGt                        <  1:477169/69‑1 (MQ=255)
                                             aaaGAGAAAGCGGAAGCAGAGTTTCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGt                        >  1:1414329/1‑68 (MQ=255)
                                             aaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat                                                   <  1:3015154/41‑1 (MQ=255)
                                             aaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTa                       >  1:1746847/1‑69 (MQ=255)
                                             aaaGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTa                       >  1:1557228/1‑69 (MQ=255)
                                                        ggAAGCAGAGGTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTta            <  1:1452149/69‑1 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGt                        <  1:2834976/47‑1 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:1545803/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:1932568/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:2085531/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:2121131/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:3219728/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:2225034/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:1038489/1‑69 (MQ=255)
                                                                  ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg  >  1:978333/1‑69 (MQ=255)
                                                                  |                                                                    
GCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGG  >  minE/288113‑288247

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: