Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,476,586 | 0 | T | C | 20.0% | 98.3 / 21.0 | 50 | intergenic (‑146/+21) | yheO/fkpA | putative PAS domain‑containing DNA‑binding transcriptional regulator/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase) |
Reads supporting (aligned to +/- strand): ref base T (19/21); new base C (5/5); total (24/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTG > NC_000913/3476463‑3476719 | gAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAtt > 1:48271/1‑149 (MQ=255) gAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCGGGTCAGGCGGCGGTTCTTAATGCTTa < 7:19113/147‑1 (MQ=255) gAGCCTGTCGTTAACGTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTa > 8:19113/1‑147 (MQ=255) ccTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAttttt < 8:21202/149‑1 (MQ=255) ccTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTAttttt < 2:82593/149‑1 (MQ=255) tCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa < 2:48271/147‑1 (MQ=255) gTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGTtttttta < 1:54524/133‑4 (MQ=255) gTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGACCAGGCGGCGGTtttttta > 2:54524/1‑130 (MQ=255) ttttCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTg > 6:14657/1‑149 (MQ=255) gggCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGc < 8:85895/123‑1 (MQ=255) gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCttt > 8:91745/1‑148 (MQ=255) gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCttt < 7:91745/148‑1 (MQ=255) gggTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTg > 3:99202/1‑149 (MQ=255) aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa < 4:42393/103‑1 (MQ=255) aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa > 3:42393/1‑103 (MQ=255) aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCa < 2:35091/106‑1 (MQ=255) aaTTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCa > 1:35091/1‑106 (MQ=255) tACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGgcg > 8:5233/1‑149 (MQ=255) cAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGgcgc < 5:14657/148‑1 (MQ=255) gTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTa < 8:74805/80‑1 (MQ=255) gTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTa > 7:74805/1‑80 (MQ=255) ttCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa < 1:182495/86‑1 (MQ=255) ttCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTa > 2:182495/1‑86 (MQ=255) ccGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATc > 5:83546/1‑149 (MQ=255) cGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATc < 4:99202/148‑1 (MQ=255) gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCtt > 2:54269/1‑129 (MQ=255) gCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCtt < 1:54269/129‑1 (MQ=255) tAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGc > 5:99162/1‑81 (MQ=255) tAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGc < 6:99162/81‑1 (MQ=255) ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc > 2:178122/1‑119 (MQ=255) ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc < 1:178122/119‑1 (MQ=255) ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAgg < 6:83546/149‑1 (MQ=255) aataaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCttt < 1:157583/88‑1 (MQ=255) aataaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCttt > 2:157583/1‑88 (MQ=255) attaagaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc > 4:97909/7‑127 (MQ=255) attaagaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTAc < 3:97909/121‑1 (MQ=255) aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTc > 4:2468/1‑129 (MQ=255) aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTc < 3:2468/129‑1 (MQ=255) aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGaa < 8:102742/143‑1 (MQ=255) aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGaa > 7:102742/1‑143 (MQ=255) agaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa > 5:128276/3‑139 (MQ=255) agaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa < 6:128276/137‑1 (MQ=255) gaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg < 4:94705/147‑1 (MQ=255) gaaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg > 3:94705/2‑148 (MQ=255) ccgccgCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg > 1:67410/1‑145 (MQ=255) ccgccgCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg < 2:67410/145‑1 (MQ=255) cgccgcCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc > 1:65919/1‑101 (MQ=255) cgccgcCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTc < 2:65919/101‑1 (MQ=255) gccgccTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg < 2:162652/129‑1 (MQ=255) gccgccTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg > 1:162652/1‑129 (MQ=255) | GAGCCTGTCGTTAACTTTTCAGGGCAGGCTCATTAATGATGCGGGTAACTAAATTAATACAGCGGAGGTTCCGCTTTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTG > NC_000913/3476463‑3476719 |
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: ‑ |