| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,345,527 | T→C | Y271H (TAT→CAT) | rpoN → | RNA polymerase, sigma 54 (sigma N) factor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,345,527 | 0 | T | C | 97.9% | 163.4 / NA | 48 | Y271H (TAT→CAT) | rpoN | RNA polymerase, sigma 54 (sigma N) factor |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (30/17); minor base G (0/1); total (30/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGTATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGACAGCATTCCGCG > NC_000913/3345449‑3345603 | gtgCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCAt > 2:625050/1‑81 (MQ=255) tgCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATg > 1:625051/1‑81 (MQ=255) tgCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATg < 1:624996/81‑1 (MQ=255) gCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGt > 2:625046/1‑81 (MQ=255) gAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCAt > 1:625049/1‑81 (MQ=255) aaaGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCa < 1:624997/79‑1 (MQ=255) aaaGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCAtt > 2:625083/1‑81 (MQ=255) aagaagCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTc > 2:625054/1‑81 (MQ=255) gaagCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCa < 2:624998/81‑1 (MQ=255) agCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTc > 1:625053/1‑77 (MQ=255) gCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGAt > 1:625059/1‑81 (MQ=255) gTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGAtgtg > 1:625043/1‑81 (MQ=255) gTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGAtgtg > 1:625077/1‑81 (MQ=255) aaTCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCt < 1:625000/80‑1 (MQ=255) aaTCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTg > 2:625063/1‑81 (MQ=255) aaTCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTg < 1:624999/81‑1 (MQ=255) aTCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTgg > 2:625048/1‑81 (MQ=255) tGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGc > 1:625057/1‑81 (MQ=255) aTCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGt < 2:625001/77‑1 (MQ=255) aTCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCg > 2:625052/1‑80 (MQ=255) aTCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCg > 2:625078/1‑80 (MQ=255) gCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGGGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATaa < 1:625002/81‑1 (MQ=255) gCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATaa > 1:625079/1‑81 (MQ=255) tCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACg > 2:625084/1‑81 (MQ=255) cGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACgg > 2:625067/1‑81 (MQ=255) gATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGt > 1:625047/1‑81 (MQ=255) gATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGt > 2:625003/1‑81 (MQ=255) gATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGt > 2:625062/1‑81 (MQ=255) aTCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGt > 2:625065/1‑80 (MQ=255) ccGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCAc > 2:625056/1‑81 (MQ=255) ccccGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGAc < 2:625004/81‑1 (MQ=255) cGGGCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGt < 1:625005/81‑1 (MQ=255) gggCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTa < 2:625006/81‑1 (MQ=255) ggCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACgg > 1:625068/1‑78 (MQ=255) gCAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGa < 2:625007/81‑1 (MQ=255) cAGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGaa < 1:625008/81‑1 (MQ=255) aGTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAAc < 1:625009/81‑1 (MQ=255) gTCGATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACt > 2:625060/1‑81 (MQ=255) gATCCAGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCa < 1:625010/80‑1 (MQ=255) aGACTGGCGAACCTGAGGATGTAATTCCAGATGTGCTGGTGCGTAAGCATAAAGGTCACTGGACGGTAGAACTCAACAGTg < 2:625012/81‑1 (MQ=255) aGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGt > 1:625069/1‑80 (MQ=255) aGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGt > 2:625073/1‑80 (MQ=255) aGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTg < 1:625011/81‑1 (MQ=255) aGACTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTg < 2:625013/81‑1 (MQ=255) aCTGGCGAACCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGAc < 2:625014/81‑1 (MQ=255) aCCTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGACAGCATTcc > 1:625075/1‑81 (MQ=255) cTGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGACAGCATTCcgc > 2:625070/1‑81 (MQ=255) tGAGCATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGACAGCATTCcgcg > 1:625061/1‑81 (MQ=255) | GTGCTGAAAGAAGCCGTCAATCTGATCCAGTCGCTCGATCCGCGCCCCGGGCAGTCGATCCAGACTGGCGAACCTGAGTATGTCATTCCAGATGTGCTGGTGCGTAAGCATAACGGTCACTGGACGGTAGAACTCAACAGTGACAGCATTCCGCG > NC_000913/3345449‑3345603 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |