Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,615,995 | T→C | 25.0% | V40A (GTC→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,995 | 0 | T | C | 25.0% | 41.0 / 13.1 | 28 | V40A (GTC→GCC) ‡ | PP_2288 | tail fiber protein |
Reads supporting (aligned to +/- strand): ref base T (9/12); new base C (4/3); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01 |
TTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615862‑2616132 | 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:334245/1‑150 (MQ=255) acacCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATc < 2:85433/150‑1 (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:110504‑M2/99‑16 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat < 2:364147‑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) 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) tCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGa > 1:247241/1‑121 (MQ=255) tCACGGACTTCACGTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGa < 2:247241/121‑1 (MQ=255) | TTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615862‑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: ‑ |