Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 524,167 | 0 | A | G | 42.9% | 11.5 / 12.4 | 29 | intergenic (+34/‑782) | tyrS/PP_16SC | tyrosine‑‑tRNA ligase/16S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (1/11); new base G (9/0); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.40e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGA > NC_002947/524020‑524291 | gTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGTCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGa > 7:16350/1‑148 (MQ=255) ttGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GCCggg > 2:148928/1‑147 (MQ=255) aGTGGTCGATAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCTACCG‑GCCT‑ccc < 2:76394/147‑1 (MQ=255) aGTGGTCGACAAGGACTTCATGTTTGGCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGCTGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCTACCG‑GCCT‑ccc < 8:101803/147‑1 (MQ=255) gggcGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGGCGTTTGGCCGGGTTACCCTCACGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCG‑CCCCGC‑t < 4:14430/145‑1 (MQ=255) aCAAGGACTTCATGTTTGTCCTTGGTGCTACCCCTGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCGTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑tt < 5:127962/148‑1 (MQ=255) aaGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGG‑GGCGGGAG‑GACT‑CCCCGC‑tttt‑tttt > 4:36988/1‑148 (MQ=255) tcTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCACGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTg < 7:117904/145‑1 (MQ=255) cTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GGCGGTAG‑GCCA‑CCCCCC‑GTTG‑TTTTTtggg > 6:133642/1‑148 (MQ=255) cATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCG‑CCCCGC‑GTTG‑TTTTTGTGGgtt > 4:21189/1‑148 (MQ=255) aTGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGG‑GCCGGTAG‑GCCTCCCCCGG‑TTTG‑TTTTTTGGtttg > 6:157324/1‑147 (MQ=255) tGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGTCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTtg < 8:16350/147‑1 (MQ=255) gTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGG‑GCCGGTGG‑GCTG‑CCCCGC‑GTTT‑gttttgggttttgg > 8:8536/1‑148 (MQ=255) gTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GCCCGTTG‑GCCT‑CCCCGCGTTTG‑TTTTTGGGGttggt > 5:110955/1‑145 (MQ=255) ccATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GGCCGGAG‑GGCG‑CCCCCC‑GTTT‑TTTTTTGGGTTTGGTGGGTGGGGTTGTtggggg > 3:15631/1‑148 (MQ=255) cATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGG‑GGCGGGAG‑GGCG‑CCCCGC‑GTTG‑TTTTTGGGGGTTGGGGGTTGGGGGTGTtgggggg > 1:13304/1‑147 (MQ=255) aTGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GCCGGGAG‑GCCG‑CCCCGC‑GTTG‑TTTTTTGGGGTTGGTGGTTGGTGTTGTtggggggg > 7:122642/1‑146 (MQ=255) gtgtgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGG‑GCCGGTGG‑GCGG‑CCCCCA‑TTTG‑TTTTT‑GGGTTTTGTGGTTTGGGTTGTtggggggggg > 1:124698/1‑148 (MQ=255) cAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGCG‑GGCGGGAG‑GGCT‑CCCCCC‑GGTG‑TGTTTGGGGGTTTGGGGGTGGTGGTGGAGGGGGGGGGGGcggg > 3:34069/1‑145 (MQ=255) aagaagTCGTTTGGCCGGGTTACCCTCAAGGCTTAGTGATTGCCTGGAAGGTGTAGCTATATAGAAGC‑GGGGGGCGCGCGTCGCCCCC‑CCCCGC‑TGTG‑‑TTTTGTGGTTTGGGGGGGGGGGTGGAAGGGGGCGGGGGCAGTTTT‑‑AGAgg > 2:99641/1‑148 (MQ=255) agaagTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGG‑GGCGGGAG‑GGCG‑CCCCCC‑GTTG‑TTTTTTGGGGTTTGTGGTTGGGGTTTTTGGGGGGGGGGGGGGtttttgggagg > 6:44860/1‑141 (MQ=255) ccttgggCCGGGGTACCCTCAAGGCTGAGGGATTGCCTGGATGGTGTCGCTATATCGAAGC‑GGGGCG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGGTGCTGTTGTAGGGGGAGGGCGCTGTTTTG‑AGAGGGTGTTgg < 4:53902/143‑1 (MQ=37) tttGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGGG‑GGCGGGTG‑GGCG‑GCCCGC‑GGTGTTTTTTTGGGTTTGG‑GGGTGGGGTTGTTTGGGGGGGGGGGGGGTTTG‑TGTGGGTGtttggg > 5:45906/1‑144 (MQ=255) ggCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGACGC‑GGGGAG‑GCCGGTAG‑GCCG‑CCCCGC‑TTTG‑TTTTTGGGGTGTGGTGGGTGCTGTTGTAGGGGGCGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGt < 7:68536/147‑1 (MQ=39) ggCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGG‑GCCGGT‑G‑GCCA‑CCCCCCGTTTT‑TTTTTGGGGTTTGTGGGTTGGGTTTTTGGGGGGGGGGGGGGGTTTTGTAGGGGGT‑TTGGCggg > 2:121792/1‑147 (MQ=255) cGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTAGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTgaga < 2:33728/148‑1 (MQ=255) gTGATTGCCTGTATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTg < 2:50888/148‑1 (MQ=255) gTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTg < 1:50101/148‑1 (MQ=255) aTGGGGTAGCTATATAGAAGC‑GGGGCG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGGTGCTGTTGTAGGGGGAGGGGGCTGGTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTga < 7:8536/148‑1 (MQ=255) ggTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTgaga < 1:148928/148‑1 (MQ=255) | GTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGC‑GGGGAG‑GCCGGTAG‑GCCT‑CCCCGC‑TTTG‑TTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTG‑AGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGA > NC_002947/524020‑524291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |