Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 678,913 | A→G | intergenic (‑64/+20) | hscC ← / ← rihA | molecular chaperone HscC/pyrimidine‑specific ribonucleoside hydrolase RihA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 678,913 | 0 | A | G | 93.8% | 95.7 / ‑0.8 | 32 | intergenic (‑64/+20) | hscC/rihA | molecular chaperone HscC/pyrimidine‑specific ribonucleoside hydrolase RihA |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (19/11); total (20/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01 |
AGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAAC > NZ_CP009273/678833‑679043 | aGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATcaa > 2:73173/1‑138 (MQ=255) tCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACg < 2:763669/84‑1 (MQ=255) tCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACg > 1:763669/1‑84 (MQ=255) tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCTGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:224962/1‑139 (MQ=255) tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:256259/1‑139 (MQ=255) tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 2:18819/1‑139 (MQ=255) tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 1:829883/1‑139 (MQ=255) tGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCaa > 1:633160/1‑139 (MQ=255) cATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAATCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa < 2:847354/139‑1 (MQ=255) cATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa < 2:98116/139‑1 (MQ=255) ttATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACa > 2:602755/1‑136 (MQ=255) tGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTAc > 1:90932/1‑139 (MQ=255) gaaagaaaAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTa < 1:253780/137‑1 (MQ=255) gaaagaaaAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTa < 2:633160/137‑1 (MQ=255) aaGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCAtt < 1:73173/139‑1 (MQ=255) aaaCAGATTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTa > 2:285310/1‑139 (MQ=255) aTTAAGACCACCCAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCt > 1:299049/1‑139 (MQ=255) cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:761864/1‑138 (MQ=255) cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:346311/1‑138 (MQ=255) cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 2:357924/1‑138 (MQ=255) cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:27311/1‑138 (MQ=255) cccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataa > 1:632998/1‑138 (MQ=255) ccAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataata < 2:619723/139‑1 (MQ=255) cAGGTAACGCGTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGataataa < 1:285310/139‑1 (MQ=255) cgcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGAc > 2:817291/1‑92 (MQ=255) cgcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGAc < 1:817291/92‑1 (MQ=255) gcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac < 1:41423/139‑1 (MQ=255) gcgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac < 2:346311/139‑1 (MQ=255) cgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac > 2:529070/1‑138 (MQ=255) cgTGGTGGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATcaac > 2:163516/1‑138 (MQ=255) gtgGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACTATCAAc < 1:436962/80‑1 (MQ=255) gtgGTCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAAc > 2:436962/1‑80 (MQ=255) | AGTTCTGCATTATCCATTGCATTTATATCTTGAAAGAAAAAGAGATAAACAGATTAAGACCACCCAGGTAACGCGTGGTGATCTTAATCAATGACGTGTGTTAAGCGTAAAATTTCAGACGATCAGCCAGTAAATCAACAAAGCCCTGACGATCAACATCGACCATTACGGTGGCATTCGGTTTATTGCCTGTCAGATAATAATAATCAAC > NZ_CP009273/678833‑679043 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |