Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,826,895 C→G 63.4% R46R (CGC→CGG mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,826,8950CG63.4% 45.2 / 49.4 60R46R (CGC→CGGmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (17/5);  new base G (24/14);  total (41/19)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCGGCGGGAGAGCCG  >  minE/1826833‑1826952
                                                              |                                                         
gagaTCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc                                                      >  1:2963186/1‑68 (MQ=255)
        gctgTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgc                                              >  1:3186286/1‑68 (MQ=255)
        gctgTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgc                                              >  1:2164316/1‑68 (MQ=255)
           gTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTgg                                           >  1:2716797/1‑68 (MQ=255)
           gTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGa                                          >  1:1261863/1‑69 (MQ=255)
           gTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGa                                          <  1:256605/69‑1 (MQ=255)
           gTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGa                                          <  1:2443226/69‑1 (MQ=255)
           gTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGa                                          >  1:1768370/1‑69 (MQ=255)
              ttACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGatat                                       >  1:1369490/1‑69 (MQ=255)
              ttACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGata                                        >  1:1827914/1‑68 (MQ=255)
              ttACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGa                                          <  1:749723/66‑1 (MQ=255)
                        ggTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGAc                              >  1:2950671/1‑68 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:3103677/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:502230/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:2608307/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:497762/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:465829/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:3212183/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:1546343/1‑69 (MQ=255)
                              ttttATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAAcgc                       >  1:1799245/1‑69 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAAcg                        >  1:1590058/1‑59 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:3155659/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:2981389/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:393185/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:371913/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:3295606/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:351999/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:3126249/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:3304938/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:26868/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:2578039/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:597146/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1972785/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:6910/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1793313/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:77683/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1429877/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1407465/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1406901/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1284208/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:1174529/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:571037/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:582960/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:630676/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:790233/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:337197/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:2991024/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:2888486/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:2813674/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:2788582/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:2770572/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               >  1:2242917/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:1952447/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:1885065/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:1645150/68‑1 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACAAAACGCGGTGCTTca              >  1:2463666/1‑68 (MQ=255)
                                       cTGTTTTATGACGACGCAAAACGGGAGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTca              >  1:40441/1‑68 (MQ=255)
                                        tGTTTTATGACGACGCAAAACGGGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTc               <  1:238726/67‑1 (MQ=255)
                                                   gacgCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCGGCGGGAGAGCCg  <  1:1800260/69‑1 (MQ=255)
                                                   gacgCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCGGCGGGAGAGCCg  <  1:1090975/69‑1 (MQ=255)
                                                              |                                                         
GAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTGACCAAACGCGGTGCTTCGGCGGGAGAGCCG  >  minE/1826833‑1826952

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: