Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,452,278 | 0 | A | C | 50.0% | 2.0 / 12.7 | 37 | S18R (AGC→CGC) | PP_5556 | bacteriophage lysis protein |
Reads supporting (aligned to +/- strand): ref base A (12/5); new base C (14/3); total (26/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.88e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACG > NC_002947/3452145‑3452392 | cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc > 1:73374/1‑136 (MQ=255) gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgcg > 2:92400/1‑136 (MQ=255) aCCAAACGTATCAATGGAGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTgg > 1:33283/1‑136 (MQ=255) aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGc > 2:307134/1‑136 (MQ=255) tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTggg > 1:158337/1‑136 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTCGgtg > 1:360445/1‑136 (MQ=255) ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCG‑CAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg < 2:182121/136‑1 (MQ=255) gTCTTGAGCTTCTCAAGCCCCTCTACCAGCGCGGGCTCGTGGTGCTGCTGTGCGGGGTTGGGCCACCCGCCTGGGGGCGGTGGCGCTGGTGCGGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg < 2:198424/136‑1 (MQ=255) gTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg < 1:322479/136‑1 (MQ=255) gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAggacca > 1:198424/1‑132 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCtata > 1:351103/1‑136 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCCGGCCGCCTGGGTGTGGGGGGGCACCGccgaga > 1:343756/1‑132 (MQ=255) cAAGGCCCTCTACAAGAGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGAGCTGCTGCTGGCCGCCAGCGCAGGCGCCGCCCCGGCACCCAGGGTGTGTCACGCCAACGCCaata > 1:174294/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATaa > 1:141144/1‑136 (MQ=255) ccccTCTACACGCGAGCGCTCGTGGGGCTGCAGTGAGCGCTGGGGCCCCCCGCCTGGGCGCCGGGGCGCGGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc < 1:109950/135‑1 (MQ=255) aaGCGAGCGCTCGAGGTGCAGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCACCGCCGATAAGGCGCCgcgcgc > 2:291594/1‑136 (MQ=255) agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATc > 2:22118/1‑136 (MQ=255) agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAATGCGCAGCTCGCTGATc > 2:183028/1‑136 (MQ=255) agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGTGTGGCCCGCCAACCCCAAAAAGGCGCAgcgcgcaggta > 2:278414/1‑131 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGCCCCCCCCCGCGCCGCCGGGGTGTGGCCGGCCCCCCCCGAGAAGGCGCCGCTCGCTGGTCAgg > 1:82374/1‑136 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGCGCCGCCGGGGTGTGGGAGGCCAACGCCAAAAAGACGCAGCTCGTTGATCAgg > 2:355513/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAAGGAGCAGCTCGCGGATCAGGcc > 1:184953/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCCACGCCAAAAAGGCGCAGCAGGCGGATCACGcc > 2:345464/1‑136 (MQ=255) tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCAGGCGGCGCCCGGGCAGCCTGCGTGTTGCCGGCCCACGCCGAAAAGGCGCCGCGCGCGGCTCAgggc > 2:255266/1‑134 (MQ=255) gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGGGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCAACGCCGATAAGGAGCCGCTCGCTGATCAGGGCGCTGACATCGGCAAGCa > 2:357101/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCAACGCCTATAAGGGGCCGCGCGCTGAGCCGGCCGATGACTTAGGCAAGCAGt > 2:75245/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCCGCTCGCGGATCAGGCCGATGACTTCGGGAAGCAGt > 2:89877/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCACGCGGCTCAGGCCGATGACTTCGGCAAGCAGt > 2:261179/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTAAAAAGCGCAGCTCGCAGATAAGGCCGATGACTTCGGCAAGCAGt > 1:176912/1‑136 (MQ=255) gcTAGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCGAGCGCCGCCCCGGCAGCCTCGGAGAGGCAGGCCCACCCCTATAACACGCAGCTCGCAGATCAGCCCCAAGACTTAAGCAAGCAga > 1:40474/1‑135 (MQ=255) ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGttt > 1:69668/1‑135 (MQ=255) gggCCCCCCGGCTGGGGGCGGGGGCGCTGCTGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc < 1:28521/136‑1 (MQ=255) ccgcGCCTGGCCGCCGTGGCGGGGGTGCTGGCCCACAGCGCGGGCGGCGCCCGGGCAGTCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga < 2:69488/133‑1 (MQ=255) ccccGCCGGGTCGCCGGGGCGGTGGTGCGGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga < 2:252702/135‑1 (MQ=255) gtggcggTGGCGCTGCTGCGGGCCGCTTGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGGAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg < 1:1039/130‑1 (MQ=255) gcGCCGGGGCGCTGGTGCGGGCGGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 1:37512/135‑1 (MQ=255) gcGCCGGGGCGCGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 1:89891/135‑1 (MQ=255) ggggCGCGGGTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg < 1:22118/134‑1 (MQ=255) | CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACG > NC_002947/3452145‑3452392 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |