| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,978,437 | G→A | intergenic (‑240/+81) | flhD ← / ← insB5 | DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,978,437 | 0 | G | A | 93.8% | 45.7 / ‑4.1 | 17 | intergenic (‑240/+81) | flhD/insB5 | DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (6/9); total (8/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
ACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGC > NC_000913/1978205‑1978683 | ccccaGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCaata < 1:295485/250‑1 (MQ=255) aCCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATg < 1:253309/251‑1 (MQ=255) ggTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACATTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGACTTTACACTTATCTAAGATTTTTCATAAATCGACGCAACTGTACTCGTCCCTAAACGCACATACAACGGAGGTGGGCTGCGTTTTTCAATAATGCGTGATGCAGATCAC‑CACAAACACTc > 2:367880/1‑251 (MQ=255) ccGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AACCACttagactgg > 2:127639/1‑243 (MQ=255) gCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGATGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTa < 1:55741/251‑1 (MQ=255) gTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATaa > 1:252363/1‑250 (MQ=255) ccGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataa < 2:76884‑M1/249‑6 (MQ=255) aatTAGCATTAGAATAGGTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTGTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat < 1:367880‑M1/251‑30 (MQ=255) tGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataag < 1:160684‑M1/200‑7 (MQ=255) tGCAATAAGCAGAAACACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATCCAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataag > 2:160684‑M1/1‑194 (MQ=255) gAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgcatttcaatctgttaaaaagtttttcgctacgctagcaagc > 1:379116‑M1/1‑183 (MQ=255) tttGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac > 2:149845‑M1/1‑173 (MQ=255) tAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGTGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt > 1:328880‑M1/1‑166 (MQ=255) tAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaatctgttaaaaagttt < 2:328880‑M1/215‑50 (MQ=255) aagcgacGCAACTGTACTCGTCACCACACGCTCATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGAAGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaacaatcatctgc < 2:176605‑M1/248‑141 (MQ=255) cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca > 2:304003‑M1/1‑108 (MQ=255) cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatactttggtgttgggtgtatttca < 1:304003‑M1/141‑34 (MQ=255) ggaggaGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaatcatcatctgctataaataatcactttcatgcaataccagataagctattt > 2:364629‑M1/1‑71 (MQ=255) gaggagGGCTGAGGTTTTCAATACTGCGTGATGCAGCTCACACA‑ACACACTCAGTTACTTAACATAAATGtataagtcatacttttgttatgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaagagtgaaacaggaataatcgaagtggaatgttgcgcacagtcaaagtaactcaccgtaattaatcatctgctataaataatcattttcatgcaataccagataagctatttt > 2:133758‑M1/1‑70 (MQ=255) | ACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACA‑AAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGC > NC_000913/1978205‑1978683 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |