Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,985,434 A→G 8.9% intergenic (+37/+118) yggW → / ← yggM predicted oxidoreductase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,985,4340AG 8.9% 176.1 / 7.9 67intergenic (+37/+118)yggW/yggMpredicted oxidoreductase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base G (4/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTG  >  minE/1985373‑1985498
                                                             |                                                                  
ataagacgcgacagcgtcgcatcaggcattgattgccggatgcggcgtaAACGCCTTATCCGGCCTACAt                                                            >  1:3573453‑M1/50‑70 (MQ=255)
ataagacgcgacagcgtcgcatcaggcattgattgccggatgcggcgtaAACGCCTTATCCGGCCTACAt                                                            >  1:3568879‑M1/50‑70 (MQ=255)
ataagacgcgacagcgtcgcatcaggcattgattgccggatgcggcacaAGTGCCTTATCCGGCCTacagg                                                           <  1:2524437‑M1/22‑3 (MQ=255)
ataagacgcgacagcgtcgcatcaggcattgattgccggatgcggcacaAGTGCCTTATCCGGCCTacagg                                                           <  1:3166230‑M1/22‑3 (MQ=255)
 ctgGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACAtg                                                           <  1:474033/70‑1 (MQ=255)
    gAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACAtgtg                                                        <  1:3066371/69‑1 (MQ=255)
    gAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTgg                                                       >  1:2307832/1‑70 (MQ=255)
      gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAg                                                    >  1:3183642/1‑71 (MQ=255)
         ttttCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCgg                                                 >  1:70812/1‑71 (MQ=255)
          tttCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCg                                                  <  1:966476/69‑1 (MQ=255)
            tCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTg                                             >  1:1838983/1‑71 (MQ=255)
            tCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTg                                             >  1:215397/1‑71 (MQ=255)
             cTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTg                                             <  1:1729025/70‑1 (MQ=255)
             cTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTg                                             <  1:469180/70‑1 (MQ=255)
                gCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAgg                                         <  1:2617255/71‑1 (MQ=255)
                 cTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGt                                        <  1:3564711/71‑1 (MQ=255)
                  tGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTc                                       >  1:2808273/1‑71 (MQ=255)
                  tGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTc                                       >  1:1480362/1‑71 (MQ=255)
                   gAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCt                                      >  1:1761731/1‑71 (MQ=255)
                    aGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCt                                      >  1:3269423/1‑70 (MQ=255)
                    aGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCt                                      <  1:2200282/70‑1 (MQ=255)
                        aaCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACAt                                                            >  1:748022/1‑46 (MQ=255)
                        aaCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATa                                  <  1:680132/70‑1 (MQ=255)
                         aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACAtgt                                                          <  1:243240/47‑1 (MQ=255)
                         aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGAt                                   <  1:619715/68‑1 (MQ=255)
                         aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAg                                <  1:3500078/71‑1 (MQ=255)
                          cTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGa                               <  1:1531311/71‑1 (MQ=255)
                           ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAc                              >  1:3173735/1‑71 (MQ=255)
                            aagaTTACCGGATGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGAc                              <  1:1433123/68‑1 (MQ=2)
                             gTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcgc                            >  1:1665548/1‑71 (MQ=255)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGt                                        <  1:3751649/55‑1 (MQ=255)
                                 tGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGAGGTTGTAGGt                                        <  1:596031/55‑1 (MQ=255)
                                  gCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATaa                                 >  1:730105/1‑61 (MQ=255)
                                  gCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcgc                            >  1:841225/1‑66 (MQ=255)
                                  gCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcgc                            >  1:1109265/1‑66 (MQ=255)
                                  gCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAg                       >  1:1522311/1‑71 (MQ=255)
                                   ccGGATGCGGCGTGAACGCCTTATCCAGCCGACTTGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGc                      <  1:2470911/71‑1 (MQ=255)
                                   ccGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCaa                        <  1:2150778/69‑1 (MQ=255)
                                    cggatgcggcgtaAACGCCTTATCCGGCCTACATTTGACAGCCGTTGTAGGCCTGATAAGACGCGCAAGCg                     >  1:455237‑M1/14‑71 (MQ=255)
                                    cGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGc                      <  1:2905422/70‑1 (MQ=255)
                                     ggATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCa                         <  1:3719056/66‑1 (MQ=255)
                                     ggATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCa                         <  1:983599/66‑1 (MQ=255)
                                     ggATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGc                      >  1:2160819/1‑69 (MQ=255)
                                      gATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTc                   >  1:3877854/1‑71 (MQ=255)
                                      gATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTc                   >  1:3154714/1‑71 (MQ=255)
                                        tGCGGCGTGAACGTCTTATCCACCCGACAGGTCAGGCGCGGTCTGTAGGCTGATAAGAcg                             <  1:1481298/60‑1 (MQ=2)
                                         gcggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCa                >  1:3812378/1‑71 (MQ=255)
                                           ggcgTGAACGCCTTATCCAGCCGACATGTGGCAGc                                                   <  1:1306985/35‑1 (MQ=255)
                                              gTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTc                                       <  1:1074741/43‑1 (MQ=255)
                                              gTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTc                                       <  1:3264685/43‑1 (MQ=255)
                                                 aaCGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATa                                  <  1:2127416/45‑1 (MQ=255)
                                                 aaCGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAgtgt                    >  1:2504970/1‑59 (MQ=255)
                                                 aaCGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAgtgt                    >  1:2500351/1‑59 (MQ=255)
                                                 aaCGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGt                    >  1:3255656/1‑59 (MQ=255)
                                                  aCGCCTTATCCATCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCaa                        <  1:2023994/54‑1 (MQ=255)
                                                   cGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGTGCAAGCGTCGCATCAGACGtt       >  1:45135/1‑70 (MQ=255)
                                                   cGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGtt       >  1:249496/1‑70 (MQ=255)
                                                    gCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGa     <  1:3578675/71‑1 (MQ=255)
                                                     ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATc              >  1:3743770/1‑61 (MQ=255)
                                                     ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATc              >  1:1167606/1‑61 (MQ=255)
                                                     ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATc              >  1:995175/1‑61 (MQ=255)
                                                     ccTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAg            <  1:1996159/63‑1 (MQ=255)
                                                      cTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGtt       >  1:1535208/1‑67 (MQ=255)
                                                      cTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGAtt   <  1:461836/71‑1 (MQ=255)
                                                       ttATCCGGCCTACATTTGACAGCCGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg            >  1:1539771/1‑61 (MQ=14)
                                                        tATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGAtt   >  1:2892249/1‑69 (MQ=255)
                                                         aTCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTg  <  1:816944/69‑1 (MQ=255)
                                                             aGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcg                             <  1:2095870/38‑1 (MQ=255)
                                                             aGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcg                             <  1:937388/38‑1 (MQ=255)
                                                             aGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGAcg                             <  1:1018290/38‑1 (MQ=255)
                                                             |                                                                  
GCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTG  >  minE/1985373‑1985498

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