| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 837,075 | C→T | 20.7% | D86D (GAC→GAT) | 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,075 | 0 | C | T | 20.7% | 106.8 / 22.8 | 53 | D86D (GAC→GAT) | narJ | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
| Reads supporting (aligned to +/- strand): ref base C (15/27); new base T (5/6); total (20/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.28e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
CGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGA > minE/837006‑837143 | agAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc > 1:3444475/2‑71 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc > 1:1583246/1‑70 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg > 1:2089670/1‑71 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg < 1:1256061/71‑1 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg < 1:2717555/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:2225904/1‑71 (MQ=255) cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc > 1:2221230/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:1634932/1‑48 (MQ=255) cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg > 1:1804337/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:1476639/1‑71 (MQ=255) tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt < 1:2562968/71‑1 (MQ=255) aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa < 1:1546293/37‑1 (MQ=255) aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa < 1:277858/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) tCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa > 1:2930074/1‑68 (MQ=255) ccGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAg < 1:3461646/69‑1 (MQ=255) cgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGc > 1:1244298/1‑67 (MQ=255) cgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg < 1:581286/71‑1 (MQ=255) aCCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCATTTAAACAGCCGCGa < 1:2929158/70‑1 (MQ=255) | CGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGA > minE/837006‑837143 |
| 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: ‑ |