Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 848,817 Δ1 bp coding (407/1254 nt) kch ← voltage‑gated potassium channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE848,8160G.98.6% 303.2 / ‑1.8 71coding (408/1254 nt)kchvoltage‑gated potassium channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (36/34);  minor base C (0/1);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCA  >  minE/848760‑848875
                                                        |                                                           
gTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTa                                                       <  1:2150403/62‑1 (MQ=255)
gTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTa                                                       <  1:4541958/62‑1 (MQ=255)
    aGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGt                                          <  1:1450832/71‑1 (MQ=255)
    aGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGt                                          <  1:1480253/71‑1 (MQ=255)
     gTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGt                                          <  1:3865839/70‑1 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGaa                                            >  1:2406184/1‑65 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGaa                                            >  1:4877500/1‑65 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGaa                                            >  1:2937483/1‑65 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:768872/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2579157/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2019254/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2657286/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2838154/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2885342/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:3448903/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:704415/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:1304594/1‑69 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATc                            >  1:16078/1‑69 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:3172287/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:381069/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:4116664/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:417899/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2529401/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2447899/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:862199/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:1385258/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2107327/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:1799040/1‑71 (MQ=255)
                       gAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCa                        >  1:570982/1‑70 (MQ=255)
                       gAAACTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCa                        >  1:1777568/1‑70 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCa                                 <  1:4462188/57‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAac                   <  1:51384/71‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAac                   <  1:1540039/71‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCT‑CCGCGCTACTGTGGGAGAAGTCTTTGCGCa                                 <  1:603226/57‑1 (MQ=255)
                                 aaTGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGc             >  1:1104307/1‑71 (MQ=255)
                                 aaTGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGc             >  1:3679155/1‑71 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:4411379/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:244474/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:735549/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:1705746/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:4504843/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:180055/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:3524425/1‑70 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGGGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:4580801/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:2098252/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:462883/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:1828330/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:1960641/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:178074/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATAAGCAAAAACACCAGCGTaaa        <  1:3060205/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1804501/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:4347032/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1924121/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:269029/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1222210/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3310340/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:4590371/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3083567/71‑1 (MQ=255)
                                           ttGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1967004/71‑1 (MQ=255)
                                           ctGTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3765359/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        >  1:4686083/1‑64 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:1946273/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:847574/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2401301/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2465386/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2512576/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2869144/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:3278054/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:4482252/70‑1 (MQ=255)
                                             gTCCCGGCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:4335573/70‑1 (MQ=255)
                                                   gCTGC‑GCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaaaa      >  1:3366026/1‑60 (MQ=255)
                                                        |                                                           
GTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCA  >  minE/848760‑848875

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: ‑