| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,457,925 | C→A | 25.1% | H93N (CAC→AAC) | PP_1276 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,457,925 | 0 | C | A | 25.1% | 42.3 / 13.5 | 28 | H93N (CAC→AAC) | PP_1276 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (10/11); new base A (3/4); total (13/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02 | |||||||||||
TACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTC > NC_002947/1457781‑1458074 | tACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAg > 1:21021/1‑150 (MQ=255) gggAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgc < 1:305332/150‑1 (MQ=255) ggAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGcgcg < 2:386196/150‑1 (MQ=255) ggAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc < 1:215164/150‑2 (MQ=255) aGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGc < 1:640415/150‑1 (MQ=255) cGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc < 1:346960/143‑2 (MQ=255) cGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc > 2:346960/1‑142 (MQ=255) gcgTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATcc > 1:632051/1‑150 (MQ=255) tACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCggg > 2:543492/1‑150 (MQ=255) tACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCggg < 1:543492/150‑1 (MQ=255) gAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTc > 1:263728/1‑150 (MQ=255) tACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTACcaca > 2:97914/1‑150 (MQ=255) cATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAg < 1:274104/150‑1 (MQ=255) cATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAg < 1:97914/150‑1 (MQ=255) cATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAg > 2:274104/1‑150 (MQ=255) cATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAg > 2:222745/1‑150 (MQ=255) tGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAAt > 1:474079/1‑150 (MQ=255) aggGCGTGCGGGATGTCGTCATCGGGGCTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTg > 1:451831/1‑150 (MQ=255) gtcgtcATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGgatgat < 1:768882/150‑1 (MQ=255) tcgtcATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGAtt < 1:222745/150‑1 (MQ=255) tcATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc > 1:362312/1‑41 (MQ=25) tcATCGGGGGTGAGGCGAACCGGGCGCTTGTGAACAAgcgcc < 2:362312/42‑2 (MQ=25) tCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGc < 1:645666/150‑1 (MQ=255) gCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAgg > 1:810722/1‑150 (MQ=255) cGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATGCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAggg < 2:632051/150‑1 (MQ=255) ccGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGggtggt < 2:263728/150‑1 (MQ=255) ccgcgcTTGTGAACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGtggtgg > 1:321037/3‑150 (MQ=255) cgcTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTggtgggtg < 2:451831/150‑1 (MQ=255) cacaAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTc > 2:202584/1‑150 (MQ=255) | TACGGGAGCGCCGAAAAATCAGCGTGGTTGTACGCGTACAAGGGGCAGGCCGAAGATGCGGACGTCTACATTGACTTCGATCTGGAACTACAGGAGGGCGTGCGGGATGTCGTCATCGGGGGTGAGGCGAACCGGGCGCTTGTTCACAAGCGCGGCACAACCTATGGCGGCTATGCGAAATCCGGGCGTATCCGGAAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTC > NC_002947/1457781‑1458074 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |