| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | ACB122 | 1,898,174 | A→C | 27.0% | K458T (AAG→ACG) | nasD → | Nitrite reductase [NAD(P)H] |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 1,898,174 | 0 | A | C | 27.0% | 67.8 / 5.0 | 59 | K458T (AAG→ACG) | nasD | Nitrite reductase [NAD(P)H] |
| Reads supporting (aligned to +/- strand): ref base A (18/22); new base C (2/16); total (20/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.63e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAGGAGCTGTACGG > ACB122/1898044‑1898306 | cgggCCCGGGGCGCCGAGGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGGCCGCAGGGGGGCGGGGGGCGCCGCGATCGACAGCGGCGGCAGCGACCGCGCTGCGGGCAAGGGCGGC‑ACCAAGGcc < 6:3133/133‑1 (MQ=255)ccGGCCGTGGGGGCCTTTGCGCGGCCCGACAGCGCGACCATCTGCGCATCCAACACCGTCAGCAAGGGGGCGGGGTGCGCGGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGCCAAGGGCTGC‑ACCACGGcc < 2:53939/135‑1 (MQ=255)ccGGCCCTGGGCGCCGATGCGCGGCCCGACAGCGCGACCATCGGCGCCGGCCACAACGGCAGCAAGGGGGCGGGGTGCGCCGCGAGCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGcc < 2:247604/135‑1 (MQ=255) cGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGCCAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCcg > 8:330023/1‑135 (MQ=255) ggCCCTGGGCGCCGCTGCGCTGCCCGACAGCGCGACCATCGGCACCTGCCACACCGTCAGCAAGGGGGCGGGGTGCGCCGCGAGCGACCGCGGCGGCAGCGACCTCGCTGCGGTCAAGGGCGGC‑CCCAAGGCcgc < 1:121484/135‑1 (MQ=255) gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGCCCATCTGCTCCTGCCCCAACGTCCGCAAGGGGGCGGTGTGCGCCGCGCTCGCCAGCGGCTGCCGCGACCTCGCGGCGGTCCAGGGCTGC‑CCCCAGGCcgcg > 1:356558/1‑135 (MQ=255) gCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGCCCTCGCGGCGGTCAAGGGCTGC‑CCCACGGCcgcg > 8:157338/1‑135 (MQ=255) cTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGccc > 2:46011/1‑135 (MQ=255) gggCGCCGATGCGCGGCCCGACAGCGCGACCTTCTGCGCCTGCCCCAACGGCCGCAAGGGGGCGGGGGGCGCCGCGATCGACAGCGGCGGCAGCGACCTCGCTGCGGGCTAGGGCTGC‑ACCAAGGCCGCGACCgg < 7:22088/135‑1 (MQ=255) cGATGGGCTGCCCGACAGCGCGACCATCTCCGCCCCCCCCAACGGCAGCACGGGGGCGGGGGGCGCCGCTATCGACCGCGGCTGCAGCGACCTCGCGGCGGTCAAGGGCTCC‑ACCACGCCCGCGACCGGCTGCgg < 2:322396/135‑1 (MQ=255) cGATGCGCGGCCCGACAGCGCGACCAGCTGCTCCTGCCACACCGTCAGCAAGGGGGCGGTGTGCGCCGCGAGCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCgg < 7:67424/135‑1 (MQ=255) tGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACACCGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACCGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCggtgg > 1:303412/1‑135 (MQ=255) tGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGCCAGCGGCTGCAGCGACCGCGCTGCGGTCAAGGGCTGC‑ACCAATGCCGCGACCGGCTGCGGTGGTTgcg > 2:264939/1‑135 (MQ=255) cccGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGcgcg < 4:108004/135‑1 (MQ=255) ccGACCGCGCGACATTCGGCTCCGCCCACAACGCCCGCAGGGGGGGGGGGTGCGCCGCGATCGCCAGCGGCGGCAGCGACCGCGCGGCGGCCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCgg < 5:365670/135‑1 (MQ=255) gACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGCCCGGCTGCGGTGGTTGCGCGGCCCTGCTcaa > 8:248825/1‑135 (MQ=255) gACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaa > 1:102040/1‑135 (MQ=255) gACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaa > 2:158323/1‑135 (MQ=255) ccATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTcaag > 7:345486/1‑134 (MQ=255) ccTGCCACACCGTCAGCCAGGGGGCGGTGGGCGCCGCGAGCGACAGCGGCGGCAGCGACCGCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCg < 7:329962/135‑1 (MQ=255) cTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGa < 7:195947/135‑1 (MQ=255) tGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGCCCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGa > 2:62032/1‑134 (MQ=255) acaACGTCAGCAAGGGTGCGGTTTGCGCCGCGATCGACAGCGGCTGAAACGCCCTCGCTGCCGTCAAGGGCTTC‑ACCAAGCCCGCGACAGGCTGCGGAGGTTGCGCGGCCCTGATAAAGCAAGAAATCGAGCAc > 4:9018/1‑134 (MQ=255) acaACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCAc > 8:223302/1‑134 (MQ=255) acaACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCAc > 5:48655/1‑134 (MQ=255) aaGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTca < 6:169129/106‑1 (MQ=255) aaGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTca > 5:169129/1‑106 (MQ=255) gCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCACGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCgcggc > 8:4078/1‑135 (MQ=255) gggggtgCGCCGCGAGCGCCAGCGGCTGCCGCGCCCTCGCTGCGGTCAAGGGCTGC‑ACCAAGCCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCgcggcg < 2:61194/131‑1 (MQ=255) gggggggcgcCGCGATCGACAGCGGCGGCAGCGACCTCGCTGCGGTCAAGGGCGGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCgcggcg < 4:21076/129‑1 (MQ=255) cGGGGTGCGCCGCTATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCCAGGGCTGC‑CCCACGGCCGCGCCCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCgcggcg < 2:261266/135‑1 (MQ=255) gggggcgcGGCGAGCGACTGGGGGTGCCGGCACGTGGCTGGGGACACGGGCGGC‑CCCACGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTg < 4:230025/131‑1 (MQ=255) gggggcgcCGCGATCGACAGCGGCGGCAGCGACCGCGCTGCGGTCAAGGGCTGC‑ACCACGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTg < 1:309979/131‑1 (MQ=255) gggggcgcCGCGAGCGACAGCGGCTGCAGCGACCGCGCTGCGGTCAAGGGCTGC‑ACCACGGCCGCGACCGGCGGCGGTGGGTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTg < 6:60830/131‑1 (MQ=255) gtgCGCCGGGAGGCACAGCGGCTGCTGCTACCTGGCTGCGGTAAAGGGCTCC‑CCCAAGCCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCATGCAAGTCTTCGAGCACGAACTGACTGCCCGCGGCGTGGc < 4:252183/135‑1 (MQ=255) gtgCGCCGCGATGGACCGCGGCTGCAGCTACCGCGCTGCGGGCATGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGc < 4:252177/135‑1 (MQ=255) gggcgcGGCGAGCGACCGCGCGGGCCGCGACCGCGCGGGGGTCACGGGCGTC‑CCCACGGCCGCGACCGGCGGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGc < 2:303412/133‑1 (MQ=255) gggcgcCGCTCGCGACCGCGGCGGCCGCGACCGCGCGGCGGACACGGGCGCC‑ACCACGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCACGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGc < 7:330023/133‑1 (MQ=255) gggcgcCGCGAGCGCCCGCGGGGGCCGCTCCCGCGGGGCGGCCACGGGCGCC‑ACCACGGCCGCGACCGGCGGCGGTGGGTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGc < 1:347793/133‑1 (MQ=255) tgCCCCGCGATCGACCGCGGCGGCTGCGCCCGCGGTGGGGT‑AACGGCCGCAACCACGGCCGCGACGGGCTGGGGTGGTTGCGCGACCCTGCTCAACAAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCg < 2:77079/135‑1 (MQ=255) cggcgagcgaCAGCGGGTGCAGCGACCGCGCTGCGGTCACGGGCTGC‑ACCATGCCCGCGACCGGCTGCGGGGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGCGCACGAACTGACGGCCCGCGGCGTGGCGGTgg < 1:264939/133‑1 (MQ=255) tCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAaga < 7:207620/135‑1 (MQ=255) tCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGACCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAaga > 1:77079/1‑135 (MQ=255) gcGACAGCGGCGGCTGCCACCGCGCTGCGGGCACGGGCGCC‑ACCACGGCCGCGACCGGCGGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAaga < 2:102040/134‑1 (MQ=255) gCAGCTCCCTGGCTGCGGTCATGGGCTGC‑CCCACGGCCGCGACCGGCTGCGGGGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGc < 4:346261/135‑1 (MQ=255) gCAGCGACGTCGCGGCGGTCACGGGCGGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCACGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGc < 2:211533/135‑1 (MQ=255) gCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGc < 2:71257/135‑1 (MQ=255) ccgcgcCCGCGCTGGGGGCACGGGCGGC‑CCCACGCCCGCGACCGGCGGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCa < 2:356558/133‑1 (MQ=255) gACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACttt < 2:216389/135‑1 (MQ=255) cGGTCCCGGGCGGC‑ACCACGGCCGCGACCGGCGGGGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACGGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACAccc < 6:248640/135‑1 (MQ=255) acgggCTGC‑CCCACGCCCGCGACCGGCGGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAg < 1:158323/133‑1 (MQ=255) aaGGGCTGC‑ACCACGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAg < 7:108462/135‑1 (MQ=255) aaGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCa > 4:337308/1‑134 (MQ=255) tggccggC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAACCAAGTCTTCGAGCACTAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAg < 6:349876/128‑1 (MQ=255) aGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAg < 3:337308/134‑1 (MQ=255) gcggC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAGGAGc < 4:361670/132‑1 (MQ=255) cAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACAccc > 4:104367/1‑119 (MQ=255) cAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACAccc < 3:104367/119‑1 (MQ=255) cAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAGGAGCTGTACgg > 6:149962/1‑135 (MQ=255) | CCGGCCCTGGGCGCCGATGCGCTGCCCGACAGCGCGACCATCTGCTCCTGCCACAACGTCAGCAAGGGGGCGGTGTGCGCCGCGATCGACAGCGGCTGCAGCGACCTCGCTGCGGTCAAGGGCTGC‑ACCAAGGCCGCGACCGGCTGCGGTGGTTGCGCGGCCCTGCTCAAGCAAGTCTTCGAGCACGAACTGACGGCCCGCGGCGTGGCGGTGGACAAGAGCCTGTGCGAGCACTTTGCCTACACCCGGCAGGAGCTGTACGG > ACB122/1898044‑1898306 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |