| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,267,126 | G→A | V5I (GTT→ATT) | tyrB → | tyrosine aminotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,267,126 | 0 | G | A | 97.0% | 110.7 / ‑4.0 | 33 | V5I (GTT→ATT) | tyrB | tyrosine aminotransferase |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (17/15); total (18/15) | |||||||||||
| 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.83e-01 | |||||||||||
ATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAA > NC_000913/4267050‑4267196 | aTCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGa > 1:853636/1‑81 (MQ=255) cTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTTACg > 2:853632/1‑81 (MQ=255) gCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTa > 2:853639/1‑81 (MQ=255) cttcCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTg > 2:853628/1‑81 (MQ=255) ccGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTgg < 2:853588/79‑1 (MQ=255) cGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGa < 2:853589/81‑1 (MQ=255) tttATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGAccc < 1:853590/81‑1 (MQ=255) ttATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGAccc > 1:853629/1‑80 (MQ=255) tATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGa > 1:853631/1‑81 (MQ=255) ttGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGAt < 1:853591/80‑1 (MQ=255) tgtTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCtt < 1:853593/81‑1 (MQ=255) tgtTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCt < 1:853594/80‑1 (MQ=255) tgtTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGACTCtt < 1:853592/81‑1 (MQ=255) ttttAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTAc < 1:853595/81‑1 (MQ=255) ttttAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTAc < 1:853596/81‑1 (MQ=255) aaCCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCt > 2:853654/1‑80 (MQ=255) aaCCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCt > 2:853641/1‑80 (MQ=255) caccTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATg > 2:853649/1‑81 (MQ=255) cccGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCg < 1:853597/80‑1 (MQ=255) cGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGttt > 2:853652/1‑81 (MQ=255) cGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGttt < 1:853598/81‑1 (MQ=255) tAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGtttt > 1:853644/1‑80 (MQ=255) tAAACCTGGAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGtttt > 1:853642/1‑80 (MQ=255) ggAGAACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAagaaga > 2:853651/1‑81 (MQ=255) gagaACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGAc > 2:853635/1‑81 (MQ=255) gaACCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGAccc < 1:853599/81‑1 (MQ=255) aaCCATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCt > 2:853637/1‑81 (MQ=255) cATCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCg > 1:853655/1‑80 (MQ=255) aTCGCGTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCg > 2:853647/1‑79 (MQ=255) gcgTGTTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCg < 1:853600/81‑1 (MQ=255) gtgtTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGa > 2:853638/1‑80 (MQ=255) gtgtTTCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGAc < 2:853601/81‑1 (MQ=255) tttCAAAAAATTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACaa < 2:853602/80‑1 (MQ=255) | ATCTTCGCCTTCTTCCGGTTTATTGTGTTTTAACCACCTGCCCGTAAACCTGGAGAACCATCGCGTGTTTCAAAAAGTTGACGCCTACGCTGGCGACCCGATTCTTACGCTTATGGAGCGTTTTAAAGAAGACCCTCGCAGCGACAA > NC_000913/4267050‑4267196 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |