Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,879,891 | C→A | 11.9% | intergenic (+75/+61) | PP_1686 → / ← PP_1687 | glutathione peroxidase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,879,891 | 0 | C | A | 11.9% | 114.4 / 8.0 | 42 | intergenic (+75/+61) | PP_1686/PP_1687 | glutathione peroxidase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (18/19); new base A (2/3); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01 |
TCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCA > NC_002947/1879745‑1880037 | tcCAGGATGACCAACCGGGTTGATCCCGGGTTCGGCGCTGCTATCTAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGGGGGGTTGCTTTAAAGCGCGAGGCAGTTTAGGCGACGCAGATTGCGACTTCAGGAGGCCg < 1:189358/149‑1 (MQ=37) cTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTcc < 6:158876/149‑1 (MQ=255) aCCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTg > 8:108372/1‑148 (MQ=255) ccAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTg > 7:6346/1‑147 (MQ=255) aaaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCAACTTCAGGAGAACGGCCTCCTGAAGt < 5:138382/149‑1 (MQ=39) aaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCAACTTCAGGAGAACGGCCTCCTGAAGtc < 3:113567/149‑2 (MQ=39) ccGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGc < 7:108372/149‑1 (MQ=255) ggATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCAGTTCTCCTGAAGTTGCgg < 5:325530/149‑1 (MQ=255) ccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc > 1:313061/1‑145 (MQ=255) cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagat < 2:249840/149‑1 (MQ=255) cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc > 2:88902/1‑144 (MQ=255) gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga < 5:148125/147‑1 (MQ=255) gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga > 6:148125/1‑147 (MQ=255) gTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatcaa < 4:150037/149‑1 (MQ=255) tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtg > 2:11481/1‑149 (MQ=255) aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg > 2:59120/1‑149 (MQ=255) aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg > 7:260323/1‑149 (MQ=255) aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg < 1:59120/149‑1 (MQ=255) cTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACtt > 4:45219/1‑149 (MQ=255) cTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACt > 4:103057/1‑148 (MQ=255) acaGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCa < 1:284091/148‑1 (MQ=255) aCCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtgcatg < 8:260323/148‑1 (MQ=255) ggCTGCGCTTGCTGAGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACgg < 3:103057/149‑1 (MQ=255) cgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGtgaa < 2:313061/148‑4 (MQ=255) cgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaa < 3:164430/147‑1 (MQ=255) tgcGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGccacc > 7:5865/1‑148 (MQ=255) cggcgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAAc < 1:11481/149‑1 (MQ=255) gTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGca < 3:45219/149‑1 (MQ=255) ttGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcac < 8:6346/149‑1 (MQ=255) cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac > 7:38888/1‑149 (MQ=255) cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac > 7:84965/1‑149 (MQ=255) cTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccac > 8:183300/1‑149 (MQ=255) aaaGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccaccac < 1:88902/148‑1 (MQ=255) gcgAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggccg < 7:183300/149‑1 (MQ=255) cGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacaca > 8:218703/1‑149 (MQ=255) cccgcaaCTTCAAGAGAACGGCCTCCTGAAGTCGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATg > 3:81085/7‑149 (MQ=255) gcaaCTTCAGGAGAACGGCCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGac < 4:148302/123‑1 (MQ=255) gcaaCTTCAGGAGAACGGCCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGac > 3:148302/4‑126 (MQ=255) ttCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCgg > 4:182760/1‑148 (MQ=255) gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca < 4:269437/118‑1 (MQ=255) gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca > 3:269437/1‑118 (MQ=255) gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCa > 1:109766/1‑148 (MQ=255) | TCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCA > NC_002947/1879745‑1880037 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |