Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,088,442 A→T 7.0% G110G (GGA→GGT yhaI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,088,4420AT 7.0% 202.9 / 6.4 71G110G (GGA→GGTyhaIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base T (3/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAA  >  minE/2088373‑2088510
                                                                     |                                                                    
gTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGa                                                                      >  1:1751347/1‑70 (MQ=255)
 tttATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                    <  1:2800223/71‑1 (MQ=255)
   tATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAAc                                                                   >  1:656025/1‑70 (MQ=255)
    aTTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACg                                                                  <  1:1183070/70‑1 (MQ=255)
    aTTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACg                                                                  <  1:908424/70‑1 (MQ=255)
     ttGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                                 <  1:1269827/70‑1 (MQ=255)
          tggctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCg                                                             >  1:1513797/1‑69 (MQ=255)
            gctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAg                                                          <  1:2444319/70‑1 (MQ=255)
            gctCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGt                                                         >  1:940797/1‑71 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGttt                                                       <  1:3322472/69‑1 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTgg                                                     >  1:2094172/1‑71 (MQ=255)
                gTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTgg                                                     >  1:2032576/1‑71 (MQ=255)
                     tCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCa                                                 >  1:948818/1‑70 (MQ=255)
                      cTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCaaa                                               >  1:998736/1‑71 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                                 <  1:877461/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                                 <  1:747196/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGa                                                                 <  1:2831140/51‑1 (MQ=255)
                        ggTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAAt                                              >  1:133507/1‑70 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTTGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAAtt                                             <  1:991033/70‑1 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTa                                            <  1:3614027/71‑1 (MQ=255)
                         gTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTa                                            <  1:493242/71‑1 (MQ=255)
                          tttttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTaa                                           <  1:1883243/71‑1 (MQ=255)
                            tttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAAt                                          <  1:1661962/70‑1 (MQ=255)
                            tttttCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATc                                         >  1:3649611/1‑71 (MQ=255)
                                 cTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                     >  1:1537249/1‑70 (MQ=255)
                                    cacaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                                  <  1:2860554/70‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                    >  1:1454653/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                    >  1:2225127/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGaaa                                                                    >  1:2387553/1‑35 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                      <  1:2417027/65‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                      <  1:2694569/65‑1 (MQ=255)
                                     acaGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                                  <  1:1860574/69‑1 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTa                                      <  1:3518529/63‑1 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAt                                  >  1:181967/1‑67 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTg                                >  1:285557/1‑69 (MQ=255)
                                       aGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTAGGTTCAAATTAATCTTAGAATTg                                >  1:382726/1‑69 (MQ=255)
                                        gAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTggg                              >  1:943477/1‑70 (MQ=255)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                                >  1:3542184/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                                >  1:3636099/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTTCCGTAACGAt                                                                >  1:554743/1‑35 (MQ=21)
                                         aaGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGc                            >  1:3141962/1‑71 (MQ=255)
                                            gTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGc                            <  1:1259329/68‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCg                           <  1:682525/66‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCg                           <  1:3157854/66‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATAtt                      <  1:1368469/71‑1 (MQ=255)
                                               cTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATAtt                      <  1:2993969/71‑1 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTg                                                      >  1:2793462/1‑36 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:3538681/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:2721553/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:1624977/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:2808577/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:2703910/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGat                         >  1:2306034/1‑65 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTc                    <  1:1741141/70‑1 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCg                   >  1:2929257/1‑71 (MQ=255)
                                                  cTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCg                   >  1:695694/1‑71 (MQ=255)
                                                       agcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt              <  1:2947200/71‑1 (MQ=255)
                                                       agcaATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt              <  1:1493511/71‑1 (MQ=255)
                                                        cggaTCGTTTCGGTAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt              <  1:2905270/67‑1 (MQ=255)
                                                          gttCGTTTCCGTAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCt              <  1:1861340/66‑1 (MQ=255)
                                                          aaTCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt           <  1:2320696/71‑1 (MQ=255)
                                                          aaTCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt           >  1:207655/1‑71 (MQ=255)
                                                             cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                     <  1:2015858/42‑1 (MQ=255)
                                                             cGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAg                                     <  1:2011637/42‑1 (MQ=255)
                                                              gTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTttatt       <  1:2513754/71‑1 (MQ=255)
                                                               ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAAtt                                 <  1:2604305/44‑1 (MQ=255)
                                                               ttACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTttatt       >  1:1332242/1‑70 (MQ=255)
                                                                 aCGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTaa     >  1:710161/1‑70 (MQ=255)
                                                                  cGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt           >  1:2220667/1‑63 (MQ=255)
                                                                  cGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtttt           >  1:2099133/1‑63 (MQ=255)
                                                                   ggAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGa                                    >  1:3000988/1‑37 (MQ=255)
                                                                   ggAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGa                                    >  1:1217034/1‑37 (MQ=255)
                                                                   ggAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGa                                    >  1:905792/1‑37 (MQ=255)
                                                                   ggAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGa                                    >  1:595692/1‑37 (MQ=255)
                                                                   ggAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACaa  <  1:2637351/71‑1 (MQ=255)
                                                                    gAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCtt             <  1:1601267/59‑1 (MQ=255)
                                                                     |                                                                    
GTTTATTGGCTGGCTCGTGCTTCTGGTTTTTTTCTGCACAGAAGGTACTTCTGGCAGCAATCGTTACGGAAACGATCCGAAGTTTGGTTCAAATTAATCTTAGAATTGGGGCGATATTTCGCCCCTTTTTATTAACAA  >  minE/2088373‑2088510

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: