| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000964 | 2,579,550 | 0 | T | G | 50.0% | 3.4 / 16.8 | 22 | G56G (GGA→GGC) ‡ | pstA | phosphate ABC transporter (permease) |
| Reads supporting (aligned to +/- strand): ref base T (5/6); new base G (11/0); total (16/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.24e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGTATAACCCCAGCCTGTTAAGTTGACAAACATCAATAATCCAAACATCCCAATGACAATAGACGGAAGAGATGAAAGCACCTCTATACATGTTCTAATGAAGTCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTTCATTACACTTTCCTCCTAGATGACAAGTAGCGAATAAGCAGGATGAACAAGAATG > NC_000964/2579322‑2579769 | cgaataaacacAGCCTGATAAGTTGACAAACATCAATAATCCAAACATCCCAATGAAAATAGACGGAAGAGATGAAAGCACCTCTATACATGTTCTAGTGAAGTCTGATACTTTGTTATTAGGCGCGTATTCGTCCATGAAGACTCCGCCGCCAACACCGTCCGGGACTGTAATCAGCATTGCAATAAACAAAATATAAAAGGAGTTAATCAGCTCATCTCGTATGCCTCCTCCAGCTTCAATTgcgcaa > 2:85624/11‑248 (MQ=255)cGTATAACCCCAGCCTGTTAAGTTGACAAACATCAATAATCCAAACATCCCAATGACAATAGACGGAAGAGATGAAAGCACCTCTATACATGCTCTAATGAAGTCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGccgccttctgcttgaaaaaaaaac > 1:250612/1‑232 (MQ=255) aaaGTCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGcacaca < 2:46733/250‑1 (MQ=255) aaGTCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAATAAAACTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTaa > 1:171710/1‑224 (MQ=255) tCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAAct < 1:118403/173‑1 (MQ=255) tCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAAct > 2:118403/1‑173 (MQ=255) cagatACTTTGTTATTAGGCGCGTATTCAGCCATTAAGACTCCGCCGCCTACAACGAGCGGTATTGTAATCAGCATTATAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTAGGCCTGCTCCAGCTTAAATGGTGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGACATACGCCGATTATAAGAATATAGGAAAACAAACCTTCTAAAATGACAGCGATAATGGATGCACAcaatct < 2:132029/246‑4 (MQ=255) aGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCGCACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAAAAAACCTACTAAAATGGCAGCGAAAATGGCTGCACAcaaacca > 1:159820/1‑250 (MQ=255) ttattaGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAATAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGt > 1:241789/1‑251 (MQ=255) aGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCCCACAAACCACACATTCCGGTCGcc > 2:109137/1‑250 (MQ=255) cgcgTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGCCTCCTGCTGCAATGACGCTCGACTTTGTTGTAATAAATTg > 1:124696/1‑147 (MQ=255) tATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAACACGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACAGTTTATAATAATATAGGAAAAAAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACAAAAAATTCCGGTCGCCAGTTTAtg > 1:316422/1‑250 (MQ=255) aaGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAACAAACGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACATCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCg > 1:286568/1‑251 (MQ=255) gAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTTCAGCTTGAATGGCGATCGACTTTGTTGTAGTAAATTGAAAACTGAGATGAGAAACACCGTTTATAATAATGTAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGcaa > 1:110353/1‑210 (MQ=255) gggATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTTCATTACACTTTcctcct > 2:57104/1‑247 (MQ=255) gggATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTTCATTACACTTTCCTCCTAGAt < 1:125486/251‑1 (MQ=255) tGTAATCAGCAATGTAATAAACAAAATATAAAATGAGTTAAACAGCGGGTATCGTATGCCTACGCCAGCTGTAATGACACACGACATTGTAGTAATAAAATGAAAACTAAGCAGAGAAACTCTGTTTATAATAATATAGGAAAACAGACCTACTAAAATGGCATCGATAATGGATGCACACACACCAAACATTCCGGTCGCCAGTTTATCTGTTAGTTTGCGGTTGATTacacgtttcgcccagaagaca < 2:101024/250‑18 (MQ=255) aaGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTTATCTCGTATGCCTCCTCAAGCTGCAATAGCGCTCGACTTAGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTACATTACACTTTCCTCCTAGATTACAAGTAGCg < 2:9392/250‑1 (MQ=255) aTTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAAAAATTTGATGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACATACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTACGGCCGCCAGTTTATCTGTTATTATGCGGATCATTACACTTTCCTCCTAGATGACAAGAAGCGAAt > 1:251720/1‑250 (MQ=255) taataaACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGCCTTCTGCTTCAATAACACTCGAATTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATTATAATATAGGAAAACAAACCTACTAAAATGGCAGCGAAACTGGCTGCACACAAACCATACATTCAGGTCGCCAGATTATCTGATATCTTGCGGTTAATTACAATTTACTACTAGATGACAAGTAGCGAATAAGc > 1:261637/1‑250 (MQ=255) aataaACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAGAACACGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTTCATTACACTTTCCTCCTAGATGACAAGTAGCGAATAAGc > 1:58220/1‑249 (MQ=255) aaaCAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAAATAGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTTCGGTTCATTACACTTTCCTCCTAGATGACAAGTAGCGAATAAGCAGGAt > 1:93075/1‑251 (MQ=255) aaaaGGAGTTAAACAGCTGATCTCGTATGCCGTCTTCTGCTTGAAAAAAACTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATACCGGTCGCCAGTTTATCTGTTATTTTCCGGTTCATTACACTTTCCTCCTAGATGACAAGTAGCGACTAAGCAGGATGAACAAGAATg > 1:134483/1‑251 (MQ=255) | CGTATAACCCCAGCCTGTTAAGTTGACAAACATCAATAATCCAAACATCCCAATGACAATAGACGGAAGAGATGAAAGCACCTCTATACATGTTCTAATGAAGTCAGTTACTTTGTTATTAGGCGCGTATTCAGCCATGAAGACTCCGCCGCCCACACCGAGCGGGATTGTAATCAGCATTGTAATAAACAAAATATAAAAGGAGTTAAACAGCTGATCTCGTATGCCTCCTCCAGCTGCAATGGCGCTCGACTTTGTTGTAATAAATTGAAAACTGAGCTGAGAAACACCGTTTATAATAATATAGGAAAACAAACCTACTAAAATGGCAGCGATAATGGCTGCACACAAACCAAACATTCCGGTCGCCAGTTTATCTGTTATTTTGCGGTTCATTACACTTTCCTCCTAGATGACAAGTAGCGAATAAGCAGGATGAACAAGAATG > NC_000964/2579322‑2579769 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |