Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 850,455 | C→G | 9.5% | intergenic (‑35/‑770) | pcs ← / → hemA | phosphatidylcholine synthase/glutamyl‑tRNA reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 850,455 | 0 | C | G | 9.5% | 175.2 / 8.6 | 61 | intergenic (‑35/‑770) | pcs/hemA | phosphatidylcholine synthase/glutamyl‑tRNA reductase |
Reads supporting (aligned to +/- strand): ref base C (31/24); new base G (2/4); total (33/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01 |
TAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTCTAA > NC_002947/850311‑850603 | tAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGc < 1:2050568/150‑1 (MQ=255) tataGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATAc < 1:691754/150‑1 (MQ=255) taaataTAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAAt > 2:1590011/1‑150 (MQ=255) ataGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTTAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGttt < 1:82243/150‑1 (MQ=255) cTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTaa > 1:1285516/1‑150 (MQ=255) cGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATaa < 2:2219371/150‑1 (MQ=255) gagGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCtt < 2:2324838/150‑1 (MQ=255) gCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGaa < 1:1586626/109‑1 (MQ=255) gCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGaa > 2:1586626/1‑109 (MQ=255) cGAAATAGCACACAAGTTAACGACCAAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCGTGCATTGGAAATCGCTTTTAAGATAAGCTTAACac > 2:862081/1‑149 (MQ=255) aTAGCACACAAGTTAACGCCCCAAATAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACt < 1:79470/150‑1 (MQ=255) gCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACACTACAATGGAAATCGTTTTTAAGATAAGCATAACAAAACTGaa > 1:1352537/1‑150 (MQ=255) gTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATaaaa > 2:2102860/1‑150 (MQ=255) cccAAAGAAGCACTGAGTGCCTACTAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttct > 2:1236396/1‑150 (MQ=255) cAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCAATG‑CACGAAGCATAACCTGAATTAGAAATCATTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttg > 1:1159027/1‑150 (MQ=255) cctatctcTGAGTGCCTACGGGCTACTGCAGCTTGCTCCATTGCTCTCTCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttgttg < 2:613080/143‑1 (MQ=255) agaagCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttgttg < 2:1285516/150‑1 (MQ=255) cACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATATCCTGCATTGGAAATCATTTTTAAGATAAGATTAACAATTCTGATATATAAAAATGCTTCTTGTTGgtgtg > 1:1905550/1‑150 (MQ=255) tGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTCATAAGACTGAAATATAAAAATGCTTCTTGTTggagatcgg > 1:613080/1‑143 (MQ=255) tGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATCTGGAAATGCTTCTTGTTGGTGTGTCACCaa > 1:1756721/1‑150 (MQ=255) tGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaa < 1:1164160/150‑1 (MQ=255) gCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaag < 1:2102860/150‑1 (MQ=255) cTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAAAGCTTCTTGTTGGTGTGTCACCaagta > 1:1522999/1‑150 (MQ=255) cGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCt < 2:916236/98‑1 (MQ=255) cGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCt > 1:916236/1‑98 (MQ=255) agaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAgtgt < 2:855492/150‑1 (MQ=255) gaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc < 1:607148/123‑1 (MQ=255) gaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc > 2:607148/1‑123 (MQ=255) aTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGc < 2:29371/150‑1 (MQ=255) aCTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcat < 1:704871/150‑1 (MQ=255) tGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcatca < 1:1236396/150‑1 (MQ=255) gCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAACGTGGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcatcat > 1:21869/1‑150 (MQ=255) aaaTTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc > 1:1027079/1‑150 (MQ=255) aaaTTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAAAGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc < 2:1522999/150‑1 (MQ=255) aaacgatttccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTtcat < 1:891846‑M1/127‑1 (MQ=255) cATTGCTCGCGCCTATTATATAGATTAATGGCATG‑CAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTtcatc < 2:21869/150‑1 (MQ=255) tGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTaa > 2:1581685/1‑78 (MQ=255) tGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTaa < 1:1581685/78‑1 (MQ=255) cTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTAc > 1:1418175/1‑150 (MQ=255) tCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc > 2:1163682/1‑100 (MQ=255) tCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc < 1:1163682/100‑1 (MQ=255) ccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc > 2:754693‑M1/15‑133 (MQ=255) ccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCGTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc < 1:754693‑M1/119‑1 (MQ=255) gcCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctg > 2:1363430/1‑150 (MQ=255) atgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc > 2:891846‑M1/12‑150 (MQ=255) ttatATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTaca > 2:762315/1‑148 (MQ=255) atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc < 1:1158885/144‑1 (MQ=255) atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc > 2:1158885/1‑144 (MQ=255) atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAAcc > 2:1555070/1‑150 (MQ=255) gtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCt < 2:663758‑M1/146‑1 (MQ=255) aGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCtcat > 2:38586/1‑150 (MQ=255) GCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATGAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCtcatc > 1:152096‑M1/1‑150 (MQ=255) TTCCTGCCATG‑AACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGAAAAGAATAACAGTACGGAGGTAGAGAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCaa > 1:862081‑M1/1‑150 (MQ=255) gCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATaaaa < 2:1905550/150‑1 (MQ=255) aTG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACtt > 1:705056/1‑150 (MQ=255) tG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAAAAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTc > 1:1698437/1‑150 (MQ=255) aCGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaagtaagt < 2:1899131/94‑1 (MQ=255) aCGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaagtaagt > 1:1899131/1‑94 (MQ=255) aCGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTCTaa > 1:848088/1‑150 (MQ=255) aCGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTCTa < 2:2159780/149‑1 (MQ=255) aCGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTCTa > 1:2159780/1‑149 (MQ=255) | TAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAATAAAACTTCTAA > NC_002947/850311‑850603 |
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: ‑ |