Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,369,793 | G→A | 100% | E35E (GAG→GAA) | pspE → | thiosulfate:cyanide sulfurtransferase (rhodanese) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,369,793 | 0 | G | A | 98.5% | 214.6 / ‑4.3 | 67 | E35E (GAG→GAA) | pspE | thiosulfate:cyanide sulfurtransferase (rhodanese) |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (36/30); total (36/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01 |
AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC > NC_000913/1369696‑1369889 | aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACa < 2:984183/100‑1 (MQ=255) aaaGGCTTACTTGCGCTGGCACTGGTGTTTTCACTGCTCGTTTTCGCAGCTGTACACTGGATCAATGTTCGTGTTCCAGAGCAGTATCAGCAAGAacac < 2:586813/99‑1 (MQ=255) ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTc < 2:2268807/100‑1 (MQ=255) cTTACTTGCTCTGGCACTGGTGTTTTCACTGCGCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAg > 1:1393913/1‑100 (MQ=255) ttACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAgg < 1:609550/100‑1 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAggggg < 1:475471/100‑1 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGc < 1:2488976/100‑1 (MQ=255) tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGCACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGcc > 2:1741425/1‑100 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATc < 1:2493171/100‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAata > 1:414373/1‑100 (MQ=255) cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTc > 1:1668744/1‑100 (MQ=255) ggTGTTTTCACTGCCCGGTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcc < 2:2158664/99‑1 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc < 2:338231/100‑1 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc < 2:1257514/100‑1 (MQ=255) tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaa < 2:383077/100‑1 (MQ=255) cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag > 2:1247561/1‑100 (MQ=255) ccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaa > 1:435902/1‑100 (MQ=255) cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa > 1:606713/1‑100 (MQ=255) cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGAACCATCAATATTCCCCTGAAGGAAGTGaaa > 2:2297903/1‑100 (MQ=255) gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAag < 1:154228/100‑1 (MQ=255) tCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg < 2:1473436/100‑1 (MQ=255) tCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg > 1:2241535/1‑100 (MQ=255) cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc > 1:1106867/1‑100 (MQ=255) gccgcTGAACACTGGATCGATGTTCGTGTTCCAGATCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa < 2:399854/100‑1 (MQ=255) ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt > 1:1132787/1‑100 (MQ=255) tGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGcc > 1:1841872/1‑100 (MQ=255) gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa < 1:1832021/100‑1 (MQ=255) gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa < 1:1827790/100‑1 (MQ=255) gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa > 2:441250/1‑100 (MQ=255) aaCACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTTCCCTGAAGGAAGTGAAATAGCGCATTGCCAc < 1:1725207/100‑1 (MQ=255) acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc > 2:2118076/1‑100 (MQ=255) ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt > 1:2008671/1‑100 (MQ=255) ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt > 2:2118571/1‑100 (MQ=255) gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt > 2:1131229/1‑100 (MQ=255) aTGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGa < 1:1586762/100‑1 (MQ=255) tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt > 2:2091894/1‑100 (MQ=255) tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt < 1:1821771/100‑1 (MQ=255) gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa > 1:573557/1‑100 (MQ=255) tCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCCATATTCACCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaa < 2:2154325/99‑1 (MQ=255) tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc < 1:1685/100‑1 (MQ=255) tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc < 2:1942115/100‑1 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGACGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg < 2:1563693/100‑1 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg < 1:167271/100‑1 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg > 1:529453/1‑100 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg < 2:2008219/100‑1 (MQ=255) gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg > 2:2499485/1‑100 (MQ=255) tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac > 2:1593664/1‑100 (MQ=255) ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca > 1:72564/1‑100 (MQ=255) ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca > 2:255928/1‑100 (MQ=255) ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca > 1:209848/1‑100 (MQ=255) cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac > 1:1542587/1‑100 (MQ=255) agagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACAcc > 2:479393/1‑100 (MQ=255) gagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg < 1:1593675/100‑1 (MQ=255) gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg < 1:1035475/100‑1 (MQ=255) gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg > 2:520679/1‑100 (MQ=255) cAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATCCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGa < 1:11861/100‑1 (MQ=255) aGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa > 2:73018/1‑100 (MQ=255) tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg < 1:1131237/100‑1 (MQ=255) aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg > 2:2045109/1‑99 (MQ=255) tCAGCAAGAACACGTTCAGGGGGCCATCAATATTCGCCTGAAGGAAGTGAAAGAGCGCATTGCAACCGCCGTTCCGGATAAAAACGACACAGTGAAAgt > 2:2242626/1‑99 (MQ=255) tCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg < 2:99833/100‑1 (MQ=255) cAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCCGATAAACACGACACCGTGAAAatg > 2:669800/1‑99 (MQ=255) cAGCAAGAACACGTTCAGGGAGCCATCAATATTCCGATGAAGGAAGTGAAAGAGCGCATTGCCACTGCCCATCCGGATAAAACCGACACCATGAAAgtgt > 2:874638/1‑100 (MQ=255) aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa < 1:1240752/100‑1 (MQ=255) aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa > 1:482938/1‑100 (MQ=255) gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt > 2:1685341/1‑100 (MQ=255) cAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAtt > 1:324661/1‑100 (MQ=255) aaGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTg < 2:1706333/100‑1 (MQ=255) agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc > 1:1465891/1‑100 (MQ=255) agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc < 1:2118585/100‑1 (MQ=255) | AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC > NC_000913/1369696‑1369889 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |