Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,615,990 | A→T | 22.6% | T38T (ACA→ACT) | 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,990 | 0 | A | T | 22.6% | 60.3 / 16.4 | 31 | T38T (ACA→ACT) | PP_2288 | tail fiber protein |
Reads supporting (aligned to +/- strand): ref base A (10/14); new base T (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 = 8.44e-01 |
TCGGCTGGGGTTTTCTTTTTTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615843‑2616132 | 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 > 1:364147‑M2/1‑84 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat < 2:364147‑M2/99‑16 (MQ=255) aGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGACCACTGTGTCcacggacacctggat < 2:110504‑M2/99‑16 (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 > 1:169052/1‑67 (MQ=255) aaCTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACtt < 2:169052/67‑1 (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) | TCGGCTGGGGTTTTCTTTTTTTGGAACCCTCTATTTGGCTTACACACCCTACAGCTACGTCCAGTTGAAAGCGGACGGAGCAACCACAAACTTCCCCTTCAACTTCCCCTACCTGGACACCGCCCATATCCAGGTGTCCGTGGACACAGTGGTCACGGACTTCACTTGGGTAGATTCGTACACCATCAGGATCGCCTCCGCGCCTGTCGCTGGTGGCGTGGTCGAGATTCGCCGGATCACTCCCAAGGATTCCGCCATTGTTTCTTTCCAGGATGGCTCCACGCTCCTGG > NC_002947/2615843‑2616132 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |