Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,601,862:1 | +C | 10.0% | coding (1671/5418 nt) | 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,601,862 | 1 | . | C | 10.0% | 142.9 / 11.6 | 40 | coding (1671/5418 nt) | PP_4924 | subtilase family serine protease |
Reads supporting (aligned to +/- strand): ref base . (19/17); new base C (2/2); total (21/19) | |||||||||||
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.74e-01 |
GCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACC > NC_002947/5601713‑5602006 | gCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCg > 1:769452/1‑150 (MQ=255) gtcgCTGGCACTGATCACGAATACATCGTCGCCCGCCCGGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTgc > 2:460990/1‑150 (MQ=255) gCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTcc > 2:616918/1‑150 (MQ=255) tGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCg > 2:661058/1‑150 (MQ=255) gATAACGAATACATCGTCTCCCGCGCCGTTTTTATAGGTGTTTCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGt < 2:769452/150‑1 (MQ=255) cccgccccgcCGTGATAGGTGTTGCCGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGGCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTgg < 1:275049/150‑1 (MQ=255) ccccgccGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGAGCCGCGGTCAAGACAGCATAGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGa > 1:639164/1‑150 (MQ=255) cgccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGCAGCGCtcg > 1:45459/1‑111 (MQ=255) cgccgTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCGCAGCGCtcg < 2:45459/113‑3 (MQ=255) aTAGGTGTGGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCt > 2:213164/1‑150 (MQ=255) aGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTc > 1:171795/1‑150 (MQ=255) gctgcCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGt > 2:764101/1‑150 (MQ=255) ccgccaccAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACGCTGCCTGATGTAAGCCCTTCGGGAAGTGGGCCCg > 2:103355/1‑150 (MQ=255) accaaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGa < 1:616918/150‑1 (MQ=255) aaccaGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGa > 1:719736/1‑150 (MQ=255) caGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGt < 2:719736/150‑1 (MQ=255) caGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGt < 2:639164/150‑1 (MQ=255) ccAACTGCCCTCTGTGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGGCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATCTATGCCCTTCGGGAAGGGGGCCCGCTTGACGAAGTCAg > 2:341407/1‑150 (MQ=255) aaCTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTa < 1:661058/150‑1 (MQ=255) cTGCCCTCTGGGGCCGCGGTTAAGACATCATCGCCCGCGCTGGCGGAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGATGATGTAGGTTCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTaga < 2:696784/150‑1 (MQ=255) cccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTg < 2:727529/150‑1 (MQ=255) cccTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTg > 1:602541/1‑150 (MQ=255) gcgGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTgg < 2:602541/150‑1 (MQ=255) ggTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGc < 1:460990/150‑1 (MQ=255) tCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGt < 1:368339/133‑1 (MQ=255) tCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGt > 2:368339/1‑133 (MQ=255) tCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCt < 1:683319/150‑1 (MQ=255) cAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTc > 2:480309/1‑150 (MQ=255) gACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGAt > 1:466402/1‑150 (MQ=255) aCATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATc > 1:122120/1‑150 (MQ=255) tcGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTa > 1:165307/1‑150 (MQ=255) cccGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCg > 2:741715/1‑150 (MQ=255) ccGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGa < 1:213164/150‑1 (MQ=255) gcgcTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACg < 2:466402/150‑1 (MQ=255) gTAAAGCGTCTGTACGCCG‑AGCGCTTCGGCTTCCAGTGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGc < 1:341407/150‑1 (MQ=255) gTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCACCGTGCTGCCTTGAGCACTCTGGCTGCTGTTGGCCTTTCGGAAAGTGGGCCTGCTTGACGAAGTCCGTAGCGAGACACCCTCCTGGGGCTCGATCCGATTATCGCCGGTATCggg > 2:267620/1‑149 (MQ=255) ggaCGCCGCAGCGCT‑CGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGa > 2:592602/3‑122 (MQ=255) ggaCGCCGCAGCGCT‑CGGCGTACAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGa < 1:592602/120‑1 (MQ=255) tACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGTCCCTTCGGGAAGTGGGCCCGCTTGACGaa > 1:300252/1‑85 (MQ=255) tACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGaa < 2:300252/85‑1 (MQ=255) aCGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCAcc < 1:480309/150‑1 (MQ=255) | GCCCAATCGTCGCTGGCACTGATCACGAATACATCGTCGCCCGCCCCGCCGTGATAGGTGTTGCTGCCACCGCCACCAACCAGCCAACTGCCCTCTGGGGCCGCGGTCAAGACATCATCGCCCGCGCTGCCGTAAAGCGTCTGTACGCCG‑AGCGCTGCGGCGTCCAGCGTGCTGCCTTGAGCACTCTGGCTGATGTAGGCCCTTCGGGAAGTGGGCCCGCTTGACGAAGTCAGTAGAGTGCCACCCTCCTGGGCCTCGATCCGATTACCGACGGTATCGGCGAGAAACCTCACC > NC_002947/5601713‑5602006 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |