Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,156,422 G→A 45.0% T142T (ACC→ACT obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,156,4220GA45.0% 10.6 / 45.9 60T142T (ACC→ACTobgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (12/15);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCC  >  minE/2156355‑2156490
                                                                   |                                                                    
aCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTgg                                                                     <  1:2780189/69‑1 (MQ=255)
aCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTgg                                                                     <  1:286534/69‑1 (MQ=255)
aCATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTgg                                                                     <  1:211713/69‑1 (MQ=255)
   tACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGttt                                                                  >  1:112205/1‑69 (MQ=255)
   tACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGttt                                                                  >  1:1960508/1‑69 (MQ=255)
   tACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGttt                                                                  >  1:953166/1‑69 (MQ=255)
   tACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGttt                                                                  >  1:940303/1‑69 (MQ=255)
   tACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGttt                                                                  >  1:717045/1‑69 (MQ=255)
                        tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                             >  1:2706274/1‑69 (MQ=255)
                        tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                             >  1:1871465/1‑69 (MQ=255)
                        tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                             >  1:1870657/1‑69 (MQ=255)
                        tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTaa                                             >  1:727381/1‑69 (MQ=255)
                        tCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGGTGGTTTTCTGCCGCGGTGTACGGTTaa                                             >  1:2492928/1‑69 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:180459/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:2441776/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:1518374/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:101426/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:1914522/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:283867/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:2868152/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:2314784/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:3396198/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:987447/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:3076579/65‑1 (MQ=255)
                          aGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGtt                                               <  1:3373689/65‑1 (MQ=255)
                            cTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACgg                                          <  1:2139456/68‑1 (MQ=255)
                            cTCCAGCAGCAGCTCCCGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCcgt                                                          >  1:1095449/1‑51 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTc                                                                >  1:986472/1‑46 (MQ=38)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGtt                                               >  1:1496113/1‑63 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:735507/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:3408871/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:525244/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:989330/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:760635/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:1049433/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:2424796/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAActg                                          >  1:1867218/1‑66 (MQ=255)
                            cTCCAGCAGCAGCTCACGCTTATCGCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAAct                                           >  1:2654424/1‑66 (MQ=255)
                                     cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaa                                >  1:553358/1‑69 (MQ=255)
                                     cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaa                                >  1:3296405/1‑69 (MQ=255)
                                     cagcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGCTGTACGGTTAACGGACGATTTGaa                                >  1:3188072/1‑69 (MQ=255)
                                      agcTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGaaa                               >  1:2267726/1‑69 (MQ=255)
                                          cgcgCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACggg                           <  1:898731/69‑1 (MQ=255)
                                          cgcgCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACggg                           <  1:1316631/69‑1 (MQ=255)
                                                   cGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                    <  1:1250715/67‑1 (MQ=255)
                                                   cGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGcc                    <  1:2073098/67‑1 (MQ=255)
                                                    gCCCTGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2327599/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:764190/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:1075850/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:1753333/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:650685/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2470731/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:193438/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:369676/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2029854/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2194750/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:3339155/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:3007254/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2327774/68‑1 (MQ=255)
                                                    gCCCGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:2424056/68‑1 (MQ=255)
                                                     cccGGCGTACCGTTAGTTTTCTGCCGTGGTGTACGGTTAACTGACGACTTGAAACGGGTATTGCCCa                  <  1:3427035/67‑1 (MQ=255)
                                                         gCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGt            <  1:3355/69‑1 (MQ=255)
                                                         gCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGt            <  1:3064786/69‑1 (MQ=255)
                                                                   ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCg             >  1:1513432/1‑58 (MQ=255)
                                                                   ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc  >  1:1547316/1‑69 (MQ=255)
                                                                   ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc  >  1:1537801/1‑69 (MQ=255)
                                                                   ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgcc  >  1:2429258/1‑69 (MQ=255)
                                                                   ggTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAgccgc   >  1:1355741/1‑68 (MQ=255)
                                                                   |                                                                    
ACATACCGACGTCAGCCAGCAGCATCAGCTCCAGCAGCAGCTCGCGCTTATCGCCCGGCGTGCCGTTGGTTTTCTGCCGCGGTGTACGGTTAACGGACGATTTGAAACGGGTATTGCCCAGACCGTGCCAGCCGCC  >  minE/2156355‑2156490

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: