Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,623,312 | 0 | A | T | 36.4% | 18.8 / 18.3 | 22 | H159Q (CAT→CAA) | PP_4941 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (8/6); new base T (4/4); total (12/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAA > NC_002947/5623195‑5623458 | cAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag < 2:266829/130‑6 (MQ=255) cAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag > 1:266829/1‑125 (MQ=255) aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag > 1:123127/1‑119 (MQ=255) aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag < 2:123127/124‑6 (MQ=255) aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTGCCGGCtgcag > 1:122968/1‑122 (MQ=255) aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTGCCGGCtgcag < 2:122968/124‑3 (MQ=255) gcgcGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCa > 1:87067/1‑150 (MQ=255) atcaccatcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCa < 2:257623/150‑1 (MQ=255) caccatcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATc < 2:87067/150‑1 (MQ=255) catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTgg > 2:1078/1‑108 (MQ=255) catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTgg < 1:1078/108‑1 (MQ=255) catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAg > 2:228311/1‑126 (MQ=255) catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAg < 1:228311/126‑1 (MQ=255) atcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACg > 1:270026/1‑150 (MQ=255) gACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATAGACGCCACGCTCGAACGCTTGCCGATc > 2:73006/1‑150 (MQ=255) ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgccagcca < 2:114866/135‑1 (MQ=255) ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgccagcca > 1:114866/1‑135 (MQ=255) ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAg > 1:300656/1‑150 (MQ=255) gcggCCGGCATTGCCGGCTGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTc > 1:288154/4‑124 (MQ=255) gcggCCGGCATTGCCGGCTGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTc < 2:288154/121‑1 (MQ=255) ggCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGaa > 1:128042/1‑150 (MQ=255) gCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGaaa < 2:300656/150‑1 (MQ=255) | CAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAA > NC_002947/5623195‑5623458 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |