Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,499,480 | 2 bp→AC | intergenic (+57/‑106) | trpS → / → zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,499,480 | 0 | C | A | 97.6% | 108.1 / ‑4.9 | 41 | intergenic (+57/‑107) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base A (18/22); total (19/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
* | NC_002947 | 1,499,481 | 0 | A | C | 100.0% | 109.6 / NA | 41 | intergenic (+58/‑106) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (0/0); new base C (19/22); total (19/22) |
CCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑‑TCGGGCACG‑CCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC > NC_002947/1499338‑1499608 || ccGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGc < 2:316208/148‑1 (MQ=255) cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCtt > 8:87112/1‑149 (MQ=255) ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg < 7:87112/149‑1 (MQ=255) aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGt‑tc > 6:367781/1‑148 (MQ=255) gTTCGACAACGCTGAAGCCAGTGATGGGGCCATTCAGAGCCTGAACGAAGCCTTGGCACAACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCCGCGTTGCCGGTGt‑tgc > 1:11437/1‑147 (MQ=37) cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGt‑tcgcgg > 8:359346/1‑144 (MQ=255) gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCATCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTcc > 1:74017/1‑149 (MQ=25) gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCCCAA‑CAGTCGCCGCGGTGCAAGGCCCCCGCGTTGCCGGTGT‑TC‑GCGGGCAcc‑cccccgcc > 5:289016/1‑143 (MQ=255) cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TA‑GCGGGCACGCCCGCTCCca > 2:100375/1‑148 (MQ=25) gTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCcaca < 5:1435/149‑1 (MQ=21) tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGAc < 1:308303/149‑1 (MQ=21) aaTTCATAGCCTTAACCAAGCCTTGGAACCACACCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCTAAA‑AAGGCGCCGCGGTGCAGGGCACCGGCTTGGCCGGTTTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGc < 2:11437/148‑1 (MQ=18) gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGc > 3:249591/1‑104 (MQ=255) gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGc < 4:249591/104‑1 (MQ=255) ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCg > 8:89542/1‑149 (MQ=21) gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgc > 2:50548/1‑103 (MQ=25) gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgc < 1:50548/104‑2 (MQ=25) aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGc < 7:89542/148‑1 (MQ=21) gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCCAA‑CACGCGCCCCCGTGCAATGCACCCGCGTTGGCGGTGT‑TA‑GCGGGCCCGCCCCCACCCCCAGACAGATTTCCCATTAACCGGGcta > 8:15027/1‑146 (MQ=255) aaGCCTTGCCACCACAGCTGGACTTCGGCTCCCAAACGGTGGGAGCGGGCGTCCCCGCGAAA‑CAGGCGCCGGGGTGCATGGCACGGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATc < 7:324099/149‑1 (MQ=21) ttGCCACCACCGCCGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGAc < 3:69026/120‑1 (MQ=14) gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGTCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACa < 1:100375/148‑1 (MQ=25) accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCAGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGa > 7:204099/1‑149 (MQ=21) aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 7:359346/147‑1 (MQ=21) gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGTCTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 1:91240/145‑1 (MQ=21) gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt > 2:91240/1‑145 (MQ=21) cTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAA‑CCCGCGCCGCGGTGCCTGGCACCGGGGTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc < 5:367781/149‑1 (MQ=21) aaCGGTGGGAGCGGGCGGGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc < 3:353454/129‑1 (MQ=18) aCGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACa < 5:71002/120‑1 (MQ=17) gTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TA‑ACGGGCACGCCCGCGCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc > 7:347154/1‑148 (MQ=25) cccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGAACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt < 2:312896/148‑1 (MQ=37) ccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTAt > 6:214450/1‑94 (MQ=17) ccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTAt < 5:214450/94‑1 (MQ=17) cgcgAAA‑CAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt > 7:246076/1‑115 (MQ=21) cgcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt < 8:246076/115‑1 (MQ=21) gcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 4:66978/119‑1 (MQ=21) gcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt > 3:66978/1‑119 (MQ=21) cgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa < 7:15027/149‑1 (MQ=37) aa‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc > 6:303859/1‑149 (MQ=37) gcgGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc > 7:150151/1‑143 (MQ=37) gcgGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc < 8:150151/143‑1 (MQ=37) || CCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑‑TCGGGCACG‑CCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC > NC_002947/1499338‑1499608 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |