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.2% | 172.4 / 0.8 | 53 | W30C (TGG→TGT) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base T (24/27); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01 |
TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49827‑49996 | tCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTa < 1:858736/90‑1 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaa > 1:198029/1‑90 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaa > 1:54211/1‑90 (MQ=255) aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 2:857200/1‑90 (MQ=255) aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 1:463814/1‑90 (MQ=255) aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 2:1227284/1‑90 (MQ=255) tGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGG‑AAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 2:1434178/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 2:1420508/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 1:1429230/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 1:1357064/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCACGTTATCGGGATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 1:330972/90‑1 (MQ=255) ctttGCGCCGTTAGCGGTAGTTCGCGTTATCGGCATGGAAAACGCCATGCCGTGAACCCTTCCTGCCGATCTCGCCTGTTTTAAACGCaa < 1:1462495/88‑1 (MQ=255) gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACAcc < 2:759880/90‑1 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:444801/90‑1 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGAAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:1263890/90‑1 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGAAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 1:1295166/90‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa > 1:779785/1‑84 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCaa < 2:779785/84‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 2:421864/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 1:444231/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 2:116676/89‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt > 1:116676/1‑89 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAATCAACACCtt > 1:145756/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCCTGCCGGGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 1:1434615/1‑90 (MQ=255) tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTtaaata < 2:96700/90‑1 (MQ=255) tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGttt < 2:1010384/60‑1 (MQ=255) tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGttt > 1:1010384/1‑60 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGTCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGGGGTTa > 1:918094/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:1418400/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:488525/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:523881/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGGGATTa > 1:123720/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAAGCAACACCTTAAATAAACCCGTGGTTa > 1:920751/1‑90 (MQ=255) ttATCGGCATGGAAAACGCCTTGCCGTGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg < 1:857200/90‑1 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg < 2:198029/90‑1 (MQ=255) tATCGGCATGGAAATCGCCACGCCGGGGAGCCTGCCTGCGGATCTAGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 2:826904/90‑1 (MQ=255) tATCGGCATGAAAACCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 1:86761/90‑1 (MQ=255) ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc < 2:1419037/90‑1 (MQ=255) ggCATGAAAACCCCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc < 1:187617/90‑1 (MQ=255) aaacCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt < 1:322438/90‑1 (MQ=255) aaaaCGCCATGCGGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt < 1:1241477/90‑1 (MQ=255) aaaaCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt < 2:321853/90‑1 (MQ=255) aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCt > 1:1301868/1‑90 (MQ=255) cGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCTATACCTGGGa < 2:35766/90‑1 (MQ=255) ccATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat < 1:296585/90‑1 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca > 1:1251501/1‑90 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca > 1:1454064/1‑90 (MQ=255) ccGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc < 2:182631/90‑1 (MQ=255) gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg > 2:402859/1‑90 (MQ=255) gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg > 2:529518/1‑90 (MQ=255) tGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt > 2:1448597/1‑90 (MQ=255) ggAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 2:1301868/90‑1 (MQ=255) ctgcctgcCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATGGGTCGTCCg > 2:1313117/1‑90 (MQ=255) tgcctgccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAACCGGTCGGCCGt > 1:1444168/1‑90 (MQ=255) gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 2:1418400/90‑1 (MQ=255) | TCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49827‑49996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |