| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,615,986 | A→C | 22.6% | D37A (GAC→GCC) ‡ | PP_2288 → | tail fiber protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,615,986 | 0 | A | C | 22.6% | 47.9 / 15.6 | 31 | D37A (GAC→GCC) ‡ | PP_2288 | tail fiber protein |
| Reads supporting (aligned to +/- strand): ref base A (10/14); new base C (4/3); total (14/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.71e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 | |||||||||||
CGGCTTCGGCTGGGGTTTTCTTTTTTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615838‑2616132 | cGGCTTCGGCTGGGGTTTTCTTTTTTTGGTCCCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGac > 1:13248/1‑150 (MQ=255)cGGCTTCGGCTGGGGTTTTCTTTTTTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGac > 2:371370/1‑150 (MQ=255) tCGGCTGGGGTTTTCTTTTTTTAGCTACAGCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAATCTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGt < 1:128908/150‑1 (MQ=255) tttGGAACCCTCTATTTGGCTTACACACCCTACCAAAGTATTCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTCCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTgg < 1:155475/150‑1 (MQ=255) tttGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTgg < 1:371370/150‑1 (MQ=255) ttGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGCGGGCGTACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTggg < 2:13248/150‑1 (MQ=255) ccTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTc < 1:227385/150‑1 (MQ=255) tttGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTaccc < 2:136664/150‑3 (MQ=255) acacCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATc < 2:85433/150‑1 (MQ=255) acacCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATc > 2:334245/1‑150 (MQ=255) cTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACGATCAGGATCGCCTc > 2:189962/1‑150 (MQ=255) cAGCTACGTCCAGTTGAAAGCGGACGGAGCAACAACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCAAGCTGTCCGTGTACACAGTGGTCACGGACTTCACTTGGATAGATTCGTACACCATCAGGATCGCCTCcgc > 2:144965/1‑150 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat > 1:110504‑M2/1‑84 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat < 2:364147‑M2/99‑16 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat < 2:110504‑M2/99‑16 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat > 1:364147‑M2/1‑84 (MQ=255) gacggaGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCg < 1:358470/105‑1 (MQ=255) gacggaGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCg > 2:358470/1‑105 (MQ=255) aCTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCa > 1:240369/1‑150 (MQ=255) cTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCTGTGGACACAGTGGTCACGGACTTAACTTTTGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCAc < 1:189962/150‑1 (MQ=255) ttCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACt > 1:481384/1‑150 (MQ=255) cTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggatat > 1:506046‑M2/1‑57 (MQ=255) cTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggatat < 2:506046‑M2/74‑18 (MQ=255) ttCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTg > 1:80939/1‑54 (MQ=255) ttCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTg < 2:80939/54‑1 (MQ=255) aaCTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACtt < 2:169052/67‑1 (MQ=255) aaCTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACtt > 1:169052/1‑67 (MQ=255) ccTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCa < 2:240369/150‑1 (MQ=255) ggACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGtttc < 1:58602/150‑1 (MQ=255) cAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTcc > 1:13864/1‑150 (MQ=255) ggTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCAc < 2:481384/150‑1 (MQ=255) aagtCCGTGACCACTGTGTCCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCCGGATGGCTCCACg > 1:78125/3‑150 (MQ=255) tGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTTTTTCTAGGTTGTTGACACGCTCCTgg > 1:358452/1‑150 (MQ=255) | CGGCTTCGGCTGGGGTTTTCTTTTTTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615838‑2616132 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |