| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,126 | 0 | A | C | 52.6% | ‑4.9 / 18.8 | 19 | A26A (GCA→GCC) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base A (4/5); new base C (7/3); total (11/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGC > NC_000913/1428886‑1429360 | ncgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGAGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc > 1:45049/2‑250 (MQ=255) tGAAGGCCCCGGCAGTGGCGACGGGCGCACTGACGGTTTCATCGTGGGTTTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTAGGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTCCAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCtn < 1:133480/250‑2 (MQ=255) gCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACCGTGTTATGACAGCCCGCCGGACCCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTACAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGc > 1:197290/1‑251 (MQ=255) cACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGACGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTTAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTa > 1:119838/1‑250 (MQ=255) aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACATCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGTTGAATATGGCAGTAAAAATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGCTGAAAGCAAAACGTAACAGCACCACTGTGGTGGTGAACACCCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTGACAGCATGGACGTTg < 1:212929/251‑1 (MQ=255) gCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAgtgtcatcctgcagg > 1:85562/1‑237 (MQ=255) ccggaaccggtgggCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACTGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGc < 1:148227/171‑1 (MQ=255) gnaaccggtgggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCCCGGTGGTGGGGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGatgtgg > 1:170339/10‑185 (MQ=255) tgGGGCGAATATGGCAGTAAAGATTTCAGGAGTACTGAAAGACGGCACAGGAAAAACGGTACAGAACTGCACCATTCAGATGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATACTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAtt > 1:65686/1‑251 (MQ=255) ggTGANNATNNNAGNNAAGATTTCAGGTGTACTGAAAGANGGCACAGGAAAACCGGTACAGAACTGCACNANCNAGCTNNAAGNAAAACGTAACNGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGCTGAAGNCNGGCGTTACAGCATGGACGTTGAGT‑ACGGTCANTACAGCGTTATTCTGTTGGTGGAANGATTNncgncnnnnnnnnnn > 2:159747/1‑213 (MQ=255) gCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTAc > 1:109053/1‑144 (MQ=255) gtccTGAAAGACGGCACAGGAAAACCGGTACAGAACGGCCCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTgggct < 1:121742/94‑4 (MQ=255) aggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCAACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg > 1:15394/2‑251 (MQ=255) gCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGAGGTGAACACGGAGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa > 1:193195/1‑243 (MQ=255) aaaaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGGAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGt < 1:188498/122‑1 (MQ=255) aaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGGGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATGCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGcc < 1:159747/251‑1 (MQ=255) aGAACTGCACAAACCGGCGGAAAGCAAAAGGTAGCAGCACCACGGGGGTGGTGAACGCGGTGCCCTCAGAAAATCCGGATGAATCCGTCTGTTACAGCATGGACGTTG‑GTGAGGGTAATTATGGCGTCACTGGGTTTGTGGAAGGATTCCCGCCGTCACATGCCGGGGCTATTGCCGTGTATGAAGATTCTCAGCCCGGGACGCTGAATGATTTTCTCGGTCCCAGGACGGAGGATGATGCCCGTCCGGAg < 1:198554/251‑1 (MQ=255) caATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgcgc > 1:169723/1‑251 (MQ=255) gAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTTTCATCCTGCAGGTGGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGc < 1:178861/222‑1 (MQ=255) | CCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGC > NC_000913/1428886‑1429360 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |