Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 49,910 | T→A | W30R (TGG→AGG) | 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 | A | 97.2% | 121.5 / NA | 37 | W30R (TGG→AGG) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base T (1/0); major base A (15/20); minor base C (1/0); total (17/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49828‑49996 | cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTaa > 1:218705/1‑90 (MQ=255) cAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTaa > 2:795128/1‑90 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa < 1:827998/90‑1 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa < 2:47646/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac < 1:238508/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac < 1:318354/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac > 1:92375/1‑90 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac < 1:667745/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGCCATGGAAAACACCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAac < 1:572397/90‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt > 2:704592/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCtt > 2:704058/1‑90 (MQ=255) tAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTtaaata < 2:807634/90‑1 (MQ=255) gcgTTATCGGCATGGAAAACGCCCTGCCGGGGAACCTGCCTGCCCAACTCGCCAGGGTTAAAAGCAAAACCTTAAATAAACCCGGGATTa > 2:702727/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 2:872875/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa > 1:123371/1‑67 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa < 2:123371/67‑1 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGGTTAAACGCAACACCTTAAATAAACCCGTGATTa > 1:614736/1‑90 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCGGCCTGCCGAACTCGCCCGGGTTAAACGCAAAACCTTAAATAAACCCCTGATTa > 1:862355/1‑90 (MQ=255) cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTAt > 1:179450/1‑90 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa < 1:5491/64‑1 (MQ=255) ttATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCaa > 2:5491/1‑64 (MQ=255) aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgc < 2:92375/90‑1 (MQ=255) cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa < 2:460768/90‑1 (MQ=255) aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 1:587290/90‑1 (MQ=255) aTGGAAACCGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt < 1:702727/90‑1 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt > 1:479641/1‑90 (MQ=255) tGGAAAACGCCATGCCGTGGACCCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATa < 1:639210/90‑1 (MQ=255) cATGCCGTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatc > 2:553807/1‑90 (MQ=255) gTGGAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg < 1:872875/90‑1 (MQ=255) ggAACCTGCCTGCCGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 1:245875/90‑1 (MQ=255) tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 2:346723/84‑1 (MQ=255) tgcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc > 1:346723/1‑84 (MQ=255) gcctgccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGtt > 1:654516/1‑90 (MQ=255) gccGAACTCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCa > 2:183819/1‑90 (MQ=255) cGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa < 1:493963/90‑1 (MQ=255) gCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 2:196527/90‑1 (MQ=255) gCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 2:393592/90‑1 (MQ=255) | CAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49828‑49996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |