Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,090,860 | G→A | A223A (GCG→GCA) | hisG → | ATP phosphoribosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,090,860 | 0 | G | A | 92.9% | 124.1 / 0.5 | 42 | A223A (GCG→GCA) | hisG | ATP phosphoribosyltransferase |
Reads supporting (aligned to +/- strand): ref base G (3/0); new base A (21/18); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 |
TGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGC > NC_000913/2090771‑2090954 | tgcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaa > 1:704673/1‑100 (MQ=255) cctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaaaa < 1:1100461/100‑1 (MQ=255) cctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaaaa < 1:518362/100‑1 (MQ=255) aTTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAACCAAAATAca > 2:392115/1‑100 (MQ=255) cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCatgatg > 2:1052303/1‑100 (MQ=255) aTGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGAGTGATCCAGGCACGCGAATCAAAATACATCATGATgca > 2:1103481/1‑100 (MQ=255) ctaTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCacc < 2:908106/98‑1 (MQ=255) aaTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccg < 1:135914/100‑1 (MQ=255) aTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccag > 1:908114/1‑98 (MQ=255) tGGAAGCATCCAAACACCAACTGATCGACAAACTGCTGACCCGTAATCACGGTGTGACACACACACCCGACCCAAAATACATCATGACGCACGCaccga > 2:997080/1‑99 (MQ=255) tGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccgac > 2:513449/1‑100 (MQ=255) aagaaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCAccgaccga < 1:431301/100‑1 (MQ=255) gaaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAAc > 1:693485/1‑100 (MQ=255) aTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACg > 1:424371/1‑99 (MQ=255) tCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTc < 1:992259/100‑1 (MQ=255) tCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTc < 1:315772/100‑1 (MQ=255) cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg < 2:704669/100‑1 (MQ=255) aaaCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTgg > 2:783130/1‑100 (MQ=255) aaaCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATTAAAATACATCATGATGCACGCACCGACCGAACGt > 2:1038209/1‑96 (MQ=255) cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg > 1:793346/1‑100 (MQ=255) cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg < 1:138485/100‑1 (MQ=255) tGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCAt > 1:708542/1‑100 (MQ=255) tCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGc < 1:706897/100‑1 (MQ=255) cAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctg > 2:281268/1‑100 (MQ=255) aaCTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgct > 1:896441/1‑100 (MQ=255) aCTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgctg < 1:539999/100‑1 (MQ=255) gctgACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCa > 1:438031/1‑100 (MQ=255) gACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGt > 1:947046/1‑100 (MQ=255) cccGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGc > 1:384068/1‑100 (MQ=255) ccGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGcc < 1:466951/100‑1 (MQ=255) gTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGa < 2:1193555/100‑1 (MQ=255) aTTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAAc > 2:161170/1‑100 (MQ=255) aTTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACAGAACGGCTGGATGGAGTCATCGCCCTGCTGCCAGGTGCCGAAc > 2:586030/1‑100 (MQ=255) cAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGcc < 1:970507/100‑1 (MQ=255) aGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc > 2:497084/1‑100 (MQ=255) aGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc < 2:845222/100‑1 (MQ=255) gtgtGATCCAGGCACGCGAATCTCAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAAc < 1:1103491/100‑1 (MQ=255) tgATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAt > 1:718788/1‑100 (MQ=255) gATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAtt > 1:847815/1‑100 (MQ=255) aTCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTc < 1:189045/100‑1 (MQ=255) ccAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTg < 1:721011/100‑1 (MQ=255) cAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGc > 1:62667/1‑100 (MQ=255) | TGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGC > NC_000913/2090771‑2090954 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |