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 | 93.9% | 100.0 / 1.6 | 33 | W30C (TGG→TGT) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base T (13/18); total (14/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01 |
AGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49829‑49996 | aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTaaa > 1:156840/1‑90 (MQ=255) gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAc > 2:566524/1‑90 (MQ=255) gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAAc > 2:379798/1‑90 (MQ=255) gATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAACCGCCATGCCTTGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCaa < 2:477529/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac < 2:486706/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac < 2:248643/90‑1 (MQ=255) aTTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAac < 1:732008/90‑1 (MQ=255) cggcgTGAGCGGTAGATCGCTTTATCGGCATGGAAACCCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCt < 2:146296/90‑1 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 2:620321/1‑90 (MQ=255) ggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCtt > 1:124725/1‑90 (MQ=255) ttAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCTCCTGTTTTAAAGACAAAACCTtaaat > 1:544117/1‑90 (MQ=255) ttAGCGGTAGATCGCGTTATCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTtaaat < 2:568918/90‑1 (MQ=255) aGATCGGGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt < 1:613999/90‑1 (MQ=255) aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt < 1:292839/90‑1 (MQ=255) aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt < 2:307986/90‑1 (MQ=255) aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGt < 1:395535/90‑1 (MQ=255) gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTa > 2:466369/1‑90 (MQ=255) tATCGGCATGGAAAACGCCATGCCGGGGACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg < 1:453288/90‑1 (MQ=255) cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAccc < 1:484552/76‑1 (MQ=255) cGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAAccc > 2:484552/1‑76 (MQ=255) gCATGAAAACCCCCATGCCGGGGACCCTCCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc < 2:424616/90‑1 (MQ=255) aTGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCAt > 2:75075/1‑90 (MQ=255) tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg > 2:3213/1‑74 (MQ=255) tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTg < 1:3213/74‑1 (MQ=255) aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc > 2:511237/1‑53 (MQ=255) aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc < 1:511237/53‑1 (MQ=255) aTGCCGTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca > 2:267711/1‑90 (MQ=255) gTGGAACCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCgg > 1:652895/1‑90 (MQ=255) ggACCCTGCCTGCCGATCTCGCCTGTTTTAAACGCAACACCTTAAATAACCCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtc < 1:264181/90‑1 (MQ=255) gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg < 2:587065/52‑1 (MQ=255) gccGATCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATgg > 1:587065/1‑52 (MQ=255) aTCTCGCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGAc < 2:143338/90‑1 (MQ=255) gCCTGTTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa < 1:75075/90‑1 (MQ=255) | AGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAA > NZ_CP009273/49829‑49996 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |