| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,339,315 | C→T | intergenic (‑155/‑15) | ribA ← / → pgpB | GTP cyclohydrolase 2/phosphatidylglycerophosphatase B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,339,315 | 0 | C | T | 95.5% | 73.0 / ‑3.6 | 23 | intergenic (‑155/‑15) | ribA/pgpB | GTP cyclohydrolase 2/phosphatidylglycerophosphatase B |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (9/12); total (9/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01 | |||||||||||
CATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCT > NC_000913/1339101‑1339551 | caTCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTa < 2:194785/251‑1 (MQ=255) aaTCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTg > 1:221409/1‑251 (MQ=255) cATGGGGTTGGCAGTTTGGCGTCTGCCACTCGTTTAAGCTGCAGGAAATTCTCCAGATAATGCTGGTTCTGTAGTGGCTTATTTTGCCATAACGAGAAGGGGTGCACCTAATCACGACGCGTCGATCGGTCACGGAATGGCACAAATCTGTCAATTTTTCCGGGAACGGGCGTTTTCAGTTATTATTGTGGGACTTAGCAAAAAGGAGAGGCCATGCTTTCGATTGCCAGACGTACCGCAGTGGGAGCTGc < 1:182866/251‑1 (MQ=255) ggCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGc > 2:216545/1‑251 (MQ=255) gAAATTCTCCAGATAATGCTGGTTCTTTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCtggctggc < 1:218545/251‑1 (MQ=255) aTTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGAGTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACGAAGTTGGCTACTAAAAg < 2:121849/251‑1 (MQ=255) tggCCATAACGAGAAGGGGTTCACCTA‑TCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCgg < 1:84933/248‑1 (MQ=255) gTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATTATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGTCTACTAAAAGCGGCTTTTTGGGTTACTGaaa < 1:20651/251‑1 (MQ=255) gTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGTCTACTAAAAGCGGCTTTTTGGGTTACTGaaa < 1:162748/251‑1 (MQ=255) aTCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCAccc > 2:334411/1‑251 (MQ=255) gacgCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGTGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTATGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGccc > 1:176847/1‑251 (MQ=255) aTCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTAt > 2:257080/1‑156 (MQ=255) aTCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTAt < 1:257080/156‑1 (MQ=255) cACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTac > 1:226442/1‑250 (MQ=255) acaAATCTTTCAAGTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATGAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAGCTGTCACCCAGCCCTGGGGCGTCATTACACATTTGAtttt < 2:247945/251‑1 (MQ=255) ttCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTtct < 1:356323/250‑1 (MQ=255) ggAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTCAAAAGGAGGGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGACAGTAGCCGTATGGAATTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGgtg > 1:425815/1‑251 (MQ=255) gAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTACAAAGGAGAGGCCCTGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTctggct > 1:386875/1‑129 (MQ=255) gAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGGGGGGGCTGCACTATGGCTTGTCATGCCAGTAGCCGTATGGATTTctggct < 2:386875/129‑1 (MQ=255) ggCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCg < 1:225634/251‑1 (MQ=255) aTGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAg < 2:221409/251‑1 (MQ=255) aTGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTAGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCCTGCCAGTAGACGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGATTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTCCACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTa > 1:212034/1‑249 (MQ=255) aTTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCt < 1:334411/251‑1 (MQ=255) | CATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCT > NC_000913/1339101‑1339551 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |