Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,287,319 | T→C | 16.7% | intergenic (+17/+24) | PP_2889 → / ← xthA | transmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,287,319 | 0 | T | C | 16.7% | 64.6 / 6.5 | 30 | intergenic (+17/+24) | PP_2889/xthA | transmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI |
Reads supporting (aligned to +/- strand): ref base T (13/12); new base C (3/2); total (16/14) | |||||||||||
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.17e-01 |
GGACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGG > NC_002947/3287170‑3287460 | ggAGGGGCGGCGGGGGGGTGGCGGGCGGGGGGGGGGGTACTGGTGGTATGGGGTTTATATTTTGGCGGGGGTGAGGGCGGCGGGGGGGGGGCGTGGGGGGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGgtgt < 1:1501639/150‑1 (MQ=255) ggACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGgtgt > 1:1847263/1‑150 (MQ=255) gCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGTGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCt > 1:244777/1‑150 (MQ=255) gggTCGGTGCCCAGCGGGGGGCGCGGGACTGGTCGCGTGGGTCTTATCACTTGGCCTTGGTCAGCGCGGGGGACGGGGTGCGGGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGc < 1:1704548/150‑1 (MQ=255) ggggggaCTGGTGGCATGGGGCTTATCGCTTTGCCTTGGTCAGCGGGGCGGACGATGGGGGTGGGGGGGGGGTGGGGATGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCaa < 1:1379769/144‑1 (MQ=255) cTGGTGGCATGGGCTTTTTATGTTTTCCTTGGTAGGCGGGGGGGACGGGGGGGGGGGGGCGGGGGTGGAGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTccagc < 2:773995/150‑1 (MQ=255) tGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCTCGGCGGACGATGTGCGAGGGGCGGGGAAGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGCGCCC‑GCCACCCTAATCAGCTCAACTccagcc > 2:628093/1‑150 (MQ=255) gTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTccagccag > 2:1650397/1‑150 (MQ=255) gTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGAAGGCGTTGAGGTTGGGGTTGTACGTGTGGGCCCAAGGTGTGCACCTTGGTCCC‑GTACCCCTAATCAGCTCAATTccagccag > 2:7721/1‑150 (MQ=255) aTGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCgg > 1:1749434/1‑150 (MQ=255) tGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGTGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGgc < 2:244777/150‑1 (MQ=255) ataACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTc > 1:1582351/1‑150 (MQ=255) gggCCTTGGGCAGCGCGGCGGACGTTTGTTGGGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACgg < 2:1365225/148‑1 (MQ=255) gcggcgGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCcgc > 1:166143/1‑150 (MQ=255) cggGCGAGGTGCGGGGGGCGGGGGTGGGGGAGGTGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAg < 1:1251020/150‑1 (MQ=255) gACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGACTCCAGTCACTAT‑TTCAGCTCAACTCTAGCCAGATCGCCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCt > 1:1006054/1‑150 (MQ=255) gcggtggtgGGGGGGTGGTGGAGGCGTTGAGGGGGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCATCTCCAGGTAGAGCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTc < 2:758300/147‑1 (MQ=255) gggCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTACAAATTGATGACCCAAGGTGTGCACCTTGGGCCC‑TGCTCCCTAATACGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCACCTCATATACCACGcc > 1:1107142/1‑150 (MQ=255) cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGc > 2:1841984/1‑150 (MQ=255) cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGccca < 2:289699/66‑2 (MQ=25) cGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGCGGGCCCAAGGTGCACACCTTGGGccca > 1:289699/1‑65 (MQ=25) ggggggcgAGGGGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCCTGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGcc < 2:166143/147‑1 (MQ=255) ttGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCa > 2:1816284/1‑150 (MQ=255) gTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCtgt < 2:1582351/150‑1 (MQ=255) aGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGcc < 2:1107142/150‑1 (MQ=255) gtgGTCCCAAGGTGTGCACCTTGGGCCC‑GTCATCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATg < 2:144381/150‑1 (MQ=255) gggCCCAAGGTGCACACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAt > 2:1157158/1‑150 (MQ=255) gggCCCAAGGTGCACACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAt > 2:1010698/1‑150 (MQ=255) cccAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa < 2:156420/138‑1 (MQ=255) cccAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa > 1:156420/1‑138 (MQ=255) cccAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAg < 2:1444277/150‑1 (MQ=255) ccAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAgg > 1:954906/1‑150 (MQ=255) | GGACGGGCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCC‑GTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGG > NC_002947/3287170‑3287460 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |