Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,256,199 | G→C | 8.5% | intergenic (+259/‑445) | 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,199 | 0 | G | C | 8.5% | 115.6 / 2.5 | 45 | intergenic (+259/‑445) | asd/PP_1990 | aspartate‑semialdehyde dehydrogenase/transposase |
Reads supporting (aligned to +/- strand): ref base G (22/18); new base C (2/2); total (24/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA > NC_002947/2256057‑2256342 | attgCGAAATTTCTGCTCTGACAGAGCCGTTGCAGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGccc < 6:120123/146‑1 (MQ=255) gCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGt > 5:122817/1‑148 (MQ=255) aaaTGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCacc < 8:251576/149‑3 (MQ=255) tGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCAggc > 3:217761/1‑148 (MQ=255) ctctGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg < 6:233177/146‑1 (MQ=255) gagCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATg < 6:122817/149‑1 (MQ=255) ccgttgccgtcgccgttgCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATgggg < 2:215867/149‑1 (MQ=255) ccgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGC‑AGGCCCCATCCTGGGGcccagcg < 4:217761/149‑6 (MQ=14) cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt > 8:146654/1‑112 (MQ=255) cgttgccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt < 7:146654/112‑1 (MQ=255) ccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGTCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGcgcg < 2:900/148‑1 (MQ=37) ccgttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCTT‑‑‑Gctcccgg > 3:38991/1‑142 (MQ=255) gttgcTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGtc < 4:92443/149‑1 (MQ=35) cttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGtcctc < 4:198143/148‑1 (MQ=34) ttttgcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTcc‑‑c > 7:6765/1‑149 (MQ=34) gcttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑Cgg > 8:26171/1‑147 (MQ=32) cttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGa < 7:26171/148‑1 (MQ=32) ttttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAg > 5:305313/1‑148 (MQ=32) tttgcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGc > 5:314571/1‑148 (MQ=32) gcttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGc > 5:97079/1‑149 (MQ=32) ttCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGccc > 3:217945/1‑149 (MQ=32) tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCg > 5:138586/1‑147 (MQ=32) tAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGa < 6:305313/149‑1 (MQ=32) gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGTATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg < 4:145396/149‑1 (MQ=17) gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg < 4:217945/149‑1 (MQ=32) gcgGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAggg > 2:293671/1‑149 (MQ=32) gCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCg > 1:267569/1‑149 (MQ=32) gACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCg > 1:268407/1‑149 (MQ=34) aggcggCCGTGCGCTGGGTTGGGGGCGGTTGGTGAGGGGCAGGGATCCCCGTAAAGCGGTGGGCCCCTGTTTGGGGCCT‑‑‑GCTGCGGGGTCCTCC‑‑CGGTAGCAAGCCGGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGgatgat < 8:6765/149‑1 (MQ=2) gaggCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATgg > 3:160175/1‑149 (MQ=34) tggccttgcgcTGGGCTTGCGGTGGGAGGTGACGGGCAGGGTGGCCCGTCAAGCGCTGGGCCCCTGGTTGGGGCCG‑‑‑GCAGGGGGGTCCTCC‑‑CGGGAGAAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATggg < 6:97079/142‑1 (MQ=11) ccGTGCGCTTGGTTTGGGGTGGGTGGTTCCGGGCTTGGATGCCCGTCTAGCGTGGGCCCCCTGTAGGGGGCCG‑‑‑GCAGCGGGGTCCACC‑‑CGGCAGCAAGCCGGGTGGGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGc < 2:268407/149‑1 (MQ=2) gAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGcc > 7:53447/1‑148 (MQ=34) gcgcAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGAc > 6:289682/1‑149 (MQ=34) gggagggtgGTGACGGGCCTGGATGCCCGTCAAGCGGTGGGCCCCTGGTTGGGGCCTCTCTCTCGGGGGTCCTACATCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGttttt < 1:293671/148‑1 (MQ=2) agggaggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCtt < 5:289682/149‑1 (MQ=255) agcgcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttt > 4:11131/7‑149 (MQ=255) gcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGt < 6:3948/98‑1 (MQ=255) gcttGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGt > 5:3948/4‑101 (MQ=255) gTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGc > 1:9825/1‑148 (MQ=255) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt < 7:110770/141‑1 (MQ=255) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCttttt > 8:110770/1‑141 (MQ=255) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGa > 3:278061/1‑148 (MQ=255) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGa > 1:181579/1‑148 (MQ=255) cATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACaa > 3:88009/1‑148 (MQ=255) | ATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAGTGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCT‑‑‑GCAGCGCGGTCCTCC‑‑CGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAA > NC_002947/2256057‑2256342 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |