Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,727,940 T→C 60.5% F59F (TTT→TTC kefB → potassium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,727,9400TC60.5% 29.2 / 45.8 43F59F (TTT→TTCkefBpotassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base C (10/16);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACTTTG  >  minE/2727872‑2728005
                                                                    |                                                                 
cTGTCAGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACtttt                                                                  <  1:389767/70‑1 (MQ=255)
cTGGCAGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTc                                                                 <  1:748244/71‑1 (MQ=255)
   gCAGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCgg                                                               <  1:497290/70‑1 (MQ=255)
     aGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAAc                                                            <  1:202506/71‑1 (MQ=255)
     aGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAAc                                                            <  1:3426712/71‑1 (MQ=255)
     aGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAAc                                                            >  1:2151201/1‑71 (MQ=255)
      gggATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACt                                                           <  1:2973771/71‑1 (MQ=255)
       ggATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTc                                                          >  1:1236214/1‑71 (MQ=255)
           tGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCg                                                      >  1:645429/1‑71 (MQ=255)
            gCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCg                                                      <  1:3365483/70‑1 (MQ=255)
            gCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCg                                                      <  1:1364845/70‑1 (MQ=255)
                 tGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTAt                                                >  1:135535/1‑71 (MQ=255)
                   gCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTc                                              >  1:3681987/1‑71 (MQ=255)
                         gggggCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATg                                        >  1:925179/1‑71 (MQ=255)
                         gggggCTAGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATg                                        >  1:349930/1‑71 (MQ=255)
                          ggggCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGt                                       <  1:102285/71‑1 (MQ=255)
                            ggCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCgg                                                               <  1:836928/45‑1 (MQ=255)
                               tGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCgg                                                        <  1:3264345/49‑1 (MQ=255)
                               tGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCgg                                                        <  1:2705634/49‑1 (MQ=255)
                               tGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCgg                                                        <  1:1216483/49‑1 (MQ=255)
                                  ggTTTATTAGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTatcatc                               <  1:2586353/71‑1 (MQ=255)
                                    tttATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCgg                                                        <  1:1052473/44‑1 (MQ=255)
                                    tttATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCgg                                                        <  1:3225462/44‑1 (MQ=255)
                                      tattaGCGACGTCGATGAGATCCTCCACTTTTCGGAAc                                                            >  1:2668325/1‑38 (MQ=255)
                                      tattaGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGcc                           >  1:904082/1‑71 (MQ=255)
                                        ttaGCGACGTCGATGAGATCCTCCACTTCTAGCAACTCGGCGTGGTATTCCTg                                           >  1:2341479/1‑53 (MQ=255)
                                         taGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCtt                         >  1:2331147/1‑70 (MQ=255)
                                         taGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCtt                         >  1:3015455/1‑70 (MQ=255)
                                         taGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCg                                                      <  1:2535145/41‑1 (MQ=255)
                                         taGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTgg                                                   >  1:2371195/1‑44 (MQ=255)
                                         taGCGACGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTgg                                                   >  1:1837547/1‑44 (MQ=255)
                                              aCGTCGATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTg                   <  1:1786630/71‑1 (MQ=255)
                                                gTCGATGAGATCCTGCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGaa                 <  1:2905342/71‑1 (MQ=255)
                                                gTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGa                  >  1:3593491/1‑70 (MQ=255)
                                                 tCGATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCtt                         >  1:3212610/1‑62 (MQ=255)
                                                   gATGAGATCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATc               <  1:1910112/70‑1 (MQ=255)
                                                       agaTCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTc           >  1:2972636/1‑70 (MQ=255)
                                                        gaTCCTCCACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCa         <  1:2927323/71‑1 (MQ=255)
                                                           cctccACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAAc      >  1:1772147/1‑71 (MQ=255)
                                                               cACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACtt    <  1:986356/69‑1 (MQ=255)
                                                               cACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACtt    <  1:2306706/69‑1 (MQ=255)
                                                               cACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACtt    <  1:1180362/69‑1 (MQ=255)
                                                               cACTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCATCCAAACttt   <  1:2271121/70‑1 (MQ=255)
                                                                 cTTCTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACTTTg  >  1:2679769/1‑69 (MQ=255)
                                                                   tttCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCaaa       >  1:2216570/1‑62 (MQ=255)
                                                                   tttCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCaaa       >  1:2199057/1‑62 (MQ=255)
                                                                   tttCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCaaa       >  1:944740/1‑62 (MQ=255)
                                                                    |                                                                 
CTGGCAGGGATTGCAATTGGCCCGTGGGGGCTGGGGTTTATTAGCGACGTCGATGAGATCCTCCACTTTTCGGAACTCGGCGTGGTATTCCTGATGTTTATCATCGGCCTTGAGTTGAATCCCTCCAAACTTTG  >  minE/2727872‑2728005

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: