Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,600,638 | C→A | 22.3% | G965G (GGG→GGT) | 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,600,638 | 0 | C | A | 22.3% | 27.7 / 2.8 | 18 | G965G (GGG→GGT) | PP_4924 | subtilase family serine protease |
Reads supporting (aligned to +/- strand): ref base C (9/5); new base A (2/2); total (11/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01 |
CTCCTTGATCATTGCGCGCGCCAACCATGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATACCACTGCCGCCCCGAC > NC_002947/5600502‑5600773 | ctcctTGATCATTGCGCGCGCCAACCATGACCCCAGCCACGTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTCCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGaaaa > 1:588286/1‑148 (MQ=255) cAACCATGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTCAGATCGGGGTGCCCCATGTCGAAGATTTCCGGCGCCACGGCAAAct > 1:297426/1‑150 (MQ=255) aCCATGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGACGATTGTCGAAGATTTCCGGCGCCACGGCAAActcc < 1:511696/150‑1 (MQ=255) ccATGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTACCCAATGTCGAAGATTTCCGGCGCCACGGCAAActccc > 2:511696/1‑150 (MQ=255) aTGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAActcccct > 1:525904/1‑150 (MQ=255) cccAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGACCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCgg > 1:384283/1‑150 (MQ=255) gAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGcc < 2:525904/150‑1 (MQ=255) aTTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCAc > 1:464346/1‑150 (MQ=255) agCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGACCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGt < 2:384283/150‑1 (MQ=255) cGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATg < 2:519800/150‑1 (MQ=255) gTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATCCCGGcgcgt < 2:154839/150‑2 (MQ=255) ccAGATGGGATCGATATTCGGCTTGAGATCGGGGTGAAGCATGCCACAAAAAAGCGGCTCCACGGCAAACTTCCTTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGGCATGCCGTACCGGCGCGACATTGATCTCGTc > 1:261883/1‑150 (MQ=255) aTGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGa < 2:430836/150‑1 (MQ=255) gATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATAcc > 2:333398/1‑150 (MQ=255) aTCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGTCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATACCa > 1:352506/1‑150 (MQ=255) aTCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCATATACCa > 2:324039/1‑150 (MQ=255) cGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATACCACt < 2:464346/150‑1 (MQ=255) ttGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATACCACTGCCGCCCCGAc > 2:440158/1‑150 (MQ=255) | CTCCTTGATCATTGCGCGCGCCAACCATGACCCCAGCCACTTGAGTGGCATGATTGGACGCCAGAGCGGGCAGCGTGCCATTGGTGCGTTGGGTTTGCAGCCAGATGGGATCGATATTCGGCTTGAGATCGGGGTGCCCAATGTCGAAGATTTCCGGCGCCACGGCAAACTCCCCTCCCGGCTCGAACTGCCCTATCCTCACCCCCTTGCCCGTGAAGTCATGCCATACCGGCGCGACATTGATCTCGTCCAGATACCACTGCCGCCCCGAC > NC_002947/5600502‑5600773 |
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: ‑ |