Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 49,912 | G→T | W30C (TGG→TGT) | folA → | type 3 dihydrofolate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 49,912 | 0 | G | T | 96.4% | 183.6 / ‑0.8 | 55 | W30C (TGG→TGT) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base G (2/0); new base T (27/26); total (29/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01 |
TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49827‑49996 | tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa < 2:977996/90‑1 (MQ=255) tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa > 1:911515/1‑90 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa > 2:701688/1‑90 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCGCCGATCTCGCCTGTTTTaa > 2:1173843/1‑90 (MQ=255) aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 1:1311938/1‑90 (MQ=255) aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 1:877904/1‑90 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 2:911515/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 1:299835/90‑1 (MQ=255) tGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGGGGAACCTGCCTGCCGGTCTCGCCTGTTTTAAACGCAacac > 1:844023/1‑90 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:1508269/90‑1 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:446392/90‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 1:1092251/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 2:753906/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 2:1132412/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 1:1320507/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCAGTGGGACCTTCCTGCCGATCTCGTCTGGTTTAAACGCAACACCtt > 2:1423020/1‑90 (MQ=255) aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt < 1:254681/90‑1 (MQ=255) gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg > 1:37640/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:1084272/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:123694/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:1483363/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCCGCCTGGCGATCTCGCCCGTTTTAAACGCAAACCCTTAAATAAACCCGTGATTa > 1:1555265/1‑90 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg < 2:576227/90‑1 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAACCGTGATTATgg > 2:666312/1‑90 (MQ=255) tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 2:533413/90‑1 (MQ=255) tATCGGCATGAAAACCCCCATCCGGGGGAACCCCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 1:986972/90‑1 (MQ=255) tgggCATGAAAACCCCCCTGCGGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcc < 1:1173843/88‑1 (MQ=255) cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg > 1:1239813/1‑90 (MQ=255) aggaaaaaCGCCTTGCCGGGAAACCGGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 1:1344316/86‑1 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt > 2:588766/1‑61 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 2:675854/90‑1 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt > 1:968360/1‑90 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 2:550933/90‑1 (MQ=255) aTGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:588766/61‑1 (MQ=255) tGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa < 1:587977/90‑1 (MQ=255) aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt < 2:1320507/90‑1 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca > 2:659921/1‑90 (MQ=255) tGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa > 2:23080/1‑90 (MQ=255) gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg > 2:434503/1‑90 (MQ=255) ggACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACCCCTTAAATAACCCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 1:197816/90‑1 (MQ=255) ggAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 2:713312/90‑1 (MQ=255) ggAACCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 1:378179/90‑1 (MQ=255) gaccctccctgCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcg < 2:1081447/87‑1 (MQ=255) ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg < 2:1239813/90‑1 (MQ=255) ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg > 1:442184/1‑90 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc < 2:1173665/63‑1 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc > 1:1173665/1‑63 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt > 2:1037403/1‑90 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt > 1:1515171/1‑90 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt > 2:907079/1‑90 (MQ=255) gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa > 2:627518/1‑90 (MQ=255) gCCTGTTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 2:904921/90‑1 (MQ=255) gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 1:161925/90‑1 (MQ=255) gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 2:1518945/90‑1 (MQ=255) gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTAGGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 1:884225/90‑1 (MQ=255) | TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49827‑49996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |