| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,164 | 0 | C | G | 46.7% | ‑1.5 / 18.1 | 15 | A39G (GCC→GGC) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base C (4/4); new base G (2/5); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.08e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCAC > NC_000913/1428916‑1429412 | tATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACTTAACATCACCACGGTGGTGGTGAACACGGTGGnc < 1:181883/250‑1 (MQ=255) gTGGCGACTGGCGTACTGACGGAGTCTTCGTTGGGGTCGGTGATAAATTCTTATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTTGGCTTTGGTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCGCAGGAAAACCGGTAC‑AGACCTGCACCATTCAGCGGAAAGCCAGACGTAACAGCGCCACGGTGGTGGTGAACACGGTGGGCTCAg < 1:213237/238‑1 (MQ=255) ggCGCACTGACGGTTTCATCGGGGGTGTCGGTTATAAGTTCTGATGAGCCAGGTA‑CACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACTTAACAGCACCACGGTGGGGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGc < 1:182338/249‑1 (MQ=255) ggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTTCAGCATg < 1:185238/235‑1 (MQ=255) gTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAgn < 1:59139/208‑2 (MQ=255) gtgtTATGACAGCCCGCCGGTTCAGGCGGGATTTTTTGGGGGGGGGATATGGCAGTACAGATTTCAGGGGCACTGAAAGAAGGCACAGGAAAACCGGGACCAAAAC‑GCACAATCCAGATGACAGCAAAACCGACCAGCACAACGGTGGTGGTGAACACGATGCCCTCAGAAAATCCGGATGACGCCGGGCGTTACAGCCTGGACGTTGAGCGCGGTAAGTACAGCGCTATTCTGTTGAGGGACGGATTccc > 1:100541/1‑251 (MQ=255) gggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCGAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGTAGGACGGTCAGTACAgtgtcat < 1:194167/203‑7 (MQ=255) tttttttG‑GGGGTGAATATGGCATGAAAGATTTCAGGTGTAATAAAAGACGGCACAGGAAA‑CCCGTACCAAAACTGCAAAAGCCAGTTGAAAGAAAAAAGCAACAGCACCAGGGTGGTGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTCCAGCTTGGACGTTGAGTACTGTCAGTACATGTTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTa < 1:7013/249‑1 (MQ=255) gACGGCACAGGAAAACCGGTAC‑AGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTTAACACGGTGGGTTCAGACAATCCGGATGAAGCCTGGCGTTACAGCATGGa < 1:84132/125‑1 (MQ=255) caccattcAGCTGAAAGCCAGACGTCNCAGCACNNCGNTGGTGGGGAACACGGTGGGNTCNNNGNNNCNGGATGNAGNNGGGCGTTACNGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTNCACCATCGCNNGCCGGGACCATCACCGTGTATNAAGNNTNNCNACCNggn > 2:54309/8‑192 (MQ=255) caccattcAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGc < 1:54309/227‑1 (MQ=255) gCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCGGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTCCGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGGGCCATGACGGAGGGTGATGCCCGTCCGGAGGCACTGCGCAGTTTTg > 1:162786/1‑251 (MQ=255) aaGCAAAACGTCACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGa > 1:82048/1‑49 (MQ=255) gTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGa > 1:73356/1‑250 (MQ=255) gtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGAAGTACGGCCACTCCAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCCCGCCGGGACCATCACCGTGACTGAAGATTCACAACTGGGGACGCTGCATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGGGGAAGAGGTGGCGAGTAACGCGTc > 1:114713/1‑251 (MQ=255) gAACANNATNNNCTNNGAAAATCCGGATGAAGCCGGGCGNTACAGCATGGACGTTGAGTACGGTCAGTANANCNTTATNNTGTNGGTGGAAGGANTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACNCNGTACGCTGAATGATTTTCTCGGTGCCATGACNGAGGATGATGCCCGTCCGGAGGCANTGCGNNggntnnnnnnn > 2:166891/1‑210 (MQ=255) gAACACGCTGACCTCAGAAAAGCCGGATGAAGCCGGGCGTTACAGCAGGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATGACCGTGTATGAAGATTCTCAACCCGGTACGCGGAATGATTTTCTCGTTGCCATGACGGAGGATGATGCCCGGCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcgc < 1:166891/246‑1 (MQ=255) nccTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGAGCGTAACGCGTCCGCGGTGGcac > 1:40209/2‑250 (MQ=255) | TATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAC‑AGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCAC > NC_000913/1428916‑1429412 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |