Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 152,157 | A→T | 13.8% | intergenic (‑51/+43) | PP_0143 ← / ← PP_0144 | conserved membrane protein of unknown function/Metalloprotease, insulinase family |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 152,157 | 0 | A | T | 13.8% | 81.5 / 6.6 | 29 | intergenic (‑51/+43) | PP_0143/PP_0144 | conserved membrane protein of unknown function/Metalloprotease, insulinase family |
Reads supporting (aligned to +/- strand): ref base A (4/21); new base T (2/2); total (6/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.80e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01 |
CAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAGCTCGTCATAGCCTAGCAACGGCTTTTCGATGCGCACATAGCCATCTTCCTTGAGCAGATCCTTCAATGCC > AE015451/151943‑152377 | cAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCgg > 2:864931/1‑233 (MQ=255) ggACAAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAGGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTc < 2:717192/226‑1 (MQ=255) cAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa < 1:713334/236‑1 (MQ=255) aGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAg < 2:521811/212‑1 (MQ=255) ggACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGAtt < 1:718327/235‑1 (MQ=255) gccaggccagCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgc < 2:813850/189‑1 (MQ=255) ttGCCGATCCCGGAGAGGAAACCGCCGAGCGCAAGGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCg < 2:313478/165‑1 (MQ=255) cGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCg < 1:296490/193‑1 (MQ=255) aTACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTGGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACCAGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTa > 2:515463/1‑201 (MQ=255) cGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCg < 1:209076/187‑1 (MQ=255) caCCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCg < 2:34172/177‑1 (MQ=255) aaCCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGaacg < 2:581645/142‑3 (MQ=255) ccgccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg < 1:187944/184‑1 (MQ=255) aaaggatgCCCGGGTGGCGCGTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCgtg < 2:94644/170‑1 (MQ=255) cgcTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAgg < 1:450294/235‑1 (MQ=255) gcTTGAGTAGTCCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCcgc < 1:635818/107‑1 (MQ=255) ttGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCcgc > 2:806818/1‑105 (MQ=255) gcTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAgg < 1:468029/194‑1 (MQ=255) gcGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa < 1:858671/220‑1 (MQ=255) gggCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGc > 1:818870/1‑109 (MQ=255) tGGGGGCTACGACAGAAAAGTGCGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCTCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTa < 2:532678/198‑1 (MQ=255) ggggCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGc < 2:832741/165‑1 (MQ=255) cAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGAAAg > 1:523833/1‑51 (MQ=21) cAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCTTGCCCGCGAAAg < 2:523830/51‑1 (MQ=21) ggTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTaga < 1:34544/177‑1 (MQ=255) ctctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAGCTCGTCATAGCCTAGCAACGGCTTTTCGAt < 1:41439/193‑1 (MQ=255) ttCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAGCTCGTCATAGCCTAGCAACGGCTTTTCGATGCGCACATAGCCATCTTCCTTGAGCAGATCCTTCAATGcc < 2:572096/230‑1 (MQ=255) gcgGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCaa > 2:269668/1‑148 (MQ=255) ggCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGTTACCGCCGATTTACCGCTGCCGTGTAGCACTACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGATCAAGGCCATCAACCCGTAGAGCTCGTCATAGCCTAGCAACGGCTTTTCGATGCGCACATAGCCATCTTCCTTGAGCATATCCTTCa < 1:515466/219‑1 (MQ=255) | CAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCCATCAACCCGTAGAGCTCGTCATAGCCTAGCAACGGCTTTTCGATGCGCACATAGCCATCTTCCTTGAGCAGATCCTTCAATGCC > AE015451/151943‑152377 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |