Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 744,026 G→C 38.6% G226A (GGT→GCT)  nagZ → beta N‑acetyl‑glucosaminidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE744,0260GC38.6% 38.0 / 48.7 57G226A (GGT→GCT) nagZbeta N‑acetyl‑glucosaminidase
Reads supporting (aligned to +/- strand):  ref base G (29/6);  new base C (9/13);  total (38/19)
Fisher's exact test for biased strand distribution p-value = 1.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGA  >  minE/743958‑744074
                                                                    |                                                
aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:1052595/69‑1 (MQ=255)
aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:2038634/69‑1 (MQ=255)
aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:581824/69‑1 (MQ=255)
aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:2555192/69‑1 (MQ=255)
aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:498184/69‑1 (MQ=255)
  aTAAACTCGACTCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:2657376/67‑1 (MQ=255)
       cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg                                                  <  1:143267/62‑1 (MQ=255)
       cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTccc                                           <  1:710391/69‑1 (MQ=255)
       cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTccc                                           <  1:2687083/69‑1 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:1537655/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:2110800/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:600667/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:1962219/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:1937982/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:1002015/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:446154/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:2774724/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:2875372/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:3405107/1‑69 (MQ=255)
                   aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa                               >  1:758119/1‑69 (MQ=255)
                                    cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg               <  1:871351/68‑1 (MQ=255)
                                    cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg               <  1:2839318/68‑1 (MQ=255)
                                    cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg               <  1:3286625/68‑1 (MQ=255)
                                    cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg               <  1:1479140/68‑1 (MQ=255)
                                       cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt           >  1:2300647/1‑69 (MQ=255)
                                       cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt           >  1:1319731/1‑69 (MQ=255)
                                       cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt           >  1:3357488/1‑69 (MQ=255)
                                       cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt           >  1:1759735/1‑69 (MQ=255)
                                       cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt           >  1:1516556/1‑69 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:656377/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:3399937/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:3027927/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:1421375/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:1487598/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:1621789/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:670702/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:637812/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:201991/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:203900/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:2341893/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:2083788/65‑1 (MQ=255)
                                         gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTAATGGCTGAAAGCCGTTTTGCGTCAGGaa             <  1:62069/65‑1 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:57566/1‑69 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:3311075/1‑69 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:2604424/1‑69 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:2465361/1‑69 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:1040451/1‑69 (MQ=255)
                                            aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTAATGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt      >  1:1646673/1‑69 (MQ=255)
                                               ttGATGCGCGTCCGGCGAGCGGTTCTCACTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:2590086/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:3387991/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:306386/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:855217/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:1511340/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:3257963/1‑69 (MQ=255)
                                               ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg   >  1:2150274/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:313918/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:1110897/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:1638043/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:1751925/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:235777/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:2493275/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:2663203/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa  >  1:3027576/1‑69 (MQ=255)
                                                tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGGTTTGa  >  1:3112835/1‑69 (MQ=255)
                                                                    |                                                
AAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGA  >  minE/743958‑744074

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: