Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,976,695 G→A 36.3% intergenic (+23/‑54) galP → / → sprT D‑galactose transporter/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,976,6950GA36.3% 33.3 / 30.9 55intergenic (+23/‑54)galP/sprTD‑galactose transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/20);  new base A (5/15);  total (20/35)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAC  >  minE/1976628‑1976756
                                                                   |                                                             
cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                               >  1:1161568/1‑68 (MQ=255)
cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                               >  1:953158/1‑68 (MQ=255)
cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                               >  1:2020759/1‑68 (MQ=255)
cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCAAGCTTCCTCCCATcg                                                               >  1:446394/1‑67 (MQ=255)
 tgatgaAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                               >  1:1295191/1‑67 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2289386/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2332509/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2154429/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2472992/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:206264/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:1167390/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2690819/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2981291/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:3043593/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:3293021/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:1641400/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:1454502/69‑1 (MQ=255)
      aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:645463/69‑1 (MQ=255)
          gTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGcc                                                    <  1:3321251/69‑1 (MQ=255)
            cGTAAACTGCGCGAAATATGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:2649210/69‑1 (MQ=255)
            cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:3436286/69‑1 (MQ=255)
            cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:354378/69‑1 (MQ=255)
            cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:2625165/69‑1 (MQ=255)
            cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:1548149/69‑1 (MQ=255)
             gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:2456886/68‑1 (MQ=255)
             gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:2572398/68‑1 (MQ=255)
               aaaCTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:1949984/66‑1 (MQ=255)
                     cgcgAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:3033371/54‑1 (MQ=255)
                     cgcgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:190218/60‑1 (MQ=255)
                       cgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTc                                       >  1:575430/1‑69 (MQ=255)
                       cgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTc                                       >  1:182946/1‑69 (MQ=255)
                         aaaTCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                                        <  1:2214146/50‑1 (MQ=255)
                         aaaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCAt                                     <  1:3381842/69‑1 (MQ=255)
                         aaaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCAt                                     <  1:1748310/69‑1 (MQ=255)
                          aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc                                    <  1:247650/69‑1 (MQ=255)
                          aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc                                    <  1:3259295/69‑1 (MQ=255)
                          aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc                                    <  1:2460781/69‑1 (MQ=255)
                          aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc                                    <  1:1412990/69‑1 (MQ=255)
                          aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc                                    <  1:3144606/69‑1 (MQ=255)
                           atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                                   >  1:2533900/4‑69 (MQ=255)
                           atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                                   >  1:596624/4‑69 (MQ=255)
                           atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                                   >  1:3290564/4‑69 (MQ=255)
                           atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                                   >  1:934969/4‑69 (MQ=255)
                           atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                                   >  1:1234768/4‑69 (MQ=255)
                               gcgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                                  <  1:1660064/50‑1 (MQ=255)
                                cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttct                              >  1:628232/1‑69 (MQ=255)
                                cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttct                              >  1:2857005/1‑69 (MQ=255)
                                cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttc                               >  1:1540356/1‑68 (MQ=255)
                                          tAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATc                    >  1:2868335/1‑69 (MQ=255)
                                            aTCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCct                  >  1:3018651/1‑69 (MQ=255)
                                               tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCctct                >  1:3073049/1‑68 (MQ=255)
                                               tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCctct                >  1:76762/1‑68 (MQ=255)
                                               tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTg               >  1:2345227/1‑69 (MQ=255)
                                               tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTg               >  1:3232412/1‑69 (MQ=255)
                                                 cccAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGcc             >  1:1583595/1‑69 (MQ=255)
                                                 cccAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGcc             >  1:1836475/1‑69 (MQ=255)
                                                         tcctccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaa     >  1:476875/1‑69 (MQ=255)
                                                         tcctccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaa     >  1:2114436/1‑69 (MQ=255)
                                                            tccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaaaa   <  1:3342548/68‑1 (MQ=255)
                                                             cccATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAc  <  1:1018433/68‑1 (MQ=255)
                                                                   |                                                             
CTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAC  >  minE/1976628‑1976756

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: