| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,090,124 | C→G | 16.1% | intergenic (+78/‑68) | 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,124 | 0 | C | G | 16.1% | 116.4 / 15.4 | 44 | intergenic (+78/‑68) | hisL/hisG | his operon leader peptide/ATP phosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (18/18); new base G (2/5); total (20/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.20e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01 | |||||||||||
AAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGC > NC_000913/2089976‑2090270 | aaaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCTGAAGATc < 4:188110/149‑1 (MQ=255)aaaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATc < 3:360251/149‑1 (MQ=255) 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 > 2:336032/1‑125 (MQ=255) tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg < 1:336032/125‑1 (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) | AAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGC > NC_000913/2089976‑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: ‑ |