Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,334,062 C→T 51.5% P97P (CCG→CCA wcaI ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,334,0620CT51.5% ‑3.4 / 58.9 66P97P (CCG→CCAwcaIpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base T (22/12);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACAGCGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTT  >  minE/1334001‑1334128
                                                             |                                                                  
aaaaCAGCGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             <  1:1337793/69‑1 (MQ=255)
aaaaCAGCGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             <  1:2741360/69‑1 (MQ=255)
aaaaCAGCGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             <  1:2452683/69‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:3251826/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:1891151/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:574888/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:1965818/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:2096086/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:784418/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:111062/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:330252/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:28547/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:3274524/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:1039907/68‑1 (MQ=255)
 aaaCAGCGTTGGCACCACGCCGATAATGCGATCCGGCTTCAAGCGACGTTGCGCCATCAGTGgaaaga                                                             <  1:1384024/68‑1 (MQ=255)
       cGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:2567972/69‑1 (MQ=255)
       cGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:2163858/69‑1 (MQ=255)
       cGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:2000801/69‑1 (MQ=255)
       cGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:928249/69‑1 (MQ=255)
       cGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:1678115/69‑1 (MQ=255)
          tGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAActgct                                                      <  1:752453/66‑1 (MQ=255)
           ggCGCCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGag                                                   >  1:1008239/1‑67 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:3147159/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:1489018/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:3107435/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:801333/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:2306986/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:1242975/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:1246909/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:60636/1‑68 (MQ=255)
           ggCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAAAACTGCTGACg                                                  >  1:1060575/1‑68 (MQ=255)
               ccacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCaa                                              <  1:3089190/69‑1 (MQ=255)
               ccacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCaa                                              <  1:1087869/69‑1 (MQ=255)
                cacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             >  1:2179204/1‑53 (MQ=255)
                cacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             >  1:405342/1‑53 (MQ=255)
                cacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             >  1:1220885/1‑53 (MQ=255)
                cacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGgaaaga                                                             >  1:2527552/1‑53 (MQ=255)
                cacGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCaaa                                             <  1:2196282/69‑1 (MQ=255)
                        taatGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCa                                     <  1:147874/69‑1 (MQ=255)
                        taatGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCa                                     <  1:800198/69‑1 (MQ=255)
                                 ccGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGAc                                                   >  1:374681/1‑46 (MQ=255)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:3313874/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:523646/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:78559/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:239929/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:3388568/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:3426277/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:3051850/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1013933/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1109605/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1142954/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1189074/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1317333/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1344264/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1352055/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:1692939/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2072716/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2285118/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2423215/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2579928/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2718372/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:2969640/1‑51 (MQ=39)
                                        ccAGCGACGTTGCGCCATCAGTGGAAAGAAACTGCTGACGGCAAAACTGccg                                      >  1:3165349/1‑51 (MQ=39)
                                                           aGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGtttt  <  1:1431587/69‑1 (MQ=255)
                                                           aGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGtttt  <  1:794779/69‑1 (MQ=255)
                                                            gCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGtttt  <  1:1085984/68‑1 (MQ=255)
                                                             |                                                                  
AAAACAGCGTTGGCACCACGCCAATAATGCGATCCGGCTTCCAGCGACGTTGCGCCATCAGCGGAAAGAAACTGCTGACGGCAAAACTGCCCAGATGCAACAGGCGTTTCAGGGTGCTCGGCTGTTTT  >  minE/1334001‑1334128

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: