Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,049,986 | G→T | 33.3% | noncoding (88/109 nt) | spf → | Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,049,986 | 0 | G | T | 33.3% | 31.2 / 18.5 | 21 | noncoding (88/109 nt) | spf | Spot 42 sRNA antisense regulator of galK translation, Hfq‑dependent |
Reads supporting (aligned to +/- strand): ref base G (7/7); new base T (4/3); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGC > NC_000913/4049843‑4050125 | aaaaaaaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATg > 4:802807/1‑149 (MQ=255) aaaaTACAAAAAGGTCTTTCCTAACTTAAACAAAAAGAATAAAATTAATCGGGGAAGGTTCAGGGGTTAGATGTTCTATCTTTTCAACCCTTTAATTCCCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg < 3:1000669/149‑1 (MQ=255) aaTTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg < 3:802807/149‑1 (MQ=255) gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag < 1:1693765/93‑1 (MQ=255) gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag > 2:1693765/1‑93 (MQ=255) tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg < 4:62790/99‑1 (MQ=255) tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg > 3:62790/1‑99 (MQ=255) ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa > 1:62571/1‑149 (MQ=255) tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc < 1:1001936/75‑1 (MQ=255) tCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc > 2:1001936/1‑75 (MQ=255) tCTTTCAGACCTTTTACTTCACGTAATCGCATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTCTTTTATTGGGCGAAAGAAAAGATACGTAATGCCTGCTGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa > 1:1130146/1‑149 (MQ=255) ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg < 3:100965/102‑1 (MQ=255) ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg > 4:100965/1‑102 (MQ=255) aTTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACt > 3:1213656/1‑40 (MQ=25) aTTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACt < 4:1213656/40‑1 (MQ=21) ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATc < 2:62571/149‑1 (MQ=255) cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg > 2:515620/1‑77 (MQ=255) cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg < 1:515620/77‑1 (MQ=255) ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAGTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGc > 2:892013/1‑149 (MQ=255) cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCgg < 2:504234/118‑1 (MQ=255) cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCgg > 1:504234/1‑118 (MQ=255) cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg > 2:50366/1‑146 (MQ=255) cccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAg < 1:50366/146‑1 (MQ=255) | AAAAAAATTACAAAAAGTGCTTTCTGAACTGAACAAAAAAGAGTAAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGC > NC_000913/4049843‑4050125 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |