Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,024,609 | 0 | T | G | 20.8% | 43.0 / 11.0 | 24 | intergenic (+18/+26) | rcsA/dsrB | transcriptional regulator of colanic acid capsular biosynthesis/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base T (7/12); new base G (3/2); total (10/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.15e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAGACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAACT > NC_000913/2024465‑2024747 | cAAGACCGTTTCATCGCATAAAGGTAATATTAAACATAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTcc < 1:16310/149‑1 (MQ=255) gACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCccac < 4:158351/149‑1 (MQ=255) tttCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGGAAACCACCAGTc < 7:84203/149‑1 (MQ=255) cATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGaa > 4:51107/1‑149 (MQ=255) tAATATTAAACGTAAGAGCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTa < 3:51107/149‑1 (MQ=255) tAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTtctc < 4:97814/149‑1 (MQ=255) aaaaCGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACaaa < 5:24931/148‑1 (MQ=255) aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTc > 1:95764/1‑96 (MQ=255) aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTc < 2:95764/96‑1 (MQ=255) aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCt > 8:98904/1‑128 (MQ=255) aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCt < 7:98904/128‑1 (MQ=255) aaaCAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATCCTGACTGGTGGGAAACCACCAGTCAgaa > 2:157477/1‑104 (MQ=38) aaaCAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATCCTGACTGGTGGGAAACCACCAGTCAgaa < 1:157477/106‑3 (MQ=38) gTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCt > 7:63186/1‑149 (MQ=255) aTCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGc < 8:63186/148‑1 (MQ=255) gtcCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGaa < 3:60871/148‑1 (MQ=255) gtcCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAAc > 8:19684/1‑149 (MQ=255) cTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCtt < 6:55152/149‑1 (MQ=255) aTGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGAt < 7:11595/129‑1 (MQ=255) aTGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGAt > 8:11595/1‑129 (MQ=255) acaTTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAgg > 3:84812/1‑133 (MQ=255) acaTTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAgg < 4:84812/133‑1 (MQ=255) tGACTGGTGGGAAACCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAAct > 1:86667/1‑149 (MQ=255) tGACTGGTGGGAAACCACCAGTCAGGCTGAATAAGATTACTATGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGCACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAAct > 4:90669/1‑149 (MQ=255) gtggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATtgt < 5:74577/129‑1 (MQ=255) tggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACAtt > 2:5381/1‑126 (MQ=255) tggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACAtt < 1:5381/126‑1 (MQ=255) | CAAGACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAACT > NC_000913/2024465‑2024747 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |