Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,502,076 | A→G | N347D (AAC→GAC) | rlhA → | 23S rRNA 5‑hydroxycytidine C2501 synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,502,076 | 0 | A | G | 92.3% | 31.1 / ‑3.5 | 13 | N347D (AAC→GAC) | rlhA | 23S rRNA 5‑hydroxycytidine C2501 synthase |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (7/5); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01 |
GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCC > NZ_CP009273/1501989‑1502158 | ggCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGAc < 2:86476/90‑1 (MQ=255) gCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACg > 1:204833/1‑90 (MQ=255) gCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACg > 2:220573/1‑90 (MQ=255) tGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGtt > 1:240851/1‑90 (MQ=255) cGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTc < 2:204833/90‑1 (MQ=255) catcTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAAt > 2:498488/1‑90 (MQ=255) cGAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga < 2:104306/90‑1 (MQ=255) cGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGgaga < 2:394637/90‑1 (MQ=255) gAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGaa < 2:240851/90‑1 (MQ=255) gAGCCACTGGCAAATGGCGATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGaa > 2:385601/1‑90 (MQ=255) gATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTgg > 2:410661/1‑90 (MQ=255) gATGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTgg > 2:516177/1‑90 (MQ=255) tGGCCTGGACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGcc < 2:130303/90‑1 (MQ=255) | GGCCTGCCGGTAGGCGAAGTAGTGAAAGTGGCGAAAGATCATCTCGATGTTGCCGTTACCGAGCCACTGGCAAATGGCGATGGCCTGAACGTGTTGATTAAACGTGAAGTCGTCGGTTTTCGTGCCAATACGGTCGAGAAAACCGGAGAAAATCAGTACCGCGTCTGGCC > NZ_CP009273/1501989‑1502158 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |