Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 990,954 | 0 | A | G | 61.5% | 0.8 / 11.7 | 13 | F164S (TTC→TCC) | ssuD | FMNH2‑dependent alkanesulfonate monooxygenase |
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (8/0); total (10/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGACGGTTTCTCTCTGCAATAAACGCCGCCAGACC > NZ_CP009273/990867‑991034 | agaTCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTGAAACGCCGGGGag > 1:220279/1‑90 (MQ=255) gaTCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGAAACGGCTGGTGTAACGCCCGGGagg > 2:189676/1‑89 (MQ=255) cAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTGGAACGCCCGGGagagga > 2:91223/1‑88 (MQ=255) gTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTGAATTGCCGGGGAGAGCAGgtttt > 2:470796/1‑89 (MQ=255) gccAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCgcgc < 1:91223/90‑1 (MQ=255) cAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTGAATCGCCGGGGAGAGCAGTTTTTCTCCGCGcac > 1:48571/1‑90 (MQ=255) cTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAat < 2:291291/90‑1 (MQ=255) aCATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTGAATCGCCCGGGAGAGGAGGTTTTCTCCCCGCCCATGAATATGtttt > 2:187281/1‑89 (MQ=255) aTCCGCCAAAGTAAAGTGGCGGATACGGCTGTTTAATCGCCCGGGAGAGCAGTTTTTCTCCGCCGGCATGAATATGTTTTCCGTTGgaaa > 1:505681/1‑89 (MQ=255) aTACGGCTGTTGAATCGCCGGGGAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTTCCGTTGAAATCGACGGGTTTTCTCTGCAATaa > 1:166055/1‑90 (MQ=255) aTCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGAc < 1:396972/59‑1 (MQ=255) aTCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGAc > 2:396972/1‑59 (MQ=255) aTCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGACGGTTTCTCTCTGCAATAAACGCCGCCAGAcc > 2:500755/1‑90 (MQ=255) | AGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGAATCGCCGGGAAGAGCAGTTTTGCTCCGCGCACATGAATATGTTTACCGTTGAAATCGACGGTTTCTCTCTGCAATAAACGCCGCCAGACC > NZ_CP009273/990867‑991034 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |