Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,442,348 | T→G | 11.6% | L355R (CTG→CGG) ‡ | PP_3062 → | tail length determination protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,442,348 | 0 | T | G | 11.6% | 121.8 / 10.7 | 52 | L355R (CTG→CGG) ‡ | PP_3062 | tail length determination protein |
Reads supporting (aligned to +/- strand): ref base T (25/21); new base G (3/3); total (28/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01 |
GCGAACAACTTCAAGAACTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCGGCGATATCTTC > NC_002947/3442206‑3442496 | gcgAACAACTTCAAGAACTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCAGGGCCATg < 2:29197/150‑1 (MQ=255) cgAACAACTTCAAGAACTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGc < 1:411091/150‑1 (MQ=255) aGAACTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCgaga < 2:389782/150‑1 (MQ=255) cTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCAGGGCCATggacg < 1:633744/137‑5 (MQ=255) cTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCAGGGCCATggacg > 2:633744/1‑133 (MQ=255) agaagaTTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACg > 2:157150/1‑132 (MQ=255) agaagaTTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACg < 1:157150/132‑1 (MQ=255) agaagaTTGGCGCCGGTGATATCAAGAAGCCTACTAACGGTGTGTTTGGTGAGGATCGGGCGACACTGAGTACAGCCACTGCTCGAACAGC‑‑GAAGCCCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCgg > 1:1391710/1‑150 (MQ=255) agaTTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGAATATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCggcca < 1:1021261/150‑1 (MQ=255) aagaagGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGAc < 1:1073376/150‑1 (MQ=255) aagaagGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGAc > 2:637242/1‑150 (MQ=255) aagaagGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGAc > 2:637297/1‑150 (MQ=255) ggCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCa > 2:891315/1‑150 (MQ=255) gCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGa < 1:1258445/150‑1 (MQ=255) gggCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACaa < 1:5930/150‑1 (MQ=255) gCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAgg > 1:515617/1‑150 (MQ=255) aTTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCggccaa > 2:112001/1‑103 (MQ=255) aTTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCggccaa < 1:112001/103‑1 (MQ=255) tCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAg < 1:245423/124‑1 (MQ=255) tCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAg > 2:245423/1‑124 (MQ=255) gCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGa > 2:711706/1‑106 (MQ=255) gCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGa < 1:711706/106‑1 (MQ=255) gTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCAATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGaa > 2:904193/1‑150 (MQ=255) cACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAgg < 2:515617/150‑1 (MQ=255) aCTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGAc > 2:1181142/1‑127 (MQ=255) aCTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGAc < 1:1181142/127‑1 (MQ=255) aCTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTGGACCCGGAGAAGGc > 1:648827/1‑150 (MQ=255) cTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGcc > 1:474224/1‑150 (MQ=255) aTACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCGATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGTAGATCAAAGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAgg > 1:129011/1‑150 (MQ=255) tACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAggc > 2:936950/1‑150 (MQ=255) cTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACt > 2:147388/1‑149 (MQ=255) cTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACt < 1:147388/149‑1 (MQ=255) c‑AGAAGGGCATGAACGTCATTGAGGCGTACATGGCCCTGCGCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGaa < 2:474224/150‑1 (MQ=255) ggCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTg < 2:648827/150‑1 (MQ=255) aTGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGgag > 2:719595/1‑150 (MQ=255) gAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGagaa < 1:891315/150‑1 (MQ=255) gTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAg < 2:1202539/99‑1 (MQ=255) gTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAg > 1:1202539/1‑99 (MQ=255) aTTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAg > 1:1133811/1‑100 (MQ=255) aTTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAg < 2:1133811/100‑1 (MQ=255) aTTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTg < 1:719595/150‑1 (MQ=255) aGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCa > 1:517195/1‑150 (MQ=255) aGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAAACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCa > 1:1155632/1‑150 (MQ=255) cGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCg < 2:517195/150‑1 (MQ=255) acgcATGGCCAGGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCgg > 2:431505/4‑150 (MQ=255) tCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGgag > 1:476278/1‑132 (MQ=255) tCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGgag < 2:476278/132‑1 (MQ=255) tCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCGGc > 1:385528/1‑150 (MQ=255) gcATGGCCAGGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGagaaga < 2:1128370/134‑1 (MQ=255) gcATGGCCAGGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGagaaga > 1:1128370/2‑135 (MQ=255) cTGGCCATGCGTTACGTCGAGAAGACCGCCCCGGCCAAGGCCAAGCAGCTCAAGGACGCCGAGGCCAAGATGGACAAGGAAGTCGCCCCGGCGAAGGCCAAGGCGGCCCTGGAAGCGCTGGAGAAAAATGCCTTCACCGGCGATATCTTc > 1:783903/1‑150 (MQ=255) cTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCGGCGATATCTTc > 2:600958/1‑150 (MQ=255) | GCGAACAACTTCAAGAACTGGATGGAGAAGATTGGCGCCGGTGATATCAAGAAGGCTTACAGCGATGTGGGCATTGATTATCAGGCGTCACTGAATACCGGC‑CTGC‑AGAAGGGCATGAACGTCATTGAGGCGTCCATGGCCCTGGCCATGCGTTACGTCGAGAAGACCGACCCGGCCAAGGCCAAGCAGATCAAGGACGCCCAGGCCAAGATCGACAAGGAAGTCGACCCGGAGAAGGCCAAGGCGGCACTGGAAGCGCTGGAGAAGACCCTGCGCACCGGCGATATCTTC > NC_002947/3442206‑3442496 |
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: ‑ |