| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 581,572 | 0 | T | C | 43.8% | ‑0.7 / 13.3 | 16 | N13N (AAT→AAC) | tfaD | putative tail fiber assembly protein TfaD |
| Reads supporting (aligned to +/- strand): ref base T (4/5); new base C (2/5); total (6/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACC > NC_000913/581337‑581806 | aaTACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAATGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGc < 1:167903/249‑1 (MQ=255) cGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCGGCGGCATTTTGGCCGCGCCGGGCTTCGCTCACTGGTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCAGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCGGGGAAACACTGCCCTGTCAGCGGGGCATCAGGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt < 1:94806/251‑1 (MQ=255) gggTTGCTGAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttn < 1:103057/251‑3 (MQ=255) gAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAATGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGTTAACGATTACATACGTGAGGTCAACGTTGTTAAGTCTGCCCgtgttgtctattc < 1:120634/250‑10 (MQ=255) gAGTGAATATATCGAACAGTCAGGTTAAGGGGCTGCAGCACTCTGCGCGCGCGGGTCTACTTTCGCTGTACCGACCTGAGCCGCAAACGGCGGTTGAAGGGGCAGACGATAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATTCGTTAGGTCAACGTTTTTAAGTCTGCCCGTGTTGGCTAtt < 1:158885/249‑1 (MQ=255) ggCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAAGCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGGCATCATGAATGCGATGTGCTGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTTCACGTGTCGTTTATTCCAAAAGGCTGCTGGGTGTTTATGc < 1:145066/202‑1 (MQ=255) ggcggTTGAATGGGCAGACGATAATTNCTATCTNNCCNAAGAGTCGGCCTACCAGGANGGNNNCNNNGNCACGTNGCNNTTTCAACGCNCGATCATGAATGCGATGGGTAACGATTACCTACGTGAGGTCAACGTTGTTAAGNCTGCCCGTGNNGGCTATTCAAAAATGTTGCTGNGTGNNTNNgngtangtn > 2:42699/5‑185 (MQ=255) nttGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGGGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAGTCCAGCTGGCTTCGTGGCTGTTTTCAACAGtg > 1:108569/2‑250 (MQ=255) gatAATTACTATCTCCCCAAAGAGTCGGCCTACCAGGAAGGGCGCTGGGAAACGTTGCCGTTTCAACGCGCGATCATGAATGCGATGGGTAACGATTACATACGTTAGGTCAACGTTGTTAAGTCTGCCCGTGNTGGNNNNNc > 2:110848/3‑143 (MQ=255) nnnnaCCNGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc < 2:104174/137‑1 (MQ=255) ncATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc > 1:104174/2‑141 (MQ=255) ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg > 1:48689/1‑96 (MQ=255) gatcatGAATGCGATGGGTANCGANTACANACNTGAGGTCAACGTTGTNAAGTCTGCCCGTGTTGGCTATTCAAAAATGTTGCTGGGTGTTTANNCNNNTNTtattcag < 2:120599/108‑5 (MQ=255) gATGGGCANCGAANANNNCCGTGNGGTGAATGTNGTGAAGTCTGCCCGTGTCGNt > 2:19057/1‑55 (MQ=255) gATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCggn < 1:19057/55‑2 (MQ=255) aCATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAcc < 1:23684/251‑1 (MQ=255) | GATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACC > NC_000913/581337‑581806 |
| 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: ‑ |