Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,603,056 | C→T | 87.3% | W142* (TGG→TGA) | PP_4924 ← | subtilase family serine protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,603,056 | 0 | C | T | 87.3% | 156.5 / 13.1 | 63 | W142* (TGG→TGA) | PP_4924 | subtilase family serine protease |
Reads supporting (aligned to +/- strand): ref base C (4/4); new base T (26/29); total (30/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01 |
GTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAGG > NC_002947/5602911‑5603196 | gTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCag > 5:303307/1‑149 (MQ=255) aCTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCagta < 4:287166/148‑1 (MQ=255) aCTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCagta < 4:145136/148‑1 (MQ=255) cTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCagtagt > 2:303320/1‑149 (MQ=255) ccTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCtctc > 8:268813/1‑149 (MQ=255) tCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTTATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTc < 3:296347/148‑1 (MQ=255) ccACTCCTTAGCAATTAAACGGTTTTCGCCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGa > 1:210950/1‑149 (MQ=255) aCTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGa > 2:226099/1‑147 (MQ=255) aCTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACa < 7:81957/149‑1 (MQ=255) ttAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAAt > 7:67331/1‑149 (MQ=255) aGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAAt < 6:71501/147‑1 (MQ=255) gCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATcac < 2:113488/149‑1 (MQ=255) cggAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAAc < 2:14166/149‑1 (MQ=255) ggAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACg > 4:303798/1‑149 (MQ=255) agaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTc > 2:213299/1‑148 (MQ=255) agaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTcc < 1:303320/149‑1 (MQ=255) gaaACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCa > 3:156303/1‑149 (MQ=255) aaCTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATg > 5:11721/1‑149 (MQ=255) cgcatGATGGCAATGGCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCGGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGcc > 1:237466/5‑149 (MQ=255) tactGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGccc > 2:335551/1‑149 (MQ=255) aTGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCctct > 4:189208/1‑147 (MQ=255) gcaatggcatcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGa < 2:237466/142‑1 (MQ=255) ccattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATCAAAACACCGGACTCAACGTCCATGCCCTCTCAgg < 1:226099/149‑1 (MQ=255) ccattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAgg < 4:156303/148‑1 (MQ=255) cattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCa < 7:262487/145‑1 (MQ=255) cattgccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACAGGACTCAAAGTCCATGCACTCTCAgag > 8:262487/1‑146 (MQ=255) ccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc < 6:11721/148‑1 (MQ=255) ccatCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACAAGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc > 4:5439/1‑148 (MQ=255) catcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTc > 5:326818/1‑147 (MQ=255) catcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCaa < 6:303307/149‑1 (MQ=255) atcatGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAATGTCCATGCCCTCTCAGGATATTCAAg < 2:210950/149‑1 (MQ=255) cGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc < 7:268813/149‑1 (MQ=255) ggTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc > 3:17298/1‑148 (MQ=255) aGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGaaa < 4:17298/148‑1 (MQ=255) tgAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTc > 1:362213/1‑149 (MQ=255) aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑aaa < 2:326996/93‑1 (MQ=255) aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑aaa > 1:326996/1‑93 (MQ=255) aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc < 6:326818/148‑1 (MQ=255) aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc < 3:303798/148‑1 (MQ=255) aaaggaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTcc < 3:189208/148‑1 (MQ=255) gaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCtt < 1:211040/146‑1 (MQ=255) gaaagTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGAATCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCtt > 2:211040/1‑146 (MQ=255) tAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAg > 5:183025/1‑69 (MQ=255) tAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAg < 6:183025/69‑1 (MQ=255) tCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGAc < 1:213299/149‑1 (MQ=255) tCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACctc < 8:67331/149‑1 (MQ=255) gTTAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGAAATCACATC‑AAACACAGGACTCAACGTCCATGCCCTCTCAGGATATTAAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt > 7:41806/1‑149 (MQ=255) ttAGATTCCCCGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCtttt > 2:14296/1‑149 (MQ=255) ccccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGc > 2:9115/1‑146 (MQ=255) ccccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGc < 1:9115/146‑1 (MQ=255) cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc > 6:253902/1‑107 (MQ=255) cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTc < 5:253902/107‑1 (MQ=255) cccGCTGAGGAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGACCTCTCAGCATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACATCTTTTTGGCTTAc > 5:159501/1‑149 (MQ=255) aGGAAAGCACACGCCGTCAGTAGTCCTCTCGATTTTCGACATCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGAc < 3:5439/147‑1 (MQ=255) ggAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt > 6:219837/1‑132 (MQ=255) ggAAAGCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCttt < 5:219837/132‑1 (MQ=255) gCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGcc > 7:208165/1‑71 (MQ=255) gCACATGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGcc < 8:208165/71‑1 (MQ=255) cacaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg < 1:335551/149‑1 (MQ=255) acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACTCCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg < 5:282697/148‑1 (MQ=255) acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg > 6:370510/1‑148 (MQ=255) acaTGCCGTCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAg > 6:282697/1‑148 (MQ=255) acaTGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAgg < 2:362213/149‑1 (MQ=255) | GTGACTCGACCCTGGATATCCACTCCTTAGCAATTAAACGGTTTTCGGCGGAAGAAACTACTGATGCCATTGCCATCATGCGGTGAGTGTGAAAGGAAAGTGTAAATTCATCCGTTAGATTCCCCGCTGAGGAAAGCACATGCCGCCAGTAGTCCTCTCGATTTTCGACAGCAATCACATC‑AAACACCGGACTCAACGTCCATGCCCTCTCAGGATATTCAAGGCTTTCGAAAGTTTCCTTATGAAATGCCGCGACCTCTTTTTGGCTTACCTGACCGACAATAGG > NC_002947/5602911‑5603196 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |