Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,363,413 | 0 | G | A | 63.0% | 23.5 / 24.5 | 27 | W632* (TGG→TGA) | arnA | bifunctional UDP‑4‑amino‑4‑deoxy‑L‑arabinose formyltransferase/UDP‑glucuronic acid oxidase ArnA |
Reads supporting (aligned to +/- strand): ref base G (4/6); new base A (5/12); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCCGCTGCGCCATCATTTCCCACCGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATA > NZ_CP009273/2363280‑2363530 | aTCCGCTGCGCCATCATTTCCCACCGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGcc < 2:549361/139‑1 (MQ=255) tcatTTCCCACCGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATg < 2:475006/139‑1 (MQ=255) ccACCGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGaa < 1:35810/138‑1 (MQ=255) ccGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCa > 1:462762/1‑139 (MQ=255) tGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATcgac > 2:784227/1‑139 (MQ=255) tCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTg > 2:773709/1‑139 (MQ=255) tCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGc < 2:462762/137‑1 (MQ=255) gTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGAtt < 1:464463/139‑1 (MQ=255) gTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATttc < 1:502176/138‑1 (MQ=255) gTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTgcg < 1:133186/139‑1 (MQ=255) gTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTgcg < 1:784227/139‑1 (MQ=255) cagcTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCAcc > 1:283914/1‑139 (MQ=255) cagcTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTgcg > 1:737164/1‑135 (MQ=255) cagcTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTgcg < 2:737164/135‑1 (MQ=255) gcTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGt < 2:283914/139‑1 (MQ=255) tacGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATc > 2:701757/1‑139 (MQ=255) gCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGt < 2:292355/130‑1 (MQ=255) gCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGt > 1:292355/1‑130 (MQ=255) gTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGAc < 2:446751/139‑1 (MQ=255) aaCCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTa < 1:701757/139‑1 (MQ=255) aaCCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTa < 1:600281/139‑1 (MQ=255) gAGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCa < 2:385715/139‑1 (MQ=255) gAGCAGCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTGCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCa < 2:313860/139‑1 (MQ=255) aGCATCCGCAATGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCAt > 1:732582/1‑139 (MQ=255) tGCCCACCGCTGCCTGGACTGAGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATa < 1:306115/139‑1 (MQ=255) gACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGa > 2:409340/1‑93 (MQ=255) gACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGa < 1:409340/93‑1 (MQ=255) | ATCCGCTGCGCCATCATTTCCCACCGTTTGCGGGCTTTCGCGTTGTCGAAAGTAGCAGCTACTACGGCAAAGGATATCAGGACGTAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATA > NZ_CP009273/2363280‑2363530 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |