Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,798,416 | (T)8→7 | coding (61/936 nt) | rfaS ← | LPS core biosynthesis protein RfaS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,798,409 | 0 | T | . | 80.0% | 17.3 / 3.4 | 10 | coding (68/936 nt) | rfaS | LPS core biosynthesis protein RfaS |
Reads supporting (aligned to +/- strand): ref base T (2/0); new base . (3/5); total (5/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01 |
CCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAT > NZ_CP009273/3798320‑3798498 | ccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACt < 1:662968/90‑1 (MQ=255) cAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtt < 1:833620/90‑1 (MQ=255) gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca > 1:742534/1‑88 (MQ=255) tATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGt > 1:320222/1‑90 (MQ=255) tttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAAt < 1:388356/75‑1 (MQ=255) tttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAAt > 2:388356/1‑75 (MQ=255) ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATaa < 2:621206/90‑1 (MQ=255) tttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaataa > 1:842019/1‑90 (MQ=255) tgttgtTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTc < 1:709312/90‑1 (MQ=255) tgttATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCata > 1:888308/1‑90 (MQ=255) ttATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCataat < 2:742534/90‑1 (MQ=255) ttATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagttagt > 2:529488/1‑90 (MQ=255) aTCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTc < 2:75527/90‑1 (MQ=255) ttttttttCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAt < 2:965463/90‑1 (MQ=255) | CCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAT > NZ_CP009273/3798320‑3798498 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |