| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,128 | 0 | A | G | 53.3% | ‑0.9 / 11.9 | 15 | K27R (AAA→AGA) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base A (4/3); new base G (3/5); total (7/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.19e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTT > NC_000913/1428901‑1429366 | gCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCGGGGGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCCCAATCCAGCTGCAAGCAAAACGTACCAGCACCACggtggtgg > 1:87021/1‑251 (MQ=255) cATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGgtggtggtg < 1:63825/251‑1 (MQ=255) gTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCCGCCAACACCGTGTTATGACAGCCCGCCGGTTCAAGCGGGCTTTTTTGTGGGGTGGATATGGCAGTACAGATTTCCGGTGTACTGACAGACAGCACCGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGAAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATgg > 1:80572/1‑250 (MQ=255) aGGTAACACAGTGTGATGACAGCCCGCCGTAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTTCTGAAATACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACATCACCACGGTGGTGGTTAACACTGTGGGCTCAGATAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttg < 1:74078/251‑18 (MQ=255) gCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAg > 1:65632/1‑168 (MQ=255) ggaaaGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACTGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCGTGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCt < 1:81215/249‑1 (MQ=255) naaaGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc > 1:41961/2‑250 (MQ=255) gAAAGACTGCACAGGAAAACCGGTACAGAATTGCACAATCCATGTGAAAGCAAAACGGAACAGCACCACGGTGGTGGTGAACACGCGGGCCTCAGAAAATCCGTATGAAGCCGGGCGTTACAGCAGTGAGGTTGGTTACGGTCAGTACAGCTTTAGTCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGc < 1:122941/249‑1 (MQ=255) gCACAGGAAAACCTGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTTATCCTGCAGGTTGACGGTTTGCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggan < 1:18445/251‑2 (MQ=255) aaaCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTAACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTGTCCACCAGCGCACGCCGGGACCATCACCGTGGATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGc < 1:166676/251‑1 (MQ=255) aaCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGcc < 1:9931/251‑1 (MQ=255) gTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggcc > 1:1070/1‑248 (MQ=255) aGCTGAAAGCCAGACGTAACAGCACCACGGTGGCGGTGCACACGGTGGGCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAACGTGTCATCCTGCAGGCTGACGGCTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcgtctt > 1:150431/1‑237 (MQ=255) ntGAAAGCCAGACGTAACAGCACCACGGTGGTNGNGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATNGATGTGGAGTNCGGNNAGTANANTNTCATCCTGCAGGTTGACGGTTTTCCACCANCGCACGCCGGGACCATCACCGTGNNTGNNNATNNACAAc < 2:121871/177‑1 (MQ=255) cTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCAGGATGAAGACGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGAAGGATGACGGTTTTCCACCATCGCACGACGGGACCATCACCGTGTATGAAGATTCACAAc > 1:121871/1‑178 (MQ=255) | GCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTT > NC_000913/1428901‑1429366 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |