Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,315,884 | 0 | C | T | 33.3% | 59.4 / 11.1 | 37 | S51N (AGC→AAC) | yciG | stress‑induced bacterial acidophilic repeat motifs‑containing protein YciG |
Reads supporting (aligned to +/- strand): ref base C (18/6); new base T (0/12); total (19/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.97e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAGTGTTGTGTGTCGACACTCATTAAAATTAGTCGCTAATGAGAATTAGTCAA > NC_000913/1315639‑1316126 | aaaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTttgt > 1:125559/1‑251 (MQ=255) tAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt > 2:299907/1‑251 (MQ=255) naTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACAGCCttttt > 2:261297/2‑251 (MQ=255) tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTAc < 2:123772/250‑1 (MQ=255) ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt < 2:355908/250‑1 (MQ=255) gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc < 1:237476/251‑1 (MQ=255) ctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGcgcg < 2:415451/249‑1 (MQ=255) gCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCg > 1:12780/1‑251 (MQ=255) gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGCATAAATCCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCCCTGCGAGCCATACGGCCCGCTGATGCAATTATCTTTATCAGCTGATATTAATAATGAACGG‑TGATT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGCCCGCTTTTTTTACCCGCTTCAGATGCGCGTTGCGGATc < 2:307902/251‑1 (MQ=255) ngTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa > 2:242437/2‑250 (MQ=255) tAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑Accgccg > 1:495660/1‑251 (MQ=255) gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCt > 2:450096/1‑251 (MQ=255) gTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGGGCCATACGGCCCGCAGATGCATTTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATTATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCtt < 1:302612/249‑1 (MQ=255) aaaaGCAGGTAAATAAATACTTCTTCATGGGTCTTCATATTCATTTTTTTCCCCATTAAAA‑CCCCTGCGAGCAATCCGCCCCGCAGATGAAATTATCTTTACCAGCTGAT‑TTAAATATTAAGGGGTGAAT‑CAGGTTTGGCCTGATTTTTTCCCCCCGTTTTGTTGACCGCCTTTTTTCCCCGCTTCAGTTGCGCGTTCGGGATCTTTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGt < 1:188852/250‑1 (MQ=255) ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACg > 1:329176/1‑251 (MQ=255) aCATCTTCAATGATCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCActgt > 1:394277/1‑248 (MQ=255) cATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTc > 2:446677/1‑251 (MQ=255) tCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt > 2:181216/1‑251 (MQ=255) aagtttCTTAATATTCATTTTTTTCTCCA‑TGAAAACCACTGCGACCCATCCGCCCCGTATATGCATTTGTCTTTATCGTCTGATATTATTAATTACCGG‑TAATT‑CAGTATTTCCTTTTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTCACCCGCTTCAGATGCGCGTGGGGGATCATTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg < 2:41756/246‑1 (MQ=255) aatattTTTCTCCAGTAAAA‑TCACTGCGACCCAAACGCCCCGCAGATGCAATTATCTTTATCAGCTGATTTTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTCCCACCGTTTTGTTGCCCGCCTTTTTTACCCGCTTAAGAGGCGCGTTGCGGATTATTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAacca < 1:446677/250‑1 (MQ=255) ttttCTCCAGTAAAA‑TCACTGCGAGCCATCCGCCCCGCAGATGCATTTATCTTTATCAGCTGATATTAATAATTAACGG‑TAAAT‑CAGGTTTGGCCTGATTTATTCCCCCCGTTTTGTGGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTAc < 1:243934/196‑1 (MQ=255) ntttCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTAc > 2:243934/2‑196 (MQ=255) ccAGTGAAA‑TCCCTGCGAGCCATACGGCCCGCAGATGAATTTATTTTTATCAGTTGATATTAATAATTAACGG‑TAAAC‑CAGGTTTTGCCTGATTTTTTACCCCCGTTTTTTGGACCGCTTTTTTTACCCGTTCCAGATGCGGGTTGGGGATCTTTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTc < 1:307228/251‑1 (MQ=255) aaa‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATTTTCGGCCAt > 2:226920/1‑251 (MQ=255) gCGAGCAATCCGGCCCGCATATGAAATTATCTTTATCAGCTGATATTAATAATAAACGG‑TAAAT‑CAGATTTTGCCTGTTTTATTCCCACCGTTTTGTGGCCCCCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTAtt < 1:16115/249‑1 (MQ=255) cGAGCCATACGGCCCGCAGATGCAATTATTTTTATCAGTTGATATTAATAATTAGCGG‑GAAAT‑CAGTTTTTGCCTTATTTTTTCCCCCCTTTTTTTTGCCCCCTTTTTTTCCCCGCTCCAGATGCGGGTTGCGGTTCATTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTtct < 1:199861/251‑1 (MQ=255) aatccGGCCCGAAGATGAAATTTTCTTTATCAGCTGATATTAATAATTACCGG‑GAAAT‑CAGGTTTTGCCTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCCCTTCAGTTGCGCGTTGCGGATCTTTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt < 1:450096/247‑1 (MQ=255) cccgaagatcaatttATTTTTACCAGCTAATATTAATAATTAACGG‑TAAAT‑CGGGTTTTGCCTGATTTATTACCCCCGTTTTGTGGCCCGCCTTTTTTCCCCGCTTCAGAGGCGCGTTGCGGATCTTTTTAAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa < 1:208704/197‑1 (MQ=255) cccGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa > 2:208704/1‑210 (MQ=255) cGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGCTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCt > 1:359005/1‑157 (MQ=255) aTGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGtt > 2:457774/1‑251 (MQ=255) atttatttttataggttGATATAAATAATTAA‑GG‑TGAATTCGGTTTTTGCGTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGCTGCAGATGCGCGTTGCGGAT‑ATTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCAt < 2:513229/235‑1 (MQ=255) ttATCAGCTAATATTAAAAATTACCGG‑TGATT‑CAGTATTTGCCTGTTTTATTCCCCCCGTTTTGTTGACCGCTTTTTTTCCCCGTTTAAGGTGCGGGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAata < 1:288596/248‑1 (MQ=255) ttAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGATTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGGTCATTTTTAAACTT‑ACCGCCGCTATGCTGACAGCCTTTACGGCCTGCGTCGGATGCATTCTCACGGACTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCAt > 1:522253/1‑230 (MQ=255) tactaatTAACGG‑TGACT‑CCGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTCACATT‑ACCGCCGCTATGCTGCCCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAATCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCAGTGACTTGTTTCATTAatgtcagttttgttggtcg > 2:411264/4‑233 (MQ=255) agtgatttgcattatgttgtcccaccGTTTTGTTGCCCGCCTTTTTTTCCGGCGTCAGATGCGGTTTGCGATTTTTTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAGTGTGGTGTGTCGACACTCATTAAAATTAGt < 2:135737/229‑1 (MQ=255) accaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAGTGTTGTGTGTCGACACTCATTAAAATTAGTCGCTAATGAGAATTAGTCaa > 2:488456/1‑250 (MQ=255) | AAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAA‑TCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGG‑TGAAT‑CAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATT‑ACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATGAGTGTTGTGTGTCGACACTCATTAAAATTAGTCGCTAATGAGAATTAGTCAA > NC_000913/1315639‑1316126 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |