Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,380,752 | C→A | 13.7% | intergenic (‑134/+1) | PP_t67 ← / ← PP_t66 | tRNA‑Met/tRNA‑Met |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,380,752 | 0 | C | A | 13.7% | 103.7 / 4.7 | 42 | intergenic (‑134/+1) | PP_t67/PP_t66 | tRNA‑Met/tRNA‑Met |
Reads supporting (aligned to +/- strand): ref base C (17/19); new base A (4/2); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTTTGATTTGTCGATGTTTTTTC > NC_002947/1380606‑1380900 | gctCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATACTTATCAGCAGGCAAAGAAAAAGGCCACTCTGTTA‑ACAGTGGCCTTTTCTTGAATCTgg < 1:933136/150‑1 (MQ=255) tCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGtt > 1:134069/1‑150 (MQ=255) ccGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGa > 1:787379/1‑150 (MQ=255) aGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGAt < 2:46054/150‑1 (MQ=255) ggCGGTTTTTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGaa < 2:1071201/150‑1 (MQ=255) cGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGTAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAAcc < 1:421724/150‑1 (MQ=255) cGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGACTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAAcc < 1:1280157/150‑1 (MQ=255) tgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGa > 1:1083111/1‑120 (MQ=255) tgGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGa < 2:1083111/120‑1 (MQ=255) ggTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGACCTGAGGACCTTCGGGTTAt < 1:1118783/150‑1 (MQ=255) gTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg > 2:154405/1‑144 (MQ=255) gTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCgg < 1:154405/144‑1 (MQ=255) ttGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCCCCGGGTTA‑ACAGTGGCCTTTTCTTGAATATGGTTGCGGGTGCCGGATTTGAACCGACGACCTTCGGGGTATAAg > 1:318747/1‑150 (MQ=18) ttGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAg > 1:720168/1‑150 (MQ=255) cAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCcgacga < 2:1027766/150‑1 (MQ=255) gcTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTAc < 2:787379/150‑1 (MQ=255) tAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCa < 1:829326/150‑1 (MQ=255) gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAgg > 2:1292027/1‑136 (MQ=255) gACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAgg < 1:1292027/136‑1 (MQ=255) aaTGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGa > 2:634231/1‑150 (MQ=255) ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGccc > 1:522688/1‑103 (MQ=255) ttATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGccc < 2:522688/103‑1 (MQ=255) ggCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTa < 1:634231/150‑1 (MQ=255) gCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAg < 2:318747/150‑1 (MQ=255) accccgcctgtgtaaaagcagcgcggggtttgttctgacGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTg < 1:1418501‑M2/88‑1 (MQ=255) accccgcctgtgtaaaagcagcgcggggtttgttctgacGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTg > 2:1418501‑M2/40‑127 (MQ=255) cccgcctgtgtaaaagcagcgcggggtttgttctgacGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt < 1:300921‑M2/47‑1 (MQ=255) cccgcctgtgtaaaagcagcgcggggtttgttctgacGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt > 2:300921‑M2/38‑84 (MQ=255) aagaaaaGGCCACTGT‑TAACACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTta > 2:1208717/4‑150 (MQ=17) ccgcctgtgtaaaagcagcgcggggtttgttctgacGAATATGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGgtttg > 1:918064‑M2/37‑149 (MQ=255) cACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAgc > 2:824183/1‑150 (MQ=255) aCTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGa > 1:75959/1‑74 (MQ=40) cTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTa > 2:1257876/1‑124 (MQ=255) cTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTa < 1:1257876/124‑1 (MQ=255) ttA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTggg < 2:1104369/109‑1 (MQ=255) ttA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTggg > 1:1104369/1‑109 (MQ=255) a‑aCAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCgg > 1:73309/1‑38 (MQ=255) tGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTTTGAtt > 2:76117/1‑150 (MQ=255) gCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTAt > 2:1290945/1‑57 (MQ=38) cTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGGGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGa > 2:1135234/1‑98 (MQ=255) cTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGGGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGa < 1:1135234/98‑1 (MQ=255) tCTAGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTTTGATTTGTCGATGTTTTTTc < 2:720168/150‑1 (MQ=255) | GCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGTGTGTGGTTTGTATGTCAAGCGCTAGGCTTTGATTTTTCGACGTTTTTTCGCGAATGCTTATCAGCAGGCAAAGAAAAAGGCCACTGTGTTA‑ACAGTGGCCTTTTCTTGAATCTGGTTGCGGGAGCCGGATTTGAACCGACGACCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCATCCCGCGTCAGTGGGGCGGATTCTACAGCGTGTAGTTTATATGTCAAGCGCTAGGCTTTGATTTGTCGATGTTTTTTC > NC_002947/1380606‑1380900 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |