| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,555,781 | 0 | G | C | 20.0% | 61.6 / 13.7 | 31 | intergenic (‑112/+45) | maeA/sra | malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein |
| Reads supporting (aligned to +/- strand): ref base G (9/15); new base C (3/3); total (13/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTT > NC_000913/1555634‑1555918 | ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCGTGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgcg > 7:135201/1‑148 (MQ=255) aaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCACT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAgg > 6:100360/1‑149 (MQ=255) aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTTTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg < 3:52556/149‑1 (MQ=255) aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg < 8:137826/149‑1 (MQ=255) aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg < 3:129926/149‑1 (MQ=255) cGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCGCTGGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATgg < 5:50938/149‑1 (MQ=255) tttttttGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGtt < 2:154649/147‑1 (MQ=255) ctcactcTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg > 1:74799/1‑136 (MQ=255) ctcactcTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg < 2:74799/136‑1 (MQ=255) atatCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCa < 8:78911/116‑1 (MQ=255) atatCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCa > 7:78911/1‑116 (MQ=255) ccatccCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGttttt > 2:182866/1‑104 (MQ=38) ccatccCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGttttt < 1:182866/106‑3 (MQ=38) catccCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGgtgt < 3:125804/101‑1 (MQ=39) catccCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGgtgt > 4:125804/1‑101 (MQ=39) tttACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGc > 3:143363/1‑115 (MQ=255) tttACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGc < 4:143363/115‑1 (MQ=255) aCCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCt < 8:135201/148‑1 (MQ=255) ccAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg > 1:129980/1‑80 (MQ=255) ccAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg < 2:129980/80‑1 (MQ=255) cAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAttt < 5:100360/149‑1 (MQ=255) taaattataAAGCAGATAAAAACAAAAAAAACATTGCTCAGGCAATGTTGTTTAATCTTCATTGAGGACTGATGTTTATGAATTAATTTTCAGCGTGGAGTTTTCTGACGTTTGGTTTATTTACTACGCTGGATTTGTCACATTAGGt > 4:15473/6‑148 (MQ=255) aaGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGAt < 6:39943/71‑1 (MQ=255) aaGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGAt > 5:39943/1‑71 (MQ=255) aaGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAAc > 7:106774/1‑149 (MQ=255) aaGCAGATAAAAACAAAAACACCTTTGCGCAGGCATTGGTGTTTAATCGTCATTGAGTCCTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTAtttt < 3:15473/148‑1 (MQ=255) aGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGtttt > 6:21617/1‑94 (MQ=255) aGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGtttt < 5:21617/94‑1 (MQ=255) aGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACATTCGTTAACTATTTTgc > 1:139295/1‑149 (MQ=255) cAGATAAAATCAAAAACATCATTGCGCAGGCAATTGTGTTTAATCGTCATTTAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCTTTTTCTGCCGGTTGTGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcgc < 2:139295/149‑1 (MQ=255) acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAtttt > 7:38771/1‑148 (MQ=255) acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAtttt < 8:38771/148‑1 (MQ=255) | ATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCT‑GGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTT > NC_000913/1555634‑1555918 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |