Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,564,097 | (A)7→6 | intergenic (‑92/+37) | yzgL ← / ← glgP | putative uncharacterized protein YzgL/glycogen phosphorylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,564,091 | 0 | A | . | 76.2% | 33.3 / 1.8 | 22 | intergenic (‑86/+43) | yzgL/glgP | putative uncharacterized protein YzgL/glycogen phosphorylase |
Reads supporting (aligned to +/- strand): ref base A (3/2); new base . (5/11); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.25e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01 |
CAGTATGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTC‑TTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑AAAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAC > NC_000913/3563943‑3564450 | agtaaTGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTC‑TTTCCATTC‑TTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGa < 1:1465168/201‑1 (MQ=255) gaggttCAAAATCCATTTTCTGTTTTGCATTATTCTTTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑‑AAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑ttt < 1:1465171/200‑1 (MQ=255) ttCTGTTTTGCATTATTC‑TTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACACTTTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGAC‑TTATCCTGAC‑ATCGACATAGCTg > 1:1465174/1‑305 (MQ=255) cTGTTTTGCATTATTC‑TTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCG‑‑‑‑‑AAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑ATTTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTagat < 1:1465172/269‑2 (MQ=255) tattC‑TTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑‑‑AAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGa < 1:1465173/276‑1 (MQ=255) tCCATTC‑TTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCCTTTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCtt < 1:566142/94‑1 (MQ=255) ttcttttGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑‑TTTCCG‑AATGATGGCGG‑‑AAAAAAAC‑GGGACCCTTTGGCCCCGTTCtattt < 1:566141/101‑1 (MQ=255) cTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑‑TTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCACTCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATcgc > 1:566144/1‑341 (MQ=255) cTTTTTGAATGGTGAAATTATACTCCCCGAGT‑CCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGc > 1:1465175/1‑138 (MQ=255) ggTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGACG‑‑GAAAAAAC‑GGGACCCTTTGGCCCCGTTCtaatt < 1:566143/95‑1 (MQ=255) cccGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑‑AAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCg < 1:1465176/287‑1 (MQ=255) ccccTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑‑AAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCagaagag < 1:1465177/148‑1 (MQ=255) cccTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑AAAAAAAC‑GGGACCCTTTGGCCCCGt > 1:566146/1‑65 (MQ=255) tGCCCC‑TTCTGGACACTTTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑‑TTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGAC‑ATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCa > 1:1465182/1‑293 (MQ=255) tGGACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑‑TTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑t > 1:1465184/1‑163 (MQ=255) tGGACAC‑TTTTCCGAAATGAT‑‑‑GG‑‑AGAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATc < 1:1465178/282‑1 (MQ=255) ggACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTg < 1:1465179/171‑1 (MQ=255) gACAC‑TTTTCCGAAATGATGGCGG‑‑‑AAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTTAA‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGc < 1:1465180/141‑1 (MQ=255) ttttCCGAAATGATGGCGG‑‑‑‑AAAAAC‑GGGACCCTTTGG‑CCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTCTTTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTg < 1:1465181/210‑1 (MQ=255) tttCCGAAATGATGGCGGAAAAAAAAACGGGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTAC‑ATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑‑TTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAc > 1:566148/1‑379 (MQ=255) tttCCGAAATGATGGCGGA‑AAAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTAC‑ATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACttt > 1:566147/1‑205 (MQ=255) atgatgGCGG‑‑‑‑AAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCg < 1:1465183/178‑1 (MQ=255) | CAGTATGGGGATGCCGAAAAAAGTCATTACCAGCGAGGTCAAAATCCATTTTCTGTTTTGCATTATTC‑TTTCCATTCTTTTTGAATGGTGAAATTATACTCCCCGAGTCCCCTTGCCCC‑TTCTGGACAC‑TTTTCCGAAATGATGGCGG‑‑AAAAAAAC‑GGGACCCTTTGGCCCCGTTCTATTT‑A‑TTGGTGAACTTACAATCTCACCGGATCGATATGCCAGATATGATCGGCGTACTC‑TTTGATAGTACGGTCAGAAGAGAAGTAGCCCATATTGGCAATGTTCAGCA‑TCGCTTTTGCGGTCCACTCTTCCTGAAGCTCGTAGAGTTCATCGACTTTATCCTGACAATCGACATAGCTGCGATAATCCGCCAGTACCTGGTAGTGATCGCCGAAGTTGATCAGCGAATCAACCAGATCGCGATAGCGACCCGGATCTTCCGGACTGAATACACCGCTGCCGATTTGCGTCAGCACCTGATGCAGCTCCTCATCTTTCTCGTAGTATTCAC > NC_000913/3563943‑3564450 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |