Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,976,681 A→G 27.3% intergenic (+9/‑68) 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,6810AG27.3% 73.8 / 29.0 73intergenic (+9/‑68)galP/sprTD‑galactose transporter/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/20);  new base G (5/15);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 7.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

TTGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCC  >  minE/1976617‑1976745
                                                                |                                                                
ttGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCtt                                                              >  1:2271546/1‑69 (MQ=255)
ttGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCtt                                                              >  1:349405/1‑69 (MQ=255)
ttGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCt                                                               >  1:1813095/1‑68 (MQ=255)
ttGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCt                                                               <  1:2675708/68‑1 (MQ=255)
ttGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCt                                                               <  1:749882/68‑1 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGTGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:1695242/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCt                                                               >  1:3402828/1‑66 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:11067/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:2379411/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:2833735/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:3004512/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:1414337/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:878512/1‑69 (MQ=255)
  gAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcc                                                            >  1:1191703/1‑69 (MQ=255)
     cGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcctcc                                                         >  1:3414102/1‑69 (MQ=255)
     cGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcctcc                                                         >  1:1551697/1‑69 (MQ=255)
     cGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTtcctcc                                                         >  1:1174241/1‑69 (MQ=255)
           cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                    >  1:953158/1‑68 (MQ=255)
           cTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATcg                                                    >  1:1161568/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:2981291/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:3043593/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:2472992/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:3293021/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:1454502/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:2690819/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:2332509/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:1641400/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:2154429/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:206264/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:645463/69‑1 (MQ=255)
                 aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga                                             <  1:1167390/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:1548149/69‑1 (MQ=255)
                       cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                       <  1:2625165/69‑1 (MQ=255)
                        gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                       <  1:2572398/68‑1 (MQ=255)
                        gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc                                       <  1:2456886/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:3144606/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)
                                      atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                        >  1:934969/4‑69 (MQ=255)
                                      atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                        >  1:2533900/4‑69 (MQ=255)
                                      atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                        >  1:3290564/4‑69 (MQ=255)
                                      atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                        >  1:1234768/4‑69 (MQ=255)
                                      atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt                        >  1:596624/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)
                                                                |                                                                
TTGAACGTAATCTGATGAAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCC  >  minE/1976617‑1976745

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: