Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,985,112 | C→G | 17.9% | intergenic (‑43/+27) | araG ← / ← araF | L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,985,112 | 0 | C | G | 17.9% | 152.2 / 27.2 | 56 | intergenic (‑43/+27) | araG/araF | L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein |
Reads supporting (aligned to +/- strand): ref base C (22/24); new base G (5/5); total (27/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01 |
TGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCC > NC_000913/1984969‑1985251 | tGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc > 1:48465/1‑148 (MQ=255) ccGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga < 4:561716/148‑1 (MQ=255) ccGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga < 3:1008346/148‑1 (MQ=255) tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTc < 2:48465/148‑1 (MQ=255) ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga > 4:888713/1‑123 (MQ=255) ccGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCaga < 3:888713/123‑1 (MQ=255) ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag > 4:599730/1‑111 (MQ=255) ttAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag < 3:599730/111‑1 (MQ=255) aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg < 3:1098052/112‑1 (MQ=255) aaCGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg > 4:1098052/1‑112 (MQ=255) aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg > 3:1073024/1‑103 (MQ=255) aaaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑gg < 4:1073024/103‑1 (MQ=255) gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGTATCACTGGTCGTTCACACCCCCTCTGCATCATGc < 1:858471/73‑1 (MQ=255) gagagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGc > 2:858471/1‑73 (MQ=255) tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa > 2:29202/1‑130 (MQ=255) tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTaa < 1:29202/130‑1 (MQ=255) tAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑ggg > 3:642089/1‑66 (MQ=255) tAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGA‑ggg < 4:642089/66‑1 (MQ=255) gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg > 3:848748/1‑124 (MQ=255) gCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcg < 4:848748/124‑1 (MQ=255) aCGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg > 1:331842/1‑148 (MQ=255) ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTc < 4:1042003/57‑1 (MQ=255) ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTc > 3:1042003/1‑57 (MQ=255) tGAATCACTGGTCGTTCACACCCCCTCTGCATCTTGCAGA‑GGGCAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACTTCGGTAACTTCGGTaaa < 4:83825/143‑1 (MQ=255) tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa < 2:398465/143‑1 (MQ=255) tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa > 1:398465/1‑143 (MQ=255) tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa > 3:83825/1‑143 (MQ=255) gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc < 4:397169/94‑1 (MQ=255) gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc > 3:397169/1‑94 (MQ=255) cTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGGTc > 3:409423/1‑148 (MQ=255) tGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc > 3:338005/1‑87 (MQ=255) tGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc < 4:338005/87‑1 (MQ=255) ggTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGa > 4:771479/1‑106 (MQ=255) ggTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGa < 3:771479/106‑1 (MQ=255) ggTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCaa < 2:331842/148‑1 (MQ=255) ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt > 4:210024/3‑64 (MQ=25) ttccccTCTGCATGATGCAGAGGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt < 3:210024/62‑1 (MQ=25) cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt > 1:431044/1‑74 (MQ=255) cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt < 2:840882/74‑1 (MQ=255) cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt < 2:431044/74‑1 (MQ=255) cccccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttctt > 1:840882/1‑74 (MQ=255) ccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt < 1:415596/97‑1 (MQ=255) ccTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt > 2:415596/1‑97 (MQ=255) ctctGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct > 3:788513/1‑69 (MQ=255) ctctGCATGATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct < 4:788513/69‑1 (MQ=255) ctctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTTCCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc < 4:409423/148‑1 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg > 1:200810/1‑114 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCgg < 2:200810/114‑1 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa > 1:62646/1‑118 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTaaa < 2:62646/118‑1 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt < 2:854865/139‑1 (MQ=255) tctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCttt > 1:854865/1‑139 (MQ=255) ctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt < 1:1005255/127‑1 (MQ=255) ctGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt > 2:1005255/1‑127 (MQ=255) gCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCg < 2:732214/41‑1 (MQ=255) gCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCg > 1:732214/1‑41 (MQ=255) tcatGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg < 4:145649/83‑1 (MQ=255) tcatGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTg > 3:145649/1‑83 (MQ=255) | TGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGA‑GGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCC > NC_000913/1984969‑1985251 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |