Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,792,768 | A→G | 18.0% | intergenic (‑217/+58) | kbl ← / ← yibB | glycine C‑acetyltransferase/YibB family protein, function unknown |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,792,768 | 0 | A | G | 18.0% | 119.3 / 17.6 | 50 | intergenic (‑217/+58) | kbl/yibB | glycine C‑acetyltransferase/YibB family protein, function unknown |
Reads supporting (aligned to +/- strand): ref base A (21/20); new base G (4/5); total (25/25) | |||||||||||
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.82e-01 |
TTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC > NC_000913/3792629‑3792903 | ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGGAAAAATGATGCTTTGGCATCAtt < 3:546381/148‑1 (MQ=255) aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg < 3:1071683/139‑1 (MQ=255) aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg > 4:1071683/1‑139 (MQ=255) tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt < 2:125992/115‑1 (MQ=255) tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt > 1:125992/1‑115 (MQ=255) tctctcAACGCCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAGTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca > 4:739232/1‑105 (MQ=255) tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt < 3:924451/110‑1 (MQ=255) tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt > 2:1067939/1‑124 (MQ=255) tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt < 1:1067939/124‑1 (MQ=255) tctctcAACACCAACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca < 3:739232/105‑1 (MQ=255) ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt < 2:93771/127‑1 (MQ=255) ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt > 1:93771/1‑127 (MQ=255) acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt < 1:5389/104‑1 (MQ=255) acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt > 2:5389/1‑104 (MQ=255) acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc < 1:211799/100‑1 (MQ=255) acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc > 2:211799/1‑100 (MQ=255) aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt < 3:982843/107‑1 (MQ=255) aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt > 4:982843/1‑107 (MQ=255) aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc < 4:644204/91‑1 (MQ=255) aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc > 3:644204/1‑91 (MQ=255) aaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt > 3:424388/1‑85 (MQ=39) gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg > 3:25986/1‑72 (MQ=255) gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg < 4:25986/72‑1 (MQ=255) gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg < 2:475409/99‑1 (MQ=255) gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg > 1:475409/1‑99 (MQ=255) atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa > 1:946134/1‑102 (MQ=255) atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa < 2:946134/102‑1 (MQ=255) atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg > 4:837206/1‑145 (MQ=255) atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg < 3:837206/145‑1 (MQ=255) taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt < 2:914553/89‑1 (MQ=255) taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt > 1:914553/1‑89 (MQ=255) aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg > 1:304532/1‑130 (MQ=255) aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg < 2:304532/130‑1 (MQ=255) aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa > 1:1013613/1‑91 (MQ=255) aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa < 2:1013613/91‑1 (MQ=255) aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa < 2:940836/90‑1 (MQ=255) aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa > 1:940836/1‑90 (MQ=255) aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt > 2:165574/1‑85 (MQ=255) aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt < 1:165574/85‑1 (MQ=255) aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACTCACTAATGGAGGGAAACTTATTTTCTTGaa > 2:229890/1‑83 (MQ=37) aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa < 1:229890/83‑1 (MQ=38) aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg > 3:802002/1‑111 (MQ=255) aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg < 4:802002/111‑1 (MQ=255) aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAGATc > 4:219191/1‑148 (MQ=255) aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt > 2:431772/1‑78 (MQ=255) aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt < 1:431772/78‑1 (MQ=255) atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt < 1:395811/45‑1 (MQ=255) atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt > 2:395811/1‑45 (MQ=255) gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa > 4:192038/1‑110 (MQ=255) gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa < 3:192038/110‑1 (MQ=255) aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt < 2:783164/69‑1 (MQ=255) aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt > 1:783164/1‑69 (MQ=255) aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc > 3:754307/1‑82 (MQ=255) aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc < 4:754307/82‑1 (MQ=255) aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt > 3:324479/1‑136 (MQ=255) aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTCTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCACTAAAGTGTTACTCCTGAAGCAGCGAAACaatt < 4:324479/136‑3 (MQ=255) | TTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC > NC_000913/3792629‑3792903 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |