Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 660,941 | A→T | 7.8% | intergenic (+73/‑94) | speA → / → PP_0568 | biosynthetic arginine decarboxylase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 660,941 | 0 | A | T | 7.8% | 143.5 / 3.3 | 51 | intergenic (+73/‑94) | speA/PP_0568 | biosynthetic arginine decarboxylase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (22/25); new base T (2/2); total (24/27) | |||||||||||
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.20e-01 |
GTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG > NC_002947/660801‑661087 | gTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCa > 7:224955/1‑149 (MQ=255) gAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc > 3:178987/1‑148 (MQ=255) gACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatc > 5:271227/1‑142 (MQ=255) gACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatc < 6:271227/146‑5 (MQ=255) gACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCTGCGt < 3:257799/149‑1 (MQ=255) gACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGt < 6:211/149‑1 (MQ=255) cccAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTca < 6:31535/149‑1 (MQ=255) tCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCg < 4:53168/148‑1 (MQ=255) gATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaa < 8:224955/148‑1 (MQ=255) cTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATcc < 1:60328/149‑1 (MQ=255) gcgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGc > 5:232171/1‑149 (MQ=255) ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAc > 5:277554/1‑149 (MQ=255) ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAc < 4:178987/149‑1 (MQ=255) gTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCaa > 8:43073/1‑149 (MQ=255) gTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGACCCTGGTCCATCGCACCAATGCCg < 7:43073/148‑1 (MQ=255) cgtcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatc > 6:305114/1‑89 (MQ=38) cgtcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatc < 5:305114/93‑5 (MQ=38) gtcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGa < 5:28867/149‑1 (MQ=255) tAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt > 4:97242/1‑143 (MQ=255) tAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt < 3:97242/143‑1 (MQ=255) tAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATg > 3:231442/1‑148 (MQ=255) ttGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt < 6:232171/149‑1 (MQ=255) gcaaggCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCCCTACCTTTGAATGCAATGta > 7:255618/5‑148 (MQ=25) gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCa > 4:1242/1‑148 (MQ=255) ggCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGACCTTGTCCATCGCACCAAAGCCGCTACCTATGAATGCAATCTTCa > 7:148599/1‑147 (MQ=25) gCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCTGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggag > 8:10115/1‑149 (MQ=255) tGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGa > 4:141931/1‑149 (MQ=255) cAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGata > 6:109727/1‑149 (MQ=255) cAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGa < 6:277554/147‑1 (MQ=255) aaaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCTGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGatat < 7:10115/149‑1 (MQ=255) aaaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGatat < 3:141931/149‑1 (MQ=255) gCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCaa > 7:159962/1‑149 (MQ=255) gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGat > 7:109180/1‑134 (MQ=255) gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGat < 8:109180/134‑1 (MQ=255) cTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTc < 8:159962/149‑1 (MQ=255) ttttGTTGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCa < 4:56020/112‑1 (MQ=255) ttttGTTGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCa > 3:56020/1‑112 (MQ=255) tttGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTcc > 5:35625/1‑148 (MQ=255) tgtgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCa < 3:1242/148‑1 (MQ=255) tgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATc < 2:159085/148‑1 (MQ=255) tgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATc > 1:159085/1‑148 (MQ=255) tgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATc < 4:231442/148‑1 (MQ=255) gggAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGt > 6:65024/1‑149 (MQ=255) ccgtGCGGCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTgg < 6:26724/146‑1 (MQ=255) tagCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCg > 1:240339/3‑149 (MQ=255) tCGCTATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCGAGGCTCCCATCGTTTACGTGGCGGTCGa < 6:4327/149‑1 (MQ=255) tCGCGATATGAGGTCGAAGCCCTCACCAGCTTCACTACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa < 8:255618/149‑1 (MQ=255) tatGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCt < 2:186664/149‑1 (MQ=255) tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa < 6:228381/141‑1 (MQ=255) tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa > 5:228381/1‑141 (MQ=255) cGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCgg > 1:33205/1‑149 (MQ=255) | GTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGG > NC_002947/660801‑661087 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |