Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,039,700 | G→T | 7.4% | I284I (ATC→ATA) ‡ | PP_4446 ← | group II intron‑encoding maturase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,039,700 | 0 | G | T | 7.4% | 161.4 / 2.7 | 54 | I284I (ATC→ATA) ‡ | PP_4446 | group II intron‑encoding maturase |
Reads supporting (aligned to +/- strand): ref base G (28/22); new base T (2/2); total (30/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01 |
GATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAAC > NC_002947/5039552‑5039844 | gATTACGCTGGCTACGACGTTTGAGGGCTTGCCACCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGa > 1:770198/1‑150 (MQ=255) tACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATgg < 1:714247/150‑1 (MQ=255) gCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGtt > 2:762558/1‑150 (MQ=255) tGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCa < 1:908714/150‑1 (MQ=255) tGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCa < 1:846076/150‑1 (MQ=255) ggCTACGACGTTTGAGGGCTTGCCACCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCAt < 2:770198/150‑1 (MQ=255) gCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaag < 1:762558/150‑3 (MQ=255) gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTc > 1:681196/1‑150 (MQ=255) acgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCg > 2:724381/1‑150 (MQ=255) ggCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAg > 1:117385/1‑150 (MQ=255) gccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATccc < 1:518294/150‑1 (MQ=255) ccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTTAACCCATCTGGGTTCaagg < 1:839343/127‑4 (MQ=255) ccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTTAACCCATCTGGGTTCaagg > 2:839343/1‑124 (MQ=255) cATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATg < 1:1144810/125‑1 (MQ=255) cATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATg > 2:1144810/1‑125 (MQ=255) aTAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCCGTGTTTCAGCTCATCCCGCa > 2:184526/1‑150 (MQ=255) gCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCcag < 1:724381/150‑1 (MQ=255) cAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCTATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCcagc < 2:376010/150‑1 (MQ=255) aCCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCagcagc > 2:509991/1‑150 (MQ=255) gccgAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATc > 2:723716/1‑150 (MQ=255) cccAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGcttc > 2:612854/1‑150 (MQ=255) cGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGAcc < 1:1163977/150‑1 (MQ=255) gTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCg > 2:525820/1‑150 (MQ=255) tttCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGc > 1:1013861/1‑150 (MQ=255) ccGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCa < 1:602987/141‑1 (MQ=255) ccGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCa > 2:602987/1‑141 (MQ=255) cGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTAt > 2:881219/1‑150 (MQ=255) aCCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGt < 1:237350/134‑1 (MQ=255) aCCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGt > 2:237350/1‑134 (MQ=255) aCCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGa < 1:612854/150‑1 (MQ=255) tGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCa < 1:881219/150‑1 (MQ=255) tAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCg > 1:455385/1‑150 (MQ=255) aTTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCt > 2:1014922/1‑150 (MQ=255) aaGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTg < 2:1013861/150‑1 (MQ=255) aaCCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg < 2:681196/150‑1 (MQ=255) cGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATTACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGtac > 1:579160/1‑150 (MQ=255) tCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGtactac < 2:579160/150‑1 (MQ=255) accaccCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGaa > 2:765991/1‑124 (MQ=255) accaccCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGaa < 1:765991/124‑1 (MQ=255) accaccCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTacaa > 1:1157499/1‑150 (MQ=255) ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGcttctt > 2:616406/1‑96 (MQ=255) ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGcttctt < 1:616406/96‑1 (MQ=255) ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTg > 2:271784/1‑150 (MQ=255) ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTg > 1:529033/1‑150 (MQ=255) cGTGTGAGCCATTGCCCTTGACCCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGc > 1:37158/1‑150 (MQ=255) tgAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGa > 1:900807/1‑150 (MQ=255) ccATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGt > 1:1039786/1‑125 (MQ=255) ccATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGt < 2:1039786/125‑1 (MQ=255) ccATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGcc > 2:681992/1‑150 (MQ=255) cATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGccc < 1:525820/150‑1 (MQ=255) ttGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAg < 2:529033/150‑1 (MQ=255) cccTTGAACCCAGATTGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGaa > 2:837143/1‑150 (MQ=255) cccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCCGTAATTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGCTGCCAAggca > 2:216630/1‑148 (MQ=255) aCCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAAc < 1:1014922/150‑1 (MQ=255) | GATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAAC > NC_002947/5039552‑5039844 |
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: ‑ |