Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 837,048 | T→C | 17.6% | F77F (TTT→TTC) | 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,048 | 0 | T | C | 17.6% | 111.4 / 17.6 | 51 | F77F (TTT→TTC) | narJ | molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1 |
Reads supporting (aligned to +/- strand): ref base T (22/20); new base C (4/5); total (26/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
CCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAG > minE/836980‑837117 | ccGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg > 1:3041469/1‑70 (MQ=255) ccGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg > 1:2394300/1‑70 (MQ=255) cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg > 1:1366829/1‑69 (MQ=255) cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa < 1:3822977/71‑1 (MQ=255) cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa > 1:1587823/1‑71 (MQ=255) tCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt > 1:1078563/1‑71 (MQ=255) gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt > 1:332109/1‑71 (MQ=255) gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg < 1:1482210/70‑1 (MQ=255) gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt < 1:523840/69‑1 (MQ=255) gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt < 1:2832612/69‑1 (MQ=255) gCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa > 1:2484570/1‑71 (MQ=255) cccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa > 1:398854/1‑62 (MQ=255) ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc > 1:850909/1‑70 (MQ=255) ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg < 1:953802/71‑1 (MQ=255) ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg < 1:3080495/71‑1 (MQ=255) gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGa < 1:2542945/71‑1 (MQ=255) 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:1256061/71‑1 (MQ=255) gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACcg < 1:2717555/71‑1 (MQ=255) aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa > 1:1874213/1‑45 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGcc < 1:1099489/70‑1 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGc > 1:1050852/1‑69 (MQ=255) ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCa < 1:1531943/71‑1 (MQ=255) cGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGa < 1:3029579/51‑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:1465211/71‑1 (MQ=255) cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc > 1:2221230/1‑71 (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) ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag < 1:2047565/68‑1 (MQ=255) tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag > 1:2668388/1‑70 (MQ=255) tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGAGCAGTACAagcag > 1:326620/1‑70 (MQ=255) | CCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAG > minE/836980‑837117 |
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: ‑ |