Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 813,203 | T→G | 7.0% | intergenic (+107/+24) | PP_0699 → / ← PP_t70 | LysE family transporter/tRNA‑Met |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 813,203 | 0 | T | G | 7.0% | 161.7 / 5.1 | 57 | intergenic (+107/+24) | PP_0699/PP_t70 | LysE family transporter/tRNA‑Met |
Reads supporting (aligned to +/- strand): ref base T (24/29); new base G (2/2); total (26/31) | |||||||||||
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.51e-01 |
CGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTT > NC_002947/813057‑813351 | ggggCTTGGCCTGGGTGTGGCGGTGACCGGAATACGTGGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc < 2:714/149‑1 (MQ=255) cTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCt < 2:282800/150‑1 (MQ=255) ggCTGTGGCGGTGACCGGACTACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGCCCttt < 2:1079679/150‑1 (MQ=255) ggCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCttt < 2:694812/150‑1 (MQ=255) cGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGCCCttt > 2:1303026/1‑142 (MQ=255) cGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGCCCttt < 1:1303026/142‑1 (MQ=255) tGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtgg < 1:354301/150‑1 (MQ=255) gACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt < 1:145273/150‑1 (MQ=255) cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGc > 2:615276/1‑150 (MQ=255) cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTACGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc < 1:163684/123‑1 (MQ=255) cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTACGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGc > 2:163684/1‑123 (MQ=255) gAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTa < 2:1483263/150‑1 (MQ=255) gTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTc < 2:242356/130‑1 (MQ=255) gTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTc > 1:242356/1‑130 (MQ=255) aTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAcccc > 2:641303/1‑150 (MQ=255) gcgcCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGAccc < 2:604239/150‑1 (MQ=255) cGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa < 1:871609/148‑1 (MQ=255) cGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa > 2:871609/1‑148 (MQ=255) gAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCa > 2:1350135/1‑150 (MQ=255) aaCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCAtt < 1:615276/150‑1 (MQ=255) ccTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt < 2:1503067/150‑1 (MQ=255) tgtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg > 1:707052/1‑121 (MQ=255) tgtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg < 2:707052/121‑1 (MQ=255) tgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAat < 2:660687/150‑1 (MQ=255) tgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAat < 2:1378543/150‑1 (MQ=255) ttCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgct > 1:1094134/1‑150 (MQ=255) cccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCAATATGAATGCTATGctct < 1:641303/150‑1 (MQ=255) tgctgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATTCTCTAACCaa > 2:419518/1‑150 (MQ=255) aTTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt > 1:1562559/1‑121 (MQ=255) aTTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt < 2:1562559/121‑1 (MQ=255) gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACCTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAg > 2:415236/1‑150 (MQ=255) cgcgGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGAt < 2:1454299/150‑1 (MQ=255) gCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatta < 2:1485941/150‑1 (MQ=255) gCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCatta > 2:1369897/1‑150 (MQ=255) cTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtg > 2:1389329/1‑60 (MQ=255) cTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtg < 1:1389329/60‑1 (MQ=255) tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCa > 2:695123/1‑92 (MQ=255) tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCa < 1:695123/92‑1 (MQ=255) tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtctc < 1:95196/150‑1 (MQ=255) cGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACCTGAACCCCGGACCCCAGCATTACGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGaa < 1:415236/150‑1 (MQ=255) gCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAGGCCGCGCCGCGGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAAc > 2:55775/1‑150 (MQ=255) gCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAAc > 1:546939/1‑103 (MQ=255) gCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAAc < 2:546939/103‑1 (MQ=255) cAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg > 2:20010/1‑58 (MQ=255) cAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg < 1:20010/58‑1 (MQ=255) cAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAGCCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGcccc > 1:122947/1‑150 (MQ=255) caaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGccccc > 1:827864/2‑150 (MQ=255) aaaGGAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGccccc > 2:1113549/1‑150 (MQ=255) aaaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGCCTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGccccc > 1:81044/1‑150 (MQ=255) aaaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGccccc > 2:324209/1‑150 (MQ=255) gCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCaa < 1:356077/93‑1 (MQ=255) gCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCaa > 2:356077/1‑93 (MQ=255) cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt < 2:1624016/72‑1 (MQ=255) cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt > 1:1624016/1‑72 (MQ=255) cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTa > 2:1128697/1‑104 (MQ=255) cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTa < 1:1128697/104‑1 (MQ=255) cTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCAttt < 1:1113549/150‑1 (MQ=255) ttGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAg < 1:1432269/91‑1 (MQ=255) ttGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAg > 2:1432269/1‑91 (MQ=255) | CGGGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTT > NC_002947/813057‑813351 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |