Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,559,904 | 0 | C | G | 52.9% | ‑0.4 / 11.7 | 17 | intergenic (‑302/+29) | PP_4888/purA | methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase |
Reads supporting (aligned to +/- strand): ref base C (4/4); new base G (0/9); total (4/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.94e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCA‑GGG‑CC‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC > NC_002947/5559801‑5560039 | gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTAGCA‑GGG‑GC‑GCGGGGGCACATAAGCCGGGGAGGCAGaaaaa > 2:1800486/1‑132 (MQ=25) tttACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTCA‑GGG‑CC‑CCGGCCGAAGAAACGCCAGGAACGCAGAACAgg > 2:2859708/1‑136 (MQ=37) tACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGGCCCTGGGA‑GGG‑GC‑CCGGCGGCACGTCAGCCGGGCAGGGCGATGAGGcc > 2:3072410/1‑135 (MQ=21) acaAATGAAAACGCCGCGGCCCATCTCA‑GGG‑CC‑GCGGCGTAAGAAAAGCCAGGAAGGCAGATAAGGCGAAACGGTGAAGCAGCACGATTGTCTCGTTGCGGTCCGGGCCGGTCGAAATTATGTAGATAGGCGCGcc > 2:3123389/1‑136 (MQ=255) taaacccgccgcGGCCTTTCTCC‑GGG‑CG‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCa < 2:908638/130‑1 (MQ=255) ccccggGGCCCTCCCCG‑GGG‑CG‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt < 2:47691/134‑1 (MQ=255) ggcggggCCCTTTTCGGGGG‑CG‑GGGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt < 1:2859708/131‑1 (MQ=255) gcgggggCCTTTGTCG‑GGG‑CC‑GGGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTc < 1:3318294/131‑1 (MQ=255) cgccgggCCCTTCCCC‑GGGCCC‑CCCGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt < 2:332343/131‑1 (MQ=255) gggggCCTTCCTCG‑GGC‑CG‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGa < 2:2811958/133‑1 (MQ=255) ggggCCCTTGCCG‑GGG‑CCGGCCGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGa < 2:4649847/134‑1 (MQ=255) gcccgttCTCG‑GGG‑CG‑GGGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATgc < 1:1800486/131‑1 (MQ=255) ccTTCTCG‑GGG‑CG‑GGGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCt < 2:4698741/136‑1 (MQ=255) ttCTCG‑GGG‑CG‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTg < 1:3123389/136‑1 (MQ=255) ttCTCG‑GGC‑CG‑GGGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTg < 2:4678755/136‑1 (MQ=255) tccgcgggg‑CG‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTg < 1:2625932/130‑1 (MQ=255) gccgcggg‑CCGGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTg < 1:174200/131‑1 (MQ=255) cc‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTCATGTAGTCa > 1:2080736/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt > 2:131577/1‑115 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCtt < 1:131577/115‑1 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGAt < 2:3110046/128‑1 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGAt > 1:3110046/1‑128 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 2:2145616/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 1:3864902/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 2:3758759/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 2:4130565/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 1:3075216/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 1:2890056/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 1:540181/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc < 1:1916590/136‑1 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc > 1:872305/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCATCGATGCGCTTGATGTAGTCAc > 2:2236404/1‑136 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCATCGATGCGCATGATGTAGTCaa > 2:2259244/1‑135 (MQ=255) c‑gcgGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGAAGGTCTCGTTGCGGTCCGGGCCGGTCGAATTGATGTCGATCGGCGCGCCAGCCCGCTCTTCGATGCGCTTGATGTAGTCAc > 1:3797635/1‑136 (MQ=255) c‑gcgGCGTCAGAACAGCCAGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTCCGGACCGGGCCGGTCGAAAAGATGTCGATCGGCGCGCCAACCAGCTCATCGAAGCGCTAGATGTAGTCAc > 1:4505621/1‑136 (MQ=255) c‑gccgcGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAc < 1:4077088/136‑1 (MQ=255) | GTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCA‑GGG‑CC‑GCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAC > NC_002947/5559801‑5560039 |
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: ‑ |