Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 296,504 | 0 | T | G | 58.3% | ‑2.5 / 14.2 | 12 | N198T (AAC→ACC) | intF | CP4‑6 prophage; putative phage integrase |
Reads supporting (aligned to +/- strand): ref base T (2/3); new base G (7/0); total (9/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTACTTTCGCCTC > NC_000913/296442‑296567 | tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTg > 1:359582/1‑65 (MQ=255) tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTg > 1:3624518/1‑65 (MQ=255) tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTg > 1:3698404/1‑65 (MQ=255) tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTg > 1:3404182/1‑65 (MQ=255) ttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa < 1:3570596/65‑1 (MQ=255) ttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa < 1:2961861/65‑1 (MQ=255) cagtCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCaa < 1:1360557/62‑1 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAACA‑GCCAAt > 1:2340977/1‑65 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAATA‑GACaaa > 1:2316495/1‑64 (MQ=255) tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGAGCAAACAA‑GACaaa > 1:1822098/1‑64 (MQ=255) gCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGAAACCA‑GCCCATggg > 1:1006618/1‑64 (MQ=255) tATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGAGGAAAATA‑GCCCAAGGGCCAGCCGGAg > 1:3099390/1‑65 (MQ=255) gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATA‑GCCCATGGGGCAGCCGAAGGC‑g > 1:290322/1‑65 (MQ=255) gTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGGGCAAAATA‑GCCCAAGGGCCCACCGAAgggc > 1:2324086/1‑62 (MQ=255) aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATA‑GCCAATGGGCCAGCCGCAGGC‑ggg > 1:2736008/1‑65 (MQ=255) aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAACA‑GCCAACGGGCCAGCCGAAGGC‑ggg > 1:468587/1‑65 (MQ=255) aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGTC‑ggg > 1:3585597/1‑65 (MQ=255) aaTCTGGGGTTAGCTCCGATAACGGCAGGGGGGGGAAACTACGCCAACGG‑CCCGCCGCAGTC‑ggg > 1:249686/1‑65 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTA‑GCCCCTGGTCCCGCCGACGAC‑gggggc > 1:2694855/1‑62 (MQ=255) cTGGGGTTAGCTCCGATAACGGCAGGGGGGGCAAACAA‑GCCCATGGGCCAGCCCGAGGC‑GGGccc > 1:777159/1‑64 (MQ=255) ggggTTAGCTCCGATAACGGCAGGGTGAGCAAAATA‑ACCAATGGTCCAGCCGCAGTC‑GGGCCCTg > 1:497763/1‑65 (MQ=255) ggTTAGCTCCGATAACGGCAGGGTGAGGAAACTA‑GCCAATGGGTCAGCCGAAGGCGGGGCCTTgc > 1:2251893/1‑64 (MQ=255) ggTTAGCTCCGATAACGGCAGGGTGAGGAAACTA‑GCCAATGGGGCAGCCGAAGAC‑GGGCCTTGGc > 1:3182860/1‑65 (MQ=255) ggTTAGCTCCGATAACGGCAGGGTGAGCAAAAAA‑GACAATGGGCCAGCCCAAGGG‑GGGGCCTGGc > 1:3040055/1‑65 (MQ=255) gTTAGCTCCGATAACGGCAGGGTGGGGAAACTA‑GCCAAAGGGCCAGCCGGAGGC‑GGGCCCTgggc > 1:2342641/1‑63 (MQ=255) gCTCCGATAACGGCAGGGTGAGCAAACTA‑GCCAATGGGCCCGCCGCAGCC‑GGGGCTCGGCCTCtt > 1:1080860/1‑65 (MQ=255) ccTCCGATAGCGCCAGGTTGAGAAAACTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCtt < 1:996016/64‑1 (MQ=255) tCCGATAACGGCAGGGTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGAC‑GGGGCTTGGCCTCtttt > 1:347929/1‑65 (MQ=255) tCCGATAACGGCAGGGGGGGGAAACAA‑ACCAAAGGGCCAGCCGAAGGC‑GGGCCGTGGCCTCtttt > 1:1412467/1‑65 (MQ=255) taccggcggTTTGAGCAATTTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTAc < 1:2484015/62‑1 (MQ=255) taccGGCAGGTTGAGCAACCTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTAc < 1:10081/62‑1 (MQ=255) gTTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTACTTTCGCCTc > 1:3461317/1‑65 (MQ=255) gTTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTACTTTCGCCTc > 1:2051024/1‑65 (MQ=255) | TTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTA‑GCCAATGGTCCAGCCGAAGTC‑GGGCCTTGGCCTCTTTTTTTACTTTCGCCTC > NC_000913/296442‑296567 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |