Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,412,631 G→T 27.8% R9L (CGT→CTT)  aslA → acrylsulfatase‑like enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,412,6310GT27.8% 154.1 / 55.0 104R9L (CGT→CTT) aslAacrylsulfatase‑like enzyme
Reads supporting (aligned to +/- strand):  ref base G (24/51);  new base T (19/10);  total (43/61)
Fisher's exact test for biased strand distribution p-value = 3.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 104 total aligned reads displayed.

ATCATTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCC  >  minE/2412572‑2412698
                                                           |                                                                   
atcatTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgt                                                            >  1:268897/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:1078812/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:702036/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:474355/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:125414/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:407231/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:1273810/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:1290853/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:3225222/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:2750468/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:2676818/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCGGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:1856823/1‑69 (MQ=255)
atcatTTTGACGACATGAAACCCGATCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCttgttgt                                                            >  1:3321087/1‑69 (MQ=255)
  catTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:2598946/68‑1 (MQ=255)
  catTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:548301/68‑1 (MQ=255)
  catTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:2216981/68‑1 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:2465973/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:2187271/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:1871497/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:184298/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:1730421/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:1706868/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:1432886/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:855162/1‑67 (MQ=255)
   atTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           >  1:1562884/1‑67 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAGCGTCTTGTTGTTGCt                                                        <  1:2949888/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:267874/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:2516751/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:2097736/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:2381186/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:2311928/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:2954771/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:1486912/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:488686/68‑1 (MQ=255)
     tttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:584050/68‑1 (MQ=255)
      ttGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGt                                                      >  1:3055280/1‑69 (MQ=255)
           gaCATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTcgcc                                                  <  1:438942/68‑1 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgca                                               >  1:3252984/1‑68 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgca                                               >  1:542405/1‑68 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgca                                               >  1:1190030/1‑68 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgca                                               >  1:2535913/1‑68 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgc                                                >  1:842133/1‑68 (MQ=255)
             cATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACTTCTTGTTGTTGCTGTcgccgc                                                >  1:277790/1‑68 (MQ=255)
                 aaaCCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGctc                                           >  1:2641741/1‑69 (MQ=255)
                 aaaCCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGctc                                           >  1:3398492/1‑69 (MQ=255)
                 aaaCCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGctc                                           >  1:2645670/1‑69 (MQ=255)
                 aaaCCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGctc                                           >  1:2351872/1‑69 (MQ=255)
                    ccGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTc                                        <  1:1398450/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:2819101/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:1106162/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:46255/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:426584/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:1408950/69‑1 (MQ=255)
                        aGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCctct                                    <  1:2304140/69‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:2621450/40‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:1697442/40‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:1194324/40‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCttgttgtt                                                           <  1:374733/40‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCt                                                        <  1:610179/43‑1 (MQ=255)
                              ttGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATg                                >  1:1362441/1‑67 (MQ=255)
                                        ttttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:1980313/69‑1 (MQ=255)
                                         tttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:1429789/68‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:890427/67‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:584432/67‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:1898195/67‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:3087348/67‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:2371832/67‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:465792/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:275135/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:1909196/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  >  1:2786528/1‑69 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:355039/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  >  1:2790555/1‑69 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  >  1:3262029/1‑69 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:1822228/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:1767687/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:1755522/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:605748/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  >  1:1073180/1‑69 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:2312476/69‑1 (MQ=255)
                                          ttCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCGCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:2538550/69‑1 (MQ=255)
                                           tCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:2092804/68‑1 (MQ=255)
                                           tCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGa                  <  1:1264356/68‑1 (MQ=255)
                                              ttttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCgcagca                    <  1:2599382/63‑1 (MQ=255)
                                                ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAcccc            >  1:3100164/1‑69 (MQ=255)
                                                ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAccc             >  1:1521256/1‑68 (MQ=255)
                                                ttCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATAcc              >  1:326179/1‑67 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:2195057/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:3205990/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:1550988/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:2539798/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:660999/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:1734001/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:1678766/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:873795/68‑1 (MQ=255)
                                                  cACCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:2257438/68‑1 (MQ=255)
                                                   aCCCAAACTTCTTGTTGTTGCTGTCGCCGCAGCTCTTCCTCTCATAGCCAGCGCAGCAGATACCCCg           <  1:1907313/67‑1 (MQ=255)
                                                     ccAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCa        <  1:570526/68‑1 (MQ=255)
                                                      cAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTACTCTCATGGCCAGCgcagca                    <  1:2856553/55‑1 (MQ=255)
                                                          cGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGcc  <  1:596100/69‑1 (MQ=255)
                                                           |                                                                   
ATCATTTTGATGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTCACCCAAACGTCTTGTTGTTGCTGTCGCCGCCGCTCTTCCTCTCATGGCCAGCGCAGCAGATACCCCGTCAACTGCC  >  minE/2412572‑2412698

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: