Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 296,504 | 0 | T | G | 54.5% | ‑6.5 / 13.9 | 11 | N198T (AAC→ACC) | intF | CP4‑6 prophage; putative phage integrase |
Reads supporting (aligned to +/- strand): ref base T (1/4); new base G (6/0); total (7/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAAC‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTC > NC_000913/296442‑296567 | tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTg > 1:2369168/1‑65 (MQ=255) ttGCCTTTTTGTTCTCCGCCCCGCTGCTATGTTATCTGGGGTTAGCTCCGATAACGGCAGGTTGa < 1:1620990/65‑1 (MQ=255) ttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa < 1:2997054/65‑1 (MQ=255) ttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa < 1:2448386/65‑1 (MQ=255) cTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCaa < 1:3310124/65‑1 (MQ=255) tCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGTGGAAAC‑t > 1:3194347/1‑65 (MQ=255) aCTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCCAAC‑TAGCCa > 1:162814/1‑65 (MQ=255) aCTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGGGCAAAA‑CAGCCa > 1:2135673/1‑65 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGTGGAAAC‑AAGCaaaa > 1:368639/1‑64 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGTGCGAAA‑CAGcccac > 1:1666331/1‑62 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑AAGCCAAt > 1:3020337/1‑65 (MQ=255) gCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGAGCACAC‑AAGCGGATggg > 1:3349559/1‑64 (MQ=255) gctgctATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAA‑TAGCAAATGGGCCAGc > 1:2543946/1‑65 (MQ=255) gctgctATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑TAGCCAATGGTCCAGc > 1:3305337/1‑65 (MQ=255) tATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAA‑TAGCCAAAAGGGCAGCCGaaa > 1:699768/1‑64 (MQ=255) tATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGAAAC‑TAGCCAATGGTCCAGGCGAAg > 1:3195999/1‑65 (MQ=255) gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑AAACCCATGGGCCAGCCGGAGGCg > 1:933097/1‑65 (MQ=255) gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAA‑TAGCCAATAGTCCCGCCGAAggtg > 1:1101359/1‑62 (MQ=255) gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGAAAA‑CAGCCCATGGGCCCGCCGAAggcc > 1:188668/1‑62 (MQ=255) tCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAAA‑CAGACCATGGGCCCGCCGCAGGCGGGcg > 1:309819/1‑64 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGTGCAAAC‑TAGCCAATGGGCCAGCCGAAGTCGGGccc > 1:871999/1‑64 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAC‑TAGCCAATGGGCCAGCCGGAGGCggggcc > 1:594784/1‑62 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGGGCCAAA‑TAGCCAATGGGCCAGCCGAAGGCGGGCCt > 1:405834/1‑65 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑TAGCCAATGGGCCCGCCGAAGGCgggggc > 1:3165821/1‑62 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑TAGCCAAAGGTCCCGCCGCAGGAGGGCCt > 1:2720241/1‑65 (MQ=255) ggTTAGCTCCGATAACGGCAGGGTGTGCAAAT‑TAACCAACCGGCCACCCGGAGTCCGGGCTCGGc > 1:1521315/1‑65 (MQ=255) ggTTAGCTCCGATAACGGCAGGGTGGGCAAAC‑TAGCCAATGGGCCAGCCGAAGGCGGGCCTTGGc > 1:25423/1‑65 (MQ=255) gCTCCGATAACGGCAGGGTGGGGAAAC‑CAGCCCAAGGGGCAGCCGCAGGCGGGGCTTGGGCTCtt > 1:1133598/1‑65 (MQ=255) gCTCCGAAACGGGGAGG‑TGAACAAACTTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCtt < 1:1917674/65‑1 (MQ=255) ggTTGAGCAATT‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCt < 1:819593/65‑1 (MQ=255) ggTTGAGAAACT‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCt < 1:1020873/65‑1 (MQ=255) ggTTGACAAACC‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCt < 1:1068471/65‑1 (MQ=255) gTTGAGCAAAC‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTc > 1:2167976/1‑65 (MQ=255) | TTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAAC‑TAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTC > NC_000913/296442‑296567 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |