Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 49,910 | T→G | W30G (TGG→GGG) | folA → | type 3 dihydrofolate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 49,910 | 0 | T | G | 96.2% | 160.4 / ‑1.9 | 52 | W30G (TGG→GGG) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (20/30); total (20/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01 |
GATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49825‑49996 | gATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATAGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGtt < 2:511731/90‑1 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTaa > 2:968890/1‑90 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGGTTaa > 2:548706/1‑90 (MQ=255) gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAAc > 2:710007/1‑90 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa < 1:108514/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa < 1:218462/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa < 1:631320/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa < 1:877082/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGAAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCaa < 1:464139/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac < 2:52582/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac < 2:397805/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac < 2:297961/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac < 1:777775/90‑1 (MQ=255) gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc > 1:699087/1‑90 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCt > 1:146660/1‑90 (MQ=255) cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCt < 2:724350/90‑1 (MQ=255) ggcgTTAGCGGTGGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt > 1:612585/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac > 1:548652/1‑85 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAac < 2:548652/85‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt > 2:381839/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt > 2:771039/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt > 1:215832/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCtt > 2:509476/1‑90 (MQ=255) ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaat < 2:839534/90‑1 (MQ=255) ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaat > 2:312763/1‑90 (MQ=255) tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaata < 1:777832/90‑1 (MQ=255) aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGt < 2:639490/90‑1 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCGCCGATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:668038/1‑90 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc < 1:243368/68‑1 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACAcc > 2:243368/1‑68 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg > 2:133621/1‑90 (MQ=255) ttATCGGAATGAAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg < 1:392859/90‑1 (MQ=255) tATCGGCATGAAAACCGCCATGCCGTGGACCCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 2:651042/90‑1 (MQ=255) aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 2:668038/90‑1 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 2:273945/90‑1 (MQ=255) gAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAcc < 2:359473/90‑1 (MQ=255) aaaCGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg < 2:682015/90‑1 (MQ=255) cGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGa > 2:153902/1‑90 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa < 2:15445/57‑1 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa > 1:15445/1‑57 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca > 1:365652/1‑90 (MQ=255) ccGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaatc < 2:253412/90‑1 (MQ=255) ggAACCTCCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 1:787529/90‑1 (MQ=255) ggAACCTCCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 2:918059/90‑1 (MQ=255) ctgcctgcCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCg > 2:306722/1‑90 (MQ=255) tgcctgccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGt > 1:781703/1‑90 (MQ=255) aTCTCGCCTGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc < 1:949611/90‑1 (MQ=255) aTCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc < 2:16593/90‑1 (MQ=255) aTCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc < 1:92676/90‑1 (MQ=255) tcGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCa < 1:710007/90‑1 (MQ=255) gCCGGGTTTAACCGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 1:463530/90‑1 (MQ=255) gCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 1:147111/90‑1 (MQ=255) | GATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49825‑49996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |