| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 837,066 | A→G | 23.4% | E83E (GAA→GAG) | narJ → | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 837,066 | 0 | A | G | 23.4% | 80.1 / 24.6 | 47 | E83E (GAA→GAG) | narJ | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
| Reads supporting (aligned to +/- strand): ref base A (15/21); new base G (5/6); total (20/27) | |||||||||||
| 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.90e-01 | |||||||||||
CGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACA > minE/836996‑837136 | cgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa < 1:3502846/71‑1 (MQ=255) aCAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc > 1:859730/1‑71 (MQ=255) aGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGaa < 1:3395378/63‑1 (MQ=255) agAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc > 1:3444475/2‑71 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc > 1:1583246/1‑70 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg > 1:2089670/1‑71 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg < 1:2717555/71‑1 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg < 1:1256061/71‑1 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc < 1:1099489/70‑1 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc > 1:1050852/1‑69 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa < 1:1531943/71‑1 (MQ=255) aCCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc > 1:2398564/1‑70 (MQ=255) gTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg < 1:2108400/71‑1 (MQ=255) gggcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt < 1:874864/68‑1 (MQ=255) cgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGa < 1:2391791/71‑1 (MQ=255) accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct > 1:3603164/2‑71 (MQ=255) accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct > 1:1722484/2‑71 (MQ=255) ccACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc < 1:3107988/58‑1 (MQ=255) cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc > 1:1785592/1‑71 (MQ=255) cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc > 1:2221230/1‑71 (MQ=255) cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc > 1:2225904/1‑71 (MQ=255) cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc < 1:1465211/71‑1 (MQ=255) cGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg < 1:1637635/62‑1 (MQ=255) gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt > 1:1731759/2‑57 (MQ=255) aCTGTTGCTGTTTGAACATGTGCACGGCTAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgctg < 1:2037465/67‑1 (MQ=255) tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggt < 1:2800032/55‑1 (MQ=255) tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc > 1:1267926/1‑60 (MQ=255) ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa < 1:29106/70‑1 (MQ=255) cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg > 1:1804337/1‑48 (MQ=255) cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg > 1:1634932/1‑48 (MQ=255) cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa < 1:3321884/71‑1 (MQ=255) cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg < 1:537710/52‑1 (MQ=255) gTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag > 1:3108591/1‑70 (MQ=255) tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag < 1:3216655/69‑1 (MQ=255) ttcgAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa > 1:3425150/4‑68 (MQ=255) ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag < 1:2047565/68‑1 (MQ=255) tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag > 1:2668388/1‑70 (MQ=255) tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGAGCAGTACAagcag > 1:326620/1‑70 (MQ=255) gAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcagc < 1:2181591/70‑1 (MQ=255) cATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACgg < 1:930769/71‑1 (MQ=255) tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt < 1:2562968/71‑1 (MQ=255) tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt > 1:1476639/1‑71 (MQ=255) aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa < 1:277858/37‑1 (MQ=255) aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa < 1:1546293/37‑1 (MQ=255) aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCa < 1:2165574/70‑1 (MQ=255) ccGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGAAAAGGTCGGCCTGCAgct < 1:1233270/69‑3 (MQ=255) cGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaa < 1:2922964/71‑1 (MQ=255) gAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAg < 1:145029/66‑1 (MQ=255) gAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa < 1:3263003/71‑1 (MQ=255) aaTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa < 1:1018712/70‑1 (MQ=255) aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagc > 1:2055620/1‑50 (MQ=255) aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa < 1:406679/71‑1 (MQ=255) aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa < 1:2831789/71‑1 (MQ=255) | CGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACA > minE/836996‑837136 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |