Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,090,122 | A→G | 26.3% | intergenic (+76/‑70) | 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,122 | 0 | A | G | 26.3% | 66.0 / 24.2 | 38 | intergenic (+76/‑70) | hisL/hisG | his operon leader peptide/ATP phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base A (14/14); new base G (5/5); total (19/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01 |
ATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATT > NC_000913/2089974‑2090262 | aTAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGa < 4:15062/149‑1 (MQ=255) aaaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATc < 2:783886/149‑1 (MQ=255) aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCt > 2:1842239/1‑149 (MQ=255) acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCg > 2:669987/1‑149 (MQ=255) agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggg < 1:1842239/149‑1 (MQ=255) agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggg < 1:669987/149‑1 (MQ=255) gagaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg < 1:617526/149‑1 (MQ=255) gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggggg > 4:1687677/1‑132 (MQ=255) gTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggggg < 3:1687677/132‑1 (MQ=255) aaTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg < 1:23733/125‑1 (MQ=255) aaTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg > 2:23733/1‑125 (MQ=255) tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa > 1:34472/1‑149 (MQ=255) ccatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCAc > 3:257299/1‑105 (MQ=255) ccatcatcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCAc < 4:257299/105‑1 (MQ=255) tcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg < 3:553839/106‑1 (MQ=255) tcaccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgg > 4:553839/1‑106 (MQ=255) caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtt < 4:1268627/111‑1 (MQ=255) caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtt > 3:1268627/1‑111 (MQ=255) catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCAATAAAAAAGCCCCCGGAAGGTGATCTTCCgg < 4:314686/103‑1 (MQ=25) catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCAATAAAAAAGCCCCCGGAAGGTGATCTTCCgg > 3:314686/1‑103 (MQ=25) tttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGa < 4:1083108/74‑1 (MQ=255) tttCAGGCGATCTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGa > 3:1083108/1‑74 (MQ=255) ggAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggggg > 1:1543055/1‑70 (MQ=255) ggAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCggggg < 2:1543055/70‑1 (MQ=255) ttCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg < 2:656355/61‑1 (MQ=255) ttCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg > 1:656355/1‑61 (MQ=255) cAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg < 4:191210/59‑1 (MQ=255) cAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCgggg > 3:191210/1‑59 (MQ=255) aGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCgg < 2:34472/149‑1 (MQ=255) ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt > 3:60092/1‑98 (MQ=39) ccccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGtt < 4:60092/98‑1 (MQ=255) cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga < 2:295267/82‑1 (MQ=38) cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga > 1:295267/1‑82 (MQ=38) cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa > 4:696860/1‑117 (MQ=255) cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa < 3:696860/117‑1 (MQ=255) cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt > 1:922341/1‑149 (MQ=255) cccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt > 3:1636631/1‑149 (MQ=255) gAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaa < 2:230737/112‑1 (MQ=255) gAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaa > 1:230737/1‑112 (MQ=255) gATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGAc < 2:1549171/131‑1 (MQ=255) gATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGAc > 1:1549171/1‑131 (MQ=255) | ATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGC‑ATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATT > NC_000913/2089974‑2090262 |
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: ‑ |