| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,171,882 | 0 | C | G | 47.1% | 39.0 / 22.7 | 34 | intergenic (+126/‑175) | rrfB/murB | 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding |
| Reads supporting (aligned to +/- strand): ref base C (18/0); new base G (0/16); total (18/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.54e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAAC > NC_000913/4171735‑4172026 | agagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCgg > 1:353795/1‑150 (MQ=31) tAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAggg > 1:142720/1‑150 (MQ=31) aCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGGGGCggg > 2:475302/1‑151 (MQ=14) cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGCGGCggg > 1:411103/1‑150 (MQ=14) cAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg > 2:563706/1‑150 (MQ=31) tAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt > 2:322351/1‑150 (MQ=31) aaaaCGAAAGGCTCAGTCGAAAGACTGGGCTTTTCTTTTTATCTGTTTTTTGGGGGTGACCGCTCTCCGGAGTAGAAAAAACCCGCGGGGGGGGGATTTGAACG‑TTGCAAAGCAACGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATaa < 2:30770/151‑1 (MQ=11) aaaCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa > 2:463122/1‑150 (MQ=31) gCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAgg > 1:691270/1‑150 (MQ=31) cccaggaaaaaGGCTGGCCTTTTCGTTTTATTGGTTTTTTGGGGGTGACCCCTCTCCGGAGAGGAACAACCCCGCGGGGGGGGAATTTGAAGG‑TGGCGAAGAAGCGGCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCa < 2:606182/143‑1 (MQ=2) cTCAGTAAAAAGACTGGGCTTTTCTTTTTATCTTTTTTTTGTCGGTGAACGCTCTCCTGAGTGGGAAAACCCCGCGGGGGGCGGATTTAAACG‑TTGCGAAGAAACGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGc < 1:219535/150‑1 (MQ=2) tCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAt > 2:786492/1‑151 (MQ=31) aGTCGAAAGACGGGGCTTTTTTTTTTATCTGTTGTTTGTCGGTGACCGCTCTCCTGAGTAGAAAAAACCCGCGGGGAGGGGTTTTGAACG‑TTGCGAAGCAACGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATc < 1:459767/150‑1 (MQ=11) ggcGAAAGACGGGCCTTTTCTTTTTATTTTTTTTTTGTGGGAAACCCCCCTCGGGAGGAGGACAACCCCGCGGGGGGGGAATTTGACGG‑TTGGGAAGCACGGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCa < 2:765267/145‑1 (MQ=2) aaaaGCCTGGGCTTTTCTTTTTTCCTGTTTTTTGTGGGGGAACGCTCTCTTGAGAAGGAAAAACCCGCCGGGAGCGTATTAAAACG‑TTGCAAACCAACGGCCGGGAGGGTGGCGGGGAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt < 1:541225/150‑1 (MQ=2) gTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGt > 1:798819/1‑151 (MQ=31) tGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGtt > 1:456179/1‑151 (MQ=31) gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGCGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCt > 2:470991/1‑149 (MQ=18) gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc > 1:730428/1‑151 (MQ=35) ctcCTGAGAAGGACAAACCCGCGGGGGGCGGATTTGACCT‑TTGCGAAGCAGCGGCCGGGGGGGGGGGGGGCAGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc < 2:301022/148‑1 (MQ=11) ctggagaAGGAAAAACCCGCGGGGAGGGTATTGGAACG‑TGGGAAAGCAACGGCCGGGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc < 1:606964/147‑1 (MQ=14) cagaagtagtAAAAACCCGCGGGGGGCGGATGTGAGCG‑TTGGGAAGCAAGGCGCGGGAGGGGGGCGGGAAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc < 2:346200/146‑1 (MQ=14) ggggaagaacaatcccgcggGGAGTGTATTT‑AACGTTTGCGAACCACGGGGTGGGGGGGGGGGGGGCAGGAGCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActc < 2:447706/136‑1 (MQ=9) tAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt > 1:326289/1‑151 (MQ=255) aaaCCCGCCGGGAGGGGATTTGAAGG‑TTGCAAAGCAACGGCCGGGGGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTcc < 1:785675/144‑1 (MQ=14) aTCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCat > 1:319705/1‑151 (MQ=255) cgccgGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATc > 1:205235/1‑151 (MQ=255) ggAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACa > 2:22432/1‑150 (MQ=255) tttGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc > 2:410396/1‑151 (MQ=255) tttGAACG‑TTGCAAAGAAACGGCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCa < 1:315727/147‑1 (MQ=18) ttGAACG‑TTGAGAAGCACCGGCGGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccc < 1:466333/150‑1 (MQ=18) atggttggaaaGCAACGGCCGGGGGGGGGGGGGGCGGGACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCca < 2:497566/142‑1 (MQ=21) ttGCAAAGAACCGGCCGGGGGGGGGGGGGGCGGAACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcaca < 1:507366/149‑1 (MQ=21) ggCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAAc < 1:688998/149‑1 (MQ=255) | AGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACG‑TTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAAC > NC_000913/4171735‑4172026 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |