Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,781,789 G→A 57.7% E80E (GAG→GAA yfjD → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,781,7890GA57.7% 12.3 / 37.2 52E80E (GAG→GAAyfjDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base A (4/26);  total (14/38)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGTTTTCTGCTGGC  >  minE/1781722‑1781857
                                                                   |                                                                    
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:2016217/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:1259206/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:535957/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:382211/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:1740080/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:3331871/69‑1 (MQ=255)
tGTACGGCGATGCGGGCGTGGCAATTGAGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAgg                                                                     <  1:2555262/69‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:3272812/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2077059/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2324799/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2729237/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1002544/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2796098/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2796797/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2912985/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2940175/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:2759978/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:3303441/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:3373867/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:48408/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:705083/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:708917/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:926919/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:980345/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1569957/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1051553/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1171893/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1462760/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1543873/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1548315/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:170288/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1727498/66‑1 (MQ=255)
   atggCGATGCTGGCGTGGCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTa                                                                   <  1:1816740/66‑1 (MQ=255)
    cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                                                  <  1:2908020/68‑1 (MQ=255)
    cGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                                                  <  1:3166493/68‑1 (MQ=255)
              ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                                                        <  1:2189055/68‑1 (MQ=255)
              ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                                                        <  1:1895410/68‑1 (MQ=255)
              ggCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                                                        <  1:458224/68‑1 (MQ=255)
                   ggCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTAt                                                                  <  1:70601/53‑1 (MQ=255)
                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                                                  >  1:1728477/1‑69 (MQ=255)
                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                                                  >  1:140396/1‑69 (MQ=255)
                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                                                  >  1:3147088/1‑69 (MQ=255)
                   ggCAATTGCGACTGGCGTGCTGACCTTCGTAGTGCTGGTTTTTGCTGAAGTATTGCCGAAAACCATTGc                                                  >  1:1980711/1‑69 (MQ=255)
                          gCGACTGGTGTGCTGACTTTTGTCTTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTg                                           >  1:1721356/1‑69 (MQ=255)
                          gCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTg                                           >  1:1620782/1‑69 (MQ=255)
                               tGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                                                        <  1:3401566/51‑1 (MQ=255)
                                    tgCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAAc                                                        <  1:790615/46‑1 (MQ=255)
                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTa                                         >  1:177341/1‑59 (MQ=255)
                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa                               >  1:3386557/1‑69 (MQ=255)
                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa                               >  1:1465614/1‑69 (MQ=255)
                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa                               >  1:918640/1‑69 (MQ=255)
                                      cTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGaaaaa                               >  1:2975779/1‑69 (MQ=255)
                                           ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg                           <  1:2703448/68‑1 (MQ=255)
                                           ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg                           <  1:370971/68‑1 (MQ=255)
                                           ttttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCg                           <  1:1197884/68‑1 (MQ=255)
                                            tttGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCt                         >  1:3085011/1‑69 (MQ=255)
                                                          aTTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGttt           <  1:3134857/69‑1 (MQ=255)
                                                           tttGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGtttt          >  1:659003/1‑69 (MQ=255)
                                                             tGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGTTTTct        >  1:1637279/1‑69 (MQ=255)
                                                                   ggTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGTTTTCTGCTGGc  >  1:1440909/1‑69 (MQ=255)
                                                                   |                                                                    
TGTACGGCGATGCGGGCGTGGCAATTGCGACTGGTGTGCTGACTTTTGTCGTACTGGTATTTGCTGAGGTATTGCCGAAAACCATTGCCGCGCTGTACCCGGAAAAAGTCGCTTATCCGAGTAGTTTTCTGCTGGC  >  minE/1781722‑1781857

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: