Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,562,589 A→G 9.5% G496G (GGT→GGC lldP ← L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,562,5890AG 9.5% 230.4 / 12.5 84G496G (GGT→GGClldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base A (43/33);  new base G (6/2);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

AATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCG  >  minE/2562519‑2562658
                                                                      |                                                                     
aaTCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                                                       <  1:389391/71‑1 (MQ=255)
  tCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTAcc                                                                     <  1:3769706/71‑1 (MQ=255)
    aGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCg                                                                    <  1:1210235/70‑1 (MQ=255)
     gACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGt                                                                  <  1:232027/71‑1 (MQ=255)
      aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                                                       >  1:3142298/1‑65 (MQ=255)
      aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTa                                                                       >  1:2317194/1‑65 (MQ=255)
      aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                                                                 <  1:2444038/71‑1 (MQ=255)
      aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                                                                 <  1:2736074/71‑1 (MQ=255)
      aCTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTg                                                                 <  1:2459157/71‑1 (MQ=255)
          tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGc                                                             >  1:677748/1‑71 (MQ=255)
          tttGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGc                                                             >  1:2717157/1‑71 (MQ=255)
          tcggCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGc                                                             >  1:421311/4‑71 (MQ=255)
            tGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCa                                                           >  1:1028580/1‑71 (MQ=255)
            ggCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGccg                                                           >  1:3720287/2‑70 (MQ=255)
             gCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGtta                                                                >  1:1182852/1‑63 (MQ=255)
             gCCCACCATGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGcc                                                            <  1:679284/69‑1 (MQ=255)
              cccACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCAcc                                                         <  1:2061269/71‑1 (MQ=255)
               ccaccaTGCCCGTCGCGGCGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGCCAGTTACGccgccg                                                        >  1:159661/1‑71 (MQ=255)
               ccaccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCg                                                        <  1:2547204/71‑1 (MQ=255)
               ccaccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCg                                                        <  1:1039501/71‑1 (MQ=255)
                caccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCgg                                                       <  1:659663/71‑1 (MQ=255)
                caccaGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCgg                                                       <  1:732329/71‑1 (MQ=255)
                   caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                                                    <  1:2312697/71‑1 (MQ=255)
                   caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                                                    >  1:3293905/1‑71 (MQ=255)
                   caGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTCACCGGTGGCGTCACCggtaa                                                    <  1:2461919/71‑3 (MQ=38)
                          aCCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggc                                             >  1:1933971/1‑71 (MQ=255)
                           ccgccgCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcg                                            >  1:2456872/1‑71 (MQ=255)
                           ccgccgCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcg                                            >  1:464145/1‑71 (MQ=255)
                            cgcggcGCAGGCCACGGCGATAGATTGCGGGGAGATCATCTTGccag                                                                   <  1:2784748/47‑3 (MQ=16)
                             gccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                  >  1:1912159/1‑63 (MQ=255)
                             gccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcgg                                           >  1:3036914/1‑70 (MQ=255)
                              ccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGt                                                                  >  1:1236178/1‑46 (MQ=255)
                              ccgcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcggc                                          >  1:1666617/1‑70 (MQ=255)
                                gcACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:3799266/1‑71 (MQ=255)
                                 cacaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:3749515/1‑70 (MQ=255)
                                 cacaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:2603306/1‑70 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCggtgg                                                    <  1:3764947/56‑1 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:3236348/1‑69 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:2817278/1‑69 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:2224424/1‑69 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAAc                                       >  1:2893795/1‑69 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa                                     <  1:2730025/71‑1 (MQ=255)
                                  acaGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa                                     <  1:2734771/71‑1 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGAc                                                               >  1:3203072/1‑43 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGAc                                                               >  1:787385/1‑43 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                  >  1:27069/1‑56 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                  >  1:620107/1‑56 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAg                                    >  1:3752057/1‑70 (MQ=255)
                                    aGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                   <  1:123004/71‑1 (MQ=255)
                                     ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAg                                    >  1:1264295/1‑69 (MQ=255)
                                     ggCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGc                                   >  1:1758728/1‑70 (MQ=255)
                                      gcgatagcgataGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                 <  1:2127498/71‑1 (MQ=255)
                                       cgatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                 >  1:3345587/1‑70 (MQ=255)
                                       cgatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCaa                                 >  1:3319831/1‑70 (MQ=255)
                                        gatagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggcggc                                          <  1:2875224/60‑1 (MQ=255)
                                          tagcgatagATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGa                             <  1:3778737/71‑1 (MQ=255)
                                               atagATTGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTTGCTGCCACCAGCagg                                >  1:1775033/1‑61 (MQ=255)
                                                 agATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTa                                                  >  1:673575/1‑43 (MQ=255)
                                                  gATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                      <  1:2433737/70‑1 (MQ=255)
                                                  gATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGAc                     <  1:1071296/71‑1 (MQ=255)
                                                   aTTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGAc                     >  1:3532808/1‑70 (MQ=255)
                                                    ttGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCa                          >  1:2375810/1‑64 (MQ=255)
                                                     tGCGGGGAGATCATCTTGCCAGTTACGCCGCCGCTGGTGTTtgct                                            >  1:2848776/1‑44 (MQ=12)
                                                      gCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGcc                  >  1:1731520/1‑70 (MQ=255)
                                                       cGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTggc                                             <  1:327789/42‑1 (MQ=255)
                                                         gtgagaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCaga                        <  1:495561/59‑1 (MQ=255)
                                                          ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCa                                     <  1:256257/47‑1 (MQ=255)
                                                          ggAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAAttt             >  1:640353/1‑71 (MQ=255)
                                                          ggAGATCATCTTACCGGTGACGCCACAGGTGGTATTGGCGGCa                                         <  1:1142865/43‑1 (MQ=255)
                                                           gagaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg            <  1:1245491/71‑1 (MQ=255)
                                                            agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCaga                        <  1:459161/58‑1 (MQ=255)
                                                            agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAgaga                      >  1:2529014/1‑60 (MQ=255)
                                                            agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgt           >  1:826490/1‑71 (MQ=255)
                                                            agaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgt           <  1:736675/71‑1 (MQ=255)
                                                             gaTCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgtt          >  1:955294/1‑71 (MQ=255)
                                                              atcatcTTACCGGTGACGCCGCCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg         <  1:2834853/71‑1 (MQ=255)
                                                              atcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg         >  1:308414/1‑71 (MQ=255)
                                                              atcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg         >  1:1820471/1‑71 (MQ=255)
                                                               tcatcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg         <  1:2209667/70‑1 (MQ=255)
                                                                catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGa                             >  1:3283652/1‑49 (MQ=255)
                                                                catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttgttg         >  1:3396376/1‑69 (MQ=255)
                                                                catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTtgtg       >  1:515946/1‑71 (MQ=255)
                                                                catcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTtgtg       >  1:2872222/1‑71 (MQ=255)
                                                                 atcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTtgtg       <  1:43085/70‑1 (MQ=255)
                                                                  tcTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGAc                     >  1:2254794/1‑55 (MQ=255)
                                                                     tACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATttg            <  1:546670/61‑1 (MQ=255)
                                                                     tACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCg  >  1:377805/1‑71 (MQ=255)
                                                                      aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCg  <  1:34961/70‑1 (MQ=255)
                                                                      aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCg  <  1:2400029/70‑1 (MQ=255)
                                                                      aCCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCg  <  1:1388077/70‑1 (MQ=255)
                                                                      |                                                                     
AATCAGACTCTTTGCCCACCAGGCCTACCGCCGCACAGGCGATAGCGATAGATTGCGGGGAGATCATCTTACCGGTGACGCCACCGGTGGTATTGGCGGCAACCAGCAACAGATCAGAGACGCCAATTTGTTGTGCTGCG  >  minE/2562519‑2562658

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: