Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 766,716 | 0 | G | T | 40.0% | 7.5 / 10.6 | 25 | A245S (GCG→TCG) | mngA | 2‑O‑alpha‑mannosyl‑D‑glycerate specific PTS enzyme II |
Reads supporting (aligned to +/- strand): ref base G (9/6); new base T (10/0); total (19/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.08e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCT > NC_000913/766474‑766966 | tGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAATGC‑AA‑ACCCCTGGATTTTTTCGGCTgg > 1:168121/1‑249 (MQ=255) aaCGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGGTTAACGGCGGAACAACTCCAAAATGC‑AA‑ACCCCTGGTTTTTTTCCGCTGGAGTGGCCCCCAAAGAGAGTGAGGGTTTTTAT‑gggcttccccccgctttactgac > 1:389865/1‑230 (MQ=255) cgccaatctcccGTCTGCGGAAGGGCCTGGCGGGGGGGTGTTGCTCACTCTTATGTGGCGACGGAATATCTGAAAAAAGCCGGACGAAGCCGCGGTTTAAA‑GTTTGAGTTGAAAAAACAAGGCGGAAGCGTTTTGGAAGGGTGTTTATGGGGGGAGCAATTAAATGGTGC‑AG‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGTAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGt < 1:620071/240‑1 (MQ=255) aaCGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAACAACTCAACAATGG‑AA‑CCGCCTGTATTTTTTCGGCTGAAAGAGCCATCAAGGAGAGTGGGGGGTTTAAT‑TGCATTTCCCCGCCTTTCTG‑GCCTGTTGCCCGGCCGAAT‑CGCCAtt > 1:381602/1‑250 (MQ=255) gggTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGC‑AA‑CCGCCTGTATTTTTTCGGCTGGAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCCGAAGCG‑TTGGTcc > 1:210800/1‑251 (MQ=255) tGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTACCGGCATTGAAGGGCGTTTAACGGCGGATTAAATCAATAGTGCCAA‑CC‑CCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGCGAGAGCGTTTTAAT‑GGAATTTCCGCCC‑TTTCAGTGACTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCcaacc > 1:660975/1‑248 (MQ=255) aaaaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAACAAATCAATAAAGC‑AA‑CCGCCCGTATTTTTTCGGGTGAAGCCGCCACCCAGGGGAGGGAGGGGTTTAAA‑GGGATTCCCGCGC‑TTTTAGTGCGTGTTTTCGCACCGATTACGC‑ATGCAGGAGCGGTTGTTCCACCAA‑CGC‑TTACCCTcccgacgaacgaggtgacccgc > 2:547835/1‑230 (MQ=255) aaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGCTCAACTCAACAGTCC‑CA‑CCGCCTGTATTTTTTCGGCTGAAGTCGCCACCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCCTGCCAAAGCG‑TTGATCCCACAAGCGC‑TTACCCTCAAGCGTAGCGGTGAGACGcgtac > 1:343774/1‑251 (MQ=255) ccGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CATGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGC‑CA‑CCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGAGAGTGAGCGGTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATTCAACAAGCGC‑TTACCCTCCAGCGTAGCGATGAGACGcgtacc‑gt > 1:67245/1‑250 (MQ=255) ggACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATCGTGC‑AA‑CCGCCTGTATTTTTTCGGCTGCAGGCGCCCTCAAGGAGAGTGAGCGGTTTTAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCCATT‑CGCCATGCCGGAGCG‑TTGATCCAACAAGCGC‑TTCCCCCCCAGCGGAGCGGTGAGACGcgtacc‑gtac > 1:304539/1‑250 (MQ=255) aCTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGC‑CA‑CCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGGGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGGCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGAt > 1:143615/1‑250 (MQ=255) aaaTGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGGGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGGGAGTGAGCGTTTTTAT‑GGCATTCCCGCGC‑TTTCCGTGCCTGTTGCCGAGCCGATT‑CGCCCTGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCCCAAGCGGAGCGGTGAGACGCGTACC‑GTAAAGCCAGATAcagatcggg > 1:10536/1‑243 (MQ=255) a‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCCACTCAATAGTGC‑AA‑CCGCCCGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTTAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑CTACCCTCAAGCGGAGCGGTGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACgg > 1:323110/1‑251 (MQ=255) gcTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCCATAGTGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCCTTTCCGCGC‑TTTCAGTGCCTGTTGCCGCGCCGATT‑CGCCATGCAGAAGCG‑TTGGTCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTCCC‑GTACAGCAAGATACGCAACCGGGGAAAAGTGGCAAAAAGGAGCTGa > 1:22737/1‑250 (MQ=255) ggCATTGAAGGGCGTTTAACGGCGGATCAACTTAATAGTGC‑AA‑CCGCCTGTATTTTTTCGGCTGAAGGCGCCATCCAGGGGAGGGAGCGGTTTAATTGGC‑TTCCCCCGC‑CTTCCGTGCCTGTTGCCGAGCCGGTT‑CGCCCTGCAGAAGCGGTTG‑TCCAACAAGCGCCTTACC‑TCAACCGGAGCGGTTGGACGCGGACCCGTA‑AGCCAGAGACGCAACCGGGGCAAAGTGTCAAAAAGGAGGTGAAAAAggg > 2:510124/1‑250 (MQ=255) gtgtaggggggatttgagcgcaggATCAACTCAATAGCGC‑AC‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCa < 2:168121/229‑1 (MQ=255) ttGAAGGGCGTTTAACGGCGGCTCAACTCAAAAGTGC‑AA‑ACACCTCTATTTTTGCGGCTGAAGACGCCATCCAGGAGAGGGAGCGCTTTTAT‑GGCATTCCCGCGC‑CTTCAGTGCCTGGTGCCGAGCCGATT‑CTCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGCTGAGACGCGTACC‑GTACAGCAAGAAACGCAAACGGTGAAAAGTGTCAAAAAGGAGCTGAAACAGGCACTg > 1:414153/1‑251 (MQ=255) gggcgaaaaaccaaaaaGTGC‑AATCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATGCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCGGAAACAGGCACTGTTGAGCGGAATctct < 1:570395/235‑1 (MQ=255) cTCAATAGTGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATTAAGGAGAGTGAGCGTATTAAT‑GGCATTCCCGCGA‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGAGCAAAGTGTCAACACGGCGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTAcc > 1:81902/1‑251 (MQ=255) aTAGTGC‑AA‑CAGCCTGTATTTTTGCGGCTGAAGTCGACATCAAGGCGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGCGCCGATT‑CGCCATGCAGACGCG‑TTGATCCAACAAGCGC‑TTACCCACAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACGGTTg > 2:662099/1‑251 (MQ=255) aGTGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGCGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGCGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGAt > 1:613523/1‑251 (MQ=255) a‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCCGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGCGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGCATCCCTTTTGCCGTACCGTTGATTGTcgg > 2:298389/1‑250 (MQ=255) aTTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCa < 1:2620/250‑1 (MQ=255) aTTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCAc < 1:674359/251‑1 (MQ=255) ttGCGGCTGAAGTCGCCATCAAGGAGAGTGTGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCACCCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGt < 2:381602/250‑1 (MQ=255) gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTACT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCt > 1:39650/1‑251 (MQ=255) gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAAc > 2:361407/1‑190 (MQ=255) gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTggtggcggcagcacagaa > 1:38457/1‑232 (MQ=255) gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGAGGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAAc < 1:361407/190‑1 (MQ=255) | TGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAA‑CAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGC‑AA‑CCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAAT‑GGCATTCCCGCGC‑TTTCAGTGCCTGTTGCCGAGCCGATT‑CGCCATGCAGAAGCG‑TTGATCCAACAAGCGC‑TTACCCTCAAGCGTAGCGATGAGACGCGTACC‑GTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCT > NC_000913/766474‑766966 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |