Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,815,587 | C→T | 12.8% | intergenic (+15/+76) | PP_5090 → / ← PP_5091 | cell division protein/membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,815,587 | 0 | C | T | 12.8% | 93.4 / 6.3 | 39 | intergenic (+15/+76) | PP_5090/PP_5091 | cell division protein/membrane protein |
Reads supporting (aligned to +/- strand): ref base C (19/15); new base T (2/3); total (21/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 |
AACCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCA > NC_002947/5815449‑5815733 | aaCCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGc < 1:413076/150‑1 (MQ=255) ccTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGa > 2:191429/1‑150 (MQ=255) ggTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTCACAGCGGCACACAAGCCGCTGGCTGGGGCCGGGCTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACaa < 1:478794/150‑1 (MQ=255) gcaacgcgTGCTGGTTGGGCCGTCCAGCAACCGCGAGCAGTTGACAGGGGAACAGAAGCAGGTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCAGCGCTGCGGCAAg < 2:222196/146‑1 (MQ=255) gCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCAAGGCcttcttc < 1:327858/150‑2 (MQ=39) tGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTGAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtat > 2:329840/1‑150 (MQ=255) gggCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCaa < 2:464246/133‑1 (MQ=255) gggCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCaa > 1:464246/1‑133 (MQ=255) aCCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGt < 1:290224/150‑1 (MQ=255) aCCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTAAACAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGt > 1:334954/1‑150 (MQ=255) ccGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTg > 2:230732/1‑150 (MQ=255) cgcgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg > 1:223024/1‑119 (MQ=255) cgcgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg < 2:223024/119‑1 (MQ=255) cgcgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGa > 1:472255/1‑150 (MQ=255) cAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGAAAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAag < 1:191429/150‑1 (MQ=255) aGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagc > 1:210256/1‑150 (MQ=255) gggCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagcag < 2:82556/150‑1 (MQ=255) tggctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAg < 2:142749/150‑1 (MQ=255) ggTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCa < 2:472255/150‑1 (MQ=255) gTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCAt > 2:347658/1‑150 (MQ=255) aGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTatat < 1:261174/84‑1 (MQ=255) aGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTatat > 2:261174/1‑84 (MQ=255) gCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGccc > 2:171850/1‑97 (MQ=255) gCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGccc < 1:171850/97‑1 (MQ=255) ccTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCGCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCgg > 2:223345/1‑150 (MQ=255) ccTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCgg > 1:361805/1‑150 (MQ=255) tgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAgg > 1:399986/1‑92 (MQ=255) tgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAgg < 2:399986/92‑1 (MQ=255) aGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCCACCCCCACGAAGCGAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTc > 1:233427/1‑150 (MQ=255) ggCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGACGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGtt < 1:438215/146‑1 (MQ=255) ggCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGACGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGtt > 2:438215/1‑146 (MQ=255) tCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagca > 2:312367/1‑85 (MQ=255) tCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAagca < 1:312367/85‑1 (MQ=255) tAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCa > 2:348973/1‑141 (MQ=255) tAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCa < 1:348973/141‑1 (MQ=255) ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGTATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccac > 1:303658/1‑150 (MQ=255) cTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt > 1:216375/1‑103 (MQ=255) cTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt < 2:216375/103‑1 (MQ=255) ttGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCccacca > 2:272616/1‑150 (MQ=255) | AACCTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCA > NC_002947/5815449‑5815733 |
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: ‑ |