Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,256,277 | Δ1 bp | 12.9% | intergenic (+294/‑410) | asd → / → PP_1990 | aspartate‑semialdehyde dehydrogenase/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,256,277 | 0 | T | . | 12.9% | 81.6 / 7.9 | 31 | intergenic (+294/‑410) | asd/PP_1990 | aspartate‑semialdehyde dehydrogenase/transposase |
Reads supporting (aligned to +/- strand): ref base T (13/14); new base . (2/2); total (15/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01 |
CGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAG > NC_002947/2256139‑2256412 | cgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGtc < 1:431430/150‑1 (MQ=38) gttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGcagc > 1:216049‑M2/1‑139 (MQ=255) gttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGcagc < 2:216049‑M2/143‑5 (MQ=255) cttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCGCACCGGGAGc > 1:334400/1‑150 (MQ=12) cttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGccc < 2:273593/150‑1 (MQ=32) ttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCg < 2:209741/150‑1 (MQ=32) cTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa < 2:544745/150‑1 (MQ=32) cAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAg > 2:134269/1‑150 (MQ=32) cGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCTGGGGCC‑agcgcttgac > 2:220235‑M2/1‑97 (MQ=255) cGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCTGGGGCC‑agcgcttgac < 1:220235‑M2/107‑11 (MQ=255) aTTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAggg < 1:134269/150‑1 (MQ=35) aTTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAggg < 1:254218/150‑1 (MQ=35) ttCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatg > 1:251147/1‑150 (MQ=35) aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTa < 2:90242/131‑1 (MQ=33) aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTa > 1:90242/1‑131 (MQ=33) aggaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGatgatg > 1:174510/1‑150 (MQ=35) gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGatgatg > 1:13827/1‑148 (MQ=35) gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGatgatg < 2:13827/148‑1 (MQ=35) gggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAggg > 2:326365/1‑124 (MQ=35) gggCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAggg < 1:326365/124‑1 (MQ=35) gagggaggTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCtt > 2:520157/1‑150 (MQ=255) ccGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGcagcgcttgacgggcatccatgcccgtcacctccct > 1:103139‑M2/1‑32 (MQ=255) ccGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGcagcgcttgacgggcatccatgcccgtcacctccct < 2:103139‑M2/68‑37 (MQ=255) aGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTaa < 2:334400/150‑1 (MQ=255) cgcTG‑‑GGCCCCAGGATGGGGCC‑agcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgcgcgcttagaagcaaaagccaaagcaagagc > 1:307837‑M2/1‑23 (MQ=255) cgcTG‑‑GGCCCCAGGATGGGGCC‑agcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcgtcgcctgaactaccgcgcgcttagaagcaaaagccaaagcaagagc < 2:307837‑M2/145‑123 (MQ=255) tG‑‑GGCCCCAGGATGGGGCCTgcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcatcgcctgaactaccgcgcgcttagaa > 1:75144‑M2/1‑21 (MQ=255) tG‑‑GGCCCCAGGATGGGGCCTgcgcttgacgggcatccatgcccgtcacctccctccgcaagccctgcgctcggcctcctgaagtcgcaatctgcatcgcctgaactaccgcgcgcttagaa < 2:75144‑M2/122‑102 (MQ=255) gCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGCGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAg > 2:170795/1‑143 (MQ=255) gCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAg < 1:170795/143‑1 (MQ=255) cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAg > 2:176336/1‑150 (MQ=255) | CGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑‑GGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAG > NC_002947/2256139‑2256412 |
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: ‑ |