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
                                                        |                                                           
gTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTa                                                       <  1:2150403/62‑1 (MQ=255)
gTTGAGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTa                                                       <  1:4541958/62‑1 (MQ=255)
    aGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGt                                          <  1:1450832/71‑1 (MQ=255)
    aGTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGt                                          <  1:1480253/71‑1 (MQ=255)
     gTAAAACAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGt                                          <  1:3865839/70‑1 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGaa                                            >  1:2406184/1‑65 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGaa                                            >  1:4877500/1‑65 (MQ=255)
        aaaCAGTAATGTCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGaa                                            >  1:2937483/1‑65 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:768872/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2579157/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2019254/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2657286/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2838154/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:2885342/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:3448903/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:704415/1‑69 (MQ=255)
                   tCGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTAt                             >  1:1304594/1‑69 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATc                            >  1:16078/1‑69 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:3172287/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:381069/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:4116664/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:417899/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2529401/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2447899/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:862199/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:1385258/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:2107327/1‑71 (MQ=255)
                    cGTGAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAg                          >  1:1799040/1‑71 (MQ=255)
                       gAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCa                        >  1:570982/1‑70 (MQ=255)
                       gAAACTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCa                        >  1:1777568/1‑70 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCa                                 <  1:4462188/57‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAac                   <  1:51384/71‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAac                   <  1:1540039/71‑1 (MQ=255)
                           cTAATAAATGCAAAAATTGTCCCGGCTCCGCGCTACTGTGGGAGAAGTCTTTGCGCa                                 <  1:603226/57‑1 (MQ=255)
                                 aaTGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGc             >  1:1104307/1‑71 (MQ=255)
                                 aaTGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGc             >  1:3679155/1‑71 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:4411379/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:244474/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:735549/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:1705746/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:4504843/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:180055/1‑70 (MQ=255)
                                   tGCAAAAATTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCg            >  1:3524425/1‑70 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGGGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:4580801/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:2098252/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:462883/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:1828330/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:1960641/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        <  1:178074/71‑1 (MQ=255)
                                      aaaaaTTGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATAAGCAAAAACACCAGCGTaaa        <  1:3060205/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1804501/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:4347032/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1924121/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:269029/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1222210/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3310340/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:4590371/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3083567/71‑1 (MQ=255)
                                           ttGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:1967004/71‑1 (MQ=255)
                                           ctGTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGc   <  1:3765359/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaa        >  1:4686083/1‑64 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:1946273/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:847574/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2401301/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2465386/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2512576/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:2869144/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:3278054/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:4482252/70‑1 (MQ=255)
                                             gTCCCGGCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTAAAAATGCa  <  1:4335573/70‑1 (MQ=255)
                                                   gCTGCGCGCTACTGTGGGAGAAGTCTTTGCGCAGTATCAGCAAAAACACCAGCGTaaaaa      >  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: