Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 636,843 G→A 37.9% V1273V (GTG→GTA mukB → fused chromosome partitioning proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE636,8430GA37.9% 25.7 / 29.0 50V1273V (GTG→GTAmukBfused chromosome partitioning proteins
Reads supporting (aligned to +/- strand):  ref base G (25/6);  new base A (17/2);  total (42/8)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAC  >  minE/636776‑636901
                                                                   |                                                          
aGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAgcgtgc                                                           <  1:1392818/69‑1 (MQ=255)
aGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAgcgtgc                                                           <  1:3335980/69‑1 (MQ=255)
                      cAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:65708/69‑1 (MQ=255)
                      cAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:178423/69‑1 (MQ=255)
                            gtggTTGCAGAACGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAacgca                               <  1:2324041/67‑1 (MQ=255)
                            ggggTTGCAGAACGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAacgca                               <  1:209567/69‑1 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGcc                           >  1:714900/1‑68 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGcc                           >  1:2376107/1‑68 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1995542/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:2602740/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:323119/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1828048/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1706891/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1604302/1‑69 (MQ=255)
                                 tGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCa                          >  1:1000140/1‑69 (MQ=255)
                                  gCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCAt                         <  1:3433839/69‑1 (MQ=255)
                                  gCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCAt                         <  1:3373915/69‑1 (MQ=255)
                                   cAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGa                                     <  1:491180/56‑1 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGTCATGCtg                     >  1:665822/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:444230/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3327198/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:30678/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3344189/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:3386944/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:374859/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:28173/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:921127/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:971496/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2840384/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:1429371/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2753598/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:152387/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2536804/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCtg                     >  1:2280454/1‑68 (MQ=255)
                                      aaCGTATCGTTTGGCCAGGTGAACAGCGTACGTCTGAACGTGAACGTACGTGAAACGCACGCCATGCt                      >  1:2598300/1‑68 (MQ=255)
                                         gTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTg                  <  1:2710583/69‑1 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACTTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2034366/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2532771/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:771969/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:1731215/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:507706/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2178640/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:400830/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2912774/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:3427327/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:2701611/1‑69 (MQ=255)
                                            tCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGAt               >  1:3287140/1‑69 (MQ=255)
                                                          gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:2566288/1‑68 (MQ=255)
                                                          gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:1696169/1‑68 (MQ=255)
                                                          gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:325110/1‑68 (MQ=255)
                                                          gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:3101440/1‑68 (MQ=255)
                                                          gAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAc  >  1:1013241/1‑68 (MQ=255)
                                                                   |                                                          
AGCAGAACCGTATCCGTATGCTCAACCAGGGGTTGCAGAACGTATCGTTTGGTCAGGTGAACAGCGTGCGTCTCAACGTGAACGTGCGTGAAACGCACGCCATGCTACTGGATGTGCTCTCTGAAC  >  minE/636776‑636901

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: