Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 115,797 | 0 | C | A | 41.9% | 33.8 / 10.2 | 31 | E435* (GAA→TAA) | hofB | T2SE secretion family protein; P‑loop ATPase superfamily protein |
Reads supporting (aligned to +/- strand): ref base C (10/8); new base A (0/13); total (10/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.28e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACAC > NC_000913/115649‑115946 | ctGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTc < 2:578543/149‑1 (MQ=255) tGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaa > 2:44112/1‑150 (MQ=255) ccATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaaaa > 1:866660/1‑149 (MQ=255) aGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTa < 1:2995011/151‑1 (MQ=255) gcgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGcc < 2:153325/150‑1 (MQ=255) cgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGcc < 1:2606192/149‑1 (MQ=255) gCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgc > 2:3640807/1‑151 (MQ=255) ccGGTTATGCCATGCCAGCGCCAGGGTTGTTTACTCGCCAGGCGGCACCCCCAGTACGCGGATTAACTCTTAAAAGGTGGTAAAGCTTTGCTCTCCGGCCGGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGtt < 2:204206/151‑1 (MQ=255) ggTTATGCCATGCCAGCGCAAGAGTTGCTTCCCCGCCATGCGGCCCCCCCAGTCCGGGGATTACCTCTTAAAGGGTGGTTAAGCCTTGCTCTGCGGCCGGGCGGCTGTTTAAAAAAGGCGTACGCATCCCCGCCTGTCGTGCGTGCGttt < 1:44109/150‑1 (MQ=255) ggTTATGCCAGGCCAGCGCCAGGGTTGTTTCCTCGCCAGGCGGCACCCCCAGGACGGGGATTAACTCTTCAAAGGTGGTTAAGCCTGCCTCTACGGCCAGGCGGCCGTTTACAAAAGGCGTCCGCATCCCCGCCTGTCGTGCGTGCGttt < 1:4047004/150‑1 (MQ=255) ggTTATGCAAGGCCAGCGCCAGGGAAGATTACTCGCCCTTTGGTATCCCCAGTACGCGGATTAACTCTTTAAAAGAAATAAAGCCTTGCTCTACGGCCAGGCAGACGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc > 1:1515990/1‑151 (MQ=255) atgccatgccAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTACCTCTTAAAAGGTGTTTAAGCCTTGCTCTACGGCCGGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAg < 1:235458/150‑1 (MQ=255) atgccatgccAGCGCCAGAGTTGCTTACTCGCCATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCa < 2:4093477/149‑1 (MQ=255) cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg > 2:1040281/1‑148 (MQ=255) gCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCgg > 1:3416272/1‑151 (MQ=255) cGCAATGCGGCACCCCCAGTACGCGAATTACCTCTTAAAAGGTGGTTAAGCTTTGCTCTACGGCCGGGCGGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATaa < 1:3033999/150‑1 (MQ=255) ccATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAGGGTGGTTACGCTTTGCTCTACGGCCAGGCAGCCTTTTAAAAAAAGCGTACGAATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt < 1:2515989/151‑1 (MQ=255) ccccAGTACGCGGATTAACTCTTCAAAGGTGGTTAACCCTTCCTCTACGGCCGGGCAGCCTTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGa < 2:3245981/149‑1 (MQ=255) cggggATTAACTCTTCAAAGGTGGTTAAGCTTTGCTCTGCGGCCGGGCGGCTGTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGtt < 2:2933481/146‑1 (MQ=255) cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATgg > 2:90336/1‑151 (MQ=255) gaaTTACCTCTTAAAAGGTGGTTAAGCTTTGCTCTCCGGCCGGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTAt < 2:1614202/146‑1 (MQ=255) aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCaa < 2:938653/149‑1 (MQ=255) aaGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt > 2:369840/1‑150 (MQ=255) aaGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTAc < 2:2890691/151‑1 (MQ=255) ttGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGAcc > 2:616347/1‑150 (MQ=255) cGGCCAGGCGGCCGTTTTAAAAAAGCGTACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaa < 1:1768700/147‑1 (MQ=255) ccgggcggcttttttCAAAAGGCGTCCGCAAATCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCg < 2:2773082/138‑1 (MQ=255) ccAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGt < 2:3049495/150‑1 (MQ=255) ggCAGCTGTTTAAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTa < 1:3206253/151‑1 (MQ=255) ggCAGCCTTTTACAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTa < 2:2989353/151‑1 (MQ=255) tttAAAAAGGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTa < 2:2274938/150‑1 (MQ=255) tCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTacac > 1:1706407/1‑151 (MQ=255) | CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACAC > NC_000913/115649‑115946 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |