Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,024,542 | 0 | G | A | 39.4% | 47.2 / 10.7 | 33 | H94Y (CAT→TAT) | rfaH | transcription antitermination protein |
Reads supporting (aligned to +/- strand): ref base G (14/6); new base A (0/13); total (14/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.04e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTCACTG > NC_000913/4024392‑4024684 | tAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCtgatg > 2:152754/1‑151 (MQ=255) gCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAAt < 2:383889/151‑1 (MQ=255) gCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAAt < 1:474420/151‑1 (MQ=255) aaTAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATaa < 2:342523/151‑1 (MQ=255) cATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAAccgcc < 2:569489/151‑1 (MQ=255) ggAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGTTAATAA‑TCCCCTTATCCCCCGG‑ATAAGGGGTTGCCGGACCGACAA‑TGTCTTTCGGTTT‑ATATCCCG‑ATAGCTGATGAATAAccgccg < 1:534576/149‑1 (MQ=255) agcCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACt > 2:505923/1‑151 (MQ=255) ccgtggggtggggTAAAAAGGGCCGAAAACCC‑TTCGACCGCGCCTTCGATAATAA‑CCCCCTTACCCCCGGG‑ATAGGGGGTTGCCGGACCGA‑AA‑TTTTTTTTGGTTTAATATCCCG‑ATAGCTGATAAATAACCGCCGATGGGACTATCGCTg < 2:544467/140‑1 (MQ=255) aaaaTGCCCTAAACCCC‑TCGGAACGCCCTTTCGTTAATAA‑TCTCTTTACCCCCCGA‑AAAGGGGGTGGCCGTAACGAAAA‑TGTTTTTGGTTTT‑TTATACGG‑ATAGCTGAAAAATTTCCGCCGATGGGACTATCGCTGGCGACGCGCCAAAgc < 2:412370/151‑1 (MQ=255) aaaTGGCCGGAAACCCGTTCG‑ACCCCCCTTTGGTAAAAACTCCCCTCA‑CACCCGGAATAGGGGGTTCCCGAACCTACAC‑TTTCTTGCTGTTT‑ATATCCCG‑ATAGCTGCAAAAGAACCGCCGATGGGACTATCGCTGGCGAcgcg < 1:255056/144‑1 (MQ=255) cc‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTg > 2:658405/1‑150 (MQ=255) c‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGac > 2:40725/1‑151 (MQ=255) gcgcCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCggg > 1:213487/1‑151 (MQ=255) ccttgtaaaataa‑tCACCTTACCGCCCAG‑GAAAGGGGTTGCCCGATCGACAA‑TGTCTTTCTGTTT‑ATATACCG‑ATAGCTGAGAAATAACCGTCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCggg < 1:142480/140‑1 (MQ=255) gggaaaaaa‑tCCCCTTACCCCCGGG‑ATAGGGGGTGGCGGGACCGACAT‑TTTTTTTGGTTTT‑ATATCCCG‑ATAGCTAATAAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGtt < 2:564555/144‑1 (MQ=255) taa‑tCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAg > 1:95606/1‑150 (MQ=255) actCACCTTA‑CCCCCGG‑TTAAGGGGTGTCCAGCTCAAGAAGTGTTTTTCTTTTA‑TTATACGG‑ATAGCTAAGAAAAACCCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATa < 1:694591/146‑1 (MQ=255) ccccttacccccGG‑GAAAGGGGGTGGCGGAACCAAAA‑TTTTTTTGCGTTT‑AATAACCGAATAGCTAAAAAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTg < 2:660044/142‑1 (MQ=255) ccccgggagagggggTTGCCGGCCCGAAAT‑TTTCTTTGGTTTT‑ATATACCG‑ATAGAGAAAAAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATgt < 2:616669/138‑2 (MQ=255) ccccgggataGGGGGTGCCGGGTCCAAAATTTTTTTTCGTTTT‑ATACCCGA‑ATAGCTAATAAAAAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATg < 2:72327/139‑1 (MQ=255) aaGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTgg > 1:96346/1‑151 (MQ=255) ggggTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGt > 1:431126/1‑151 (MQ=255) ggATCAAAAT‑TGTTTTTGGGTTT‑ATACCCCG‑ATAGGTAATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCca < 2:184584/151‑1 (MQ=255) tCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCcaca > 2:322846/1‑150 (MQ=255) tCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCcaca > 2:177048/1‑150 (MQ=255) cGACAA‑TGTCTTCCGGTTT‑ATATACCG‑ACAGCTGATGAATAACCGCCGATCGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACaa > 1:576565/1‑150 (MQ=255) ttgtggTTA‑ATATCCGG‑ATAGGTAAGAAAAACCCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTa < 1:581279/143‑1 (MQ=255) ttggtttt‑ATATACGG‑ATGCCGGATAAAAAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTg < 2:538840/145‑1 (MQ=255) gTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGa > 2:36772/1‑150 (MQ=255) ttt‑ATATACCG‑ATGGCTAATAAAAAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTgg < 2:610801/146‑1 (MQ=255) ttt‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAAc > 1:133754/1‑151 (MQ=255) g‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTCac > 1:250618/1‑151 (MQ=255) aTAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTCact > 1:374215/1‑151 (MQ=255) gttgataaATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTc < 2:68786/143‑1 (MQ=255) gCTGAGAAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTCactg < 2:224430/149‑1 (MQ=255) | TAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCC‑TTCGAACGCGCCTTCGGTAATAA‑TCACCTTATCACCCGG‑ATAAGGGGTTGCCGGATCGACAA‑TGTCTTTCGGTTT‑ATATACCG‑ATAGCTGATGAATAACCGCCGATGGGACTATCGCTGGCGACGCGCCAAAGCGCACGAAGTGGCTGACACCGCGGGTCGCGTTGATAGTCGTGGTATGAATCACTTCTGGGTCAAATTCCACAAACAGGTAGTTGGGGAACAATGGCTCACTG > NC_000913/4024392‑4024684 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |