Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 4,586,056:1 | +C | 100% | intergenic (+139/+76) | PP_4061 → / ← PP_4063 | conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 4,586,056 | 1 | . | C | 91.3% | 148.5 / 8.0 | 46 | intergenic (+139/+76) | PP_4061/PP_4063 | conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (4/0); new base C (19/23); total (23/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAT > AE015451/4585824‑4586263 | aCCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc < 1:1109245/236‑2 (MQ=255) ccGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc < 1:59786/236‑1 (MQ=255) ccgccgACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg < 2:897433/211‑1 (MQ=255) aCAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt < 2:818289/216‑1 (MQ=255) cgccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc < 2:1050340/233‑1 (MQ=255) gccgGGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 2:1124787/213‑1 (MQ=255) ccAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggccgg < 1:873815/228‑1 (MQ=255) aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGACGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAgg < 1:946782/179‑1 (MQ=255) gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt < 2:988747/186‑1 (MQ=255) gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTACACCCGCTCCTACCAGGGACC‑GTGCCGGCCtt > 2:425035/1‑151 (MQ=255) aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATGTTGCACTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgc > 2:1070393/1‑150 (MQ=38) gCACCGGCGGGAACAGCTGAGCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc < 2:269196/130‑2 (MQ=38) aCCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑TTGCCGGCCTTGCATgcggcg > 2:306325/1‑149 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGc < 2:948472/201‑1 (MQ=255) tgctCCATTGAAAGGTTGAGGATGACCGGGGGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGac > 2:322638/1‑155 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAG‑GGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTTGGAGCGGGTTTAcc > 1:172976/1‑145 (MQ=255) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACGCGCTCCTACCAGCGA‑CTGTGCCGGcc > 2:899678/1‑107 (MQ=37) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCaa < 1:528895/194‑1 (MQ=255) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACGCTCTCGCAGCCGCGA‑CTGATCCGGCCTTGAATGcgccgcaccc > 1:830496/1‑115 (MQ=255) cgggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGTTCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACAt < 1:681833/165‑1 (MQ=38) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACTCTCTCCTACCAGGGA‑CTGTACCTGCCTtgcgtgcgc > 1:753087/1‑106 (MQ=255) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa > 1:328342/1‑216 (MQ=255) ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGCAGCGGGTTTCCCCGCGCAGAGGCCGGCCGCGCCac > 2:59240/1‑157 (MQ=38) aGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGAtc < 2:396618/199‑1 (MQ=255) aGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGAtc > 1:396619/1‑199 (MQ=39) gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcggtcacggt > 1:595662/1‑235 (MQ=255) gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTAACAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt > 1:65218/1‑109 (MQ=37) gggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTg > 2:175958/1‑88 (MQ=21) gggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGcccg > 2:373355/1‑144 (MQ=37) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCCGGGACC‑GTGCCGGCCTTGAATGCGCCGCTGCCCCTGGTAGGAGc > 1:47071/1‑113 (MQ=37) ggTGCTTGGCAAAGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa > 2:632470/1‑203 (MQ=39) gATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAATCCGCTCCTAGCAGGGA‑CTGTGCATGCCTTGGATgc > 1:634920/1‑74 (MQ=9) ttGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTAc < 1:960940/215‑1 (MQ=255) ggccggccTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcgg < 1:1070394/186‑1 (MQ=255) ccggccTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCtt < 2:595661/233‑1 (MQ=255) cttcgcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGCATATACcg > 1:297461/8‑217 (MQ=255) ttcgcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg < 1:897187/148‑1 (MQ=255) ttcgcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg > 2:897186/7‑154 (MQ=255) gcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATAcc < 1:851196/209‑1 (MQ=255) aaaCCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcggtc < 1:632471/169‑1 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGcc > 1:942302/1‑116 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGcc < 2:942301/116‑1 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACtcgtc > 2:448834/1‑178 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACtcgtc < 1:448835/178‑1 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGtatata < 1:479241/193‑1 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTTTGCGAt > 2:1079809/1‑217 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAt < 1:1079810/217‑1 (MQ=255) cAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAAc < 2:635040/183‑1 (MQ=255) | ACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAT > AE015451/4585824‑4586263 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |