Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,090,126 | C→T | 18.6% | intergenic (+80/‑66) | hisL → / → hisG | his operon leader peptide/ATP phosphoribosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,090,126 | 0 | C | T | 18.6% | 93.1 / 16.5 | 43 | intergenic (+80/‑66) | hisL/hisG | his operon leader peptide/ATP phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base C (18/17); new base T (2/6); total (20/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01 |
CATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGC > NC_000913/2089979‑2090270 | cATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCAc < 1:201380/148‑1 (MQ=255) caGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg < 3:217256/149‑1 (MQ=255) agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCggg < 4:341959/149‑1 (MQ=255) agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCgg < 2:62035/148‑1 (MQ=255) gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCGGGGGc < 4:601715/149‑1 (MQ=255) acGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg < 3:134539/135‑1 (MQ=255) acGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg > 4:134539/1‑135 (MQ=255) ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGAGCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGtt > 4:590525/1‑149 (MQ=255) aaTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt > 3:201252/1‑149 (MQ=255) tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg < 1:336032/125‑1 (MQ=255) tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg > 2:336032/1‑125 (MQ=255) tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa > 4:623194/1‑149 (MQ=255) ttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa > 2:749586/1‑148 (MQ=255) acacCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGt > 2:125858/1‑148 (MQ=255) caccacCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGtt < 1:125858/148‑1 (MQ=255) ccaccaTCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa > 1:423517/1‑148 (MQ=255) catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga < 3:623194/149‑1 (MQ=255) atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga < 1:498591/148‑1 (MQ=255) atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg < 1:593818/101‑1 (MQ=255) atcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg > 2:593818/1‑101 (MQ=255) tcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg < 4:201252/149‑1 (MQ=255) atcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGaa > 3:252722/1‑149 (MQ=255) atcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGa < 1:135091/148‑1 (MQ=255) tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCAATAAAAAAGCCCCCGGAAGGTGATCTTCCggg < 4:439972/92‑1 (MQ=21) tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCAATAAAAAAGCCCCCGGAAGGTGATCTTCCggg > 3:439972/1‑92 (MQ=25) tAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa > 4:662684/1‑149 (MQ=255) aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATgacagac > 1:659691/1‑148 (MQ=255) gCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGacaa > 1:508249/1‑148 (MQ=255) aTTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTa < 3:590525/149‑1 (MQ=255) ttCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg < 3:202507/98‑1 (MQ=255) ttCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAgg > 4:202507/1‑98 (MQ=255) cTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaa < 2:423517/148‑1 (MQ=255) aTGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa < 3:662684/149‑1 (MQ=255) aTGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTa < 2:508249/148‑1 (MQ=255) tGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTaa > 1:555819/1‑148 (MQ=255) tGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAg < 4:252722/149‑1 (MQ=255) gAACGC‑ATGAGAAAGTCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAg < 1:749586/148‑1 (MQ=255) gAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga < 4:246167/99‑1 (MQ=255) gAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga > 3:246167/1‑99 (MQ=255) aaCGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGt > 2:15741/1‑148 (MQ=255) tGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTCTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGAc > 1:545654/1‑148 (MQ=255) agaAAGCCCCCGGAAGATAACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTc < 2:555819/148‑1 (MQ=255) ggAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGc < 1:15741/148‑1 (MQ=255) tCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgc > 2:42425/1‑148 (MQ=255) | CATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGC > NC_000913/2089979‑2090270 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |