Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 278,742 | 1 | . | A | 63.6% | 41.2 / 26.9 | 22 | intergenic (‑237/+77) | PP_0224/sctC | DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein |
Reads supporting (aligned to +/- strand): ref base . (2/6); new base A (7/7); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.80e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTC > NC_002947/278605‑278879 | aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctc < 4:52847‑M1/148‑9 (MQ=255) tCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctac < 3:70154‑M1/146‑13 (MQ=255) gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctacccacat > 8:15886‑M1/1‑130 (MQ=255) gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa < 6:75866/148‑1 (MQ=255) tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACCCATGCCTCAAGAt > 1:149534/1‑148 (MQ=255) gCTGGTCCGGTGGTCCGGTATATATACACTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATc > 7:33027/1‑148 (MQ=255) ggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctacagtgcacagatcttgaggcatgtgtgtaggagc < 5:7388‑M1/148‑46 (MQ=255) tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctacagtgcacagatcttgaggcatgtgtgtaggagcgg < 1:113022‑M1/148‑48 (MQ=255) tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctacagtgcacagatcttgaggcatgtgtgtaggagcgg < 3:30698‑M1/148‑48 (MQ=255) atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGTGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt < 6:2646/147‑1 (MQ=255) atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt < 1:15653/147‑1 (MQ=255) ttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTTCACACATGGCTCAAGATCTGTGCCCTGTAGGAGCGGGTTTACCCGCGCAAGGGCCGGGACAgggt > 5:139466/1‑146 (MQ=255) tATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctccta > 2:108765‑M1/1‑71 (MQ=255) aTGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AAacccgctcctac < 1:108765‑M1/82‑13 (MQ=255) cctaactTTTTGGAATATTAATAGAAGGGGGATTTGTTACCGAGGGCCCGGGCTCTGCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg < 5:160850/147‑1 (MQ=38) ctaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTTTGCACTGTAGGCTCCGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc > 2:129859/1‑146 (MQ=255) taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa < 6:44040/148‑1 (MQ=255) tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc > 7:89906/1‑148 (MQ=255) tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGGGCCGGGCTCTGCGCGGGTAAACCCGCTCCTACACACATCCCTCAAGATCTGGGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc < 7:42905/148‑1 (MQ=39) ttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCa > 1:12418/1‑148 (MQ=255) tttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCCGCACACATGCCTGAAGCACTGTGGACCGGAGGCTCGGGGTTCCCCGCGCTAGGGCCGGGACCGGGTGACGCGGAGGCCGCCGTCACGCGGCCCTTTCGCacccc > 7:150920/1‑144 (MQ=255) ctTCGCGGTTAAAACAGCTCCTACAAACATGCCTCAAGATCTGTCCACTGTAGTAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc < 8:110338/148‑1 (MQ=255) | AGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGT‑AACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTC > NC_002947/278605‑278879 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |