| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,429,246 | 0 | A | T | 52.6% | ‑5.8 / 19.2 | 20 | G66G (GGA→GGT) | stfR | putative prophage side tail fiber protein StfR |
| Reads supporting (aligned to +/- strand): ref base A (5/4); new base T (7/3); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCC > NC_000913/1429039‑1429441 | tgGGGCGAATATGGCAGTAAAGATTTCAGGAGTACTGAAAGACGGCACAGGAAAAACGGTACAGAACTGCACCATTCAGATGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATACTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGAtt > 1:65686/1‑251 (MQ=255) ggTGANNATNNNAGNNAAGATTTCAGGTGTACTGAAAGANGGCACAGGAAAACCGGTACAGAACTGCACNANCNAGCTNNAAGNAAAACGTAACNGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGCTGAAGNCNGGCGTTACAGCATGGACGTTGAGT‑ACGGTCANTACAGCGTTATTCTGTTGGTGGAANGATTNncgncnnnnnnnnnn > 2:159747/1‑213 (MQ=255) aggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCAACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg > 1:15394/2‑251 (MQ=255) gCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGAGGTGAACACGGAGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa > 1:193195/1‑243 (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) cacGTTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCAGGGAGGTGGAGT‑ACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGCTTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACGGATGGTGGAAGAGGTGGTGCGTAACGCGTCCGCgg < 1:172378/251‑1 (MQ=255) gcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTAAAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCAATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTATTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGt > 1:150114/2‑237 (MQ=255) cctaaGAAAATCCGGATGAAGCCGGGCATTCCAGCATGGACGTTGAGT‑ACGGTCGGTACAGCGGTATTCTGTTGGTGGATGGTTTCCCGCCGTCACATGCCGTGACCACTACCGTGTATGAAGAGTCCCAACCCGGTACGCTGAATGATTGTCTCGGTGACATGACGGAGGATGATGCCCGTCCGGAGGCCCTGCGCCGTTTTGAACTGATGGTGGACGAGGTGGGGCTTAACGCGTGCGCGGTGGCACAg < 1:16949/247‑1 (MQ=255) gaaTCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGGGGTGGCACAGagtacgg > 1:93871/2‑245 (MQ=255) ggATGAAGCCGGGCGTTACAGCATGGATGTGGAGT‑ACGGCCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCCCCATCGCACGCCGGGACCATCACCGTGTATGCAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGCAGAGGTGGAGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAgacg > 1:29840/1‑248 (MQ=255) gggCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTANNGTGTAtnn > 2:17771/1‑96 (MQ=255) gggCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTg < 1:17771/222‑1 (MQ=255) gTTACAGCATGGCTGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGCTGACGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc > 1:218268/1‑251 (MQ=255) gTTACAGCATGGATGTGGAGT‑ACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGATGTGCTGCGTCGTCGTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAGcc < 1:73144/251‑1 (MQ=255) gTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGcc > 1:99726/1‑251 (MQ=255) gtggAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAttt > 1:5175/1‑77 (MQ=255) gtggAAGGATCCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAttt < 2:5175/77‑1 (MQ=255) | TGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGT‑ACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAGAACACGGCAGCCGCGAAGAAGTCAGCC > NC_000913/1429039‑1429441 |
| 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: ‑ |