Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,123,164 | C→T | 11.1% | D453N (GAT→AAT) | cpsB ← | mannose‑1‑phosphate guanyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,123,164 | 0 | C | T | 11.1% | 98.1 / 3.6 | 36 | D453N (GAT→AAT) | cpsB | mannose‑1‑phosphate guanyltransferase |
Reads supporting (aligned to +/- strand): ref base C (14/18); new base T (2/2); total (16/20) | |||||||||||
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.12e-01 |
AATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACC > NC_000913/2123026‑2123291 | aaTTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAAtt > 6:149417/1‑149 (MQ=255) ttttattttTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAAtttt > 6:214043/1‑149 (MQ=255) tGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGtt > 2:211107/1‑149 (MQ=255) tCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGATTTTTCCCCGGGttt < 6:8378/148‑1 (MQ=255) ttATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCa < 4:90579/149‑1 (MQ=255) acacCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCaa < 5:149417/148‑1 (MQ=255) acacCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt > 3:146985/1‑149 (MQ=255) acCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt < 8:227538/147‑1 (MQ=255) acCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt > 7:227538/1‑147 (MQ=255) cgaattcattGCATGATGCGACGTTTACACCCTGCCTTAGCGATACGCGAAACTCACCACATAATAATCTTCGAGATAAGTTCCGGAGCGCACTTCAATTAAATCTATCTGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGcccca < 6:123326/145‑2 (MQ=255) gCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAg > 6:134488/1‑149 (MQ=255) attcattGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCg < 2:69080/149‑1 (MQ=255) ttGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAata < 4:76352/149‑1 (MQ=255) aTGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGATCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAAt < 3:16974/147‑1 (MQ=255) cGACGTTTACTCCCGTCCTTCGCGATCCGCGACACGCACCACATAATCCTCTTCGAGATAAGAGACGGAGCGCAGTTCAATTAAATCGTGCTTAATTTTCCACGGGTTTTCAATGCAATGCTTCGCCCGAAGCTGTTTATAAATGGAc < 4:146985/148‑1 (MQ=255) gACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTc > 6:211612/1‑149 (MQ=255) gTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGtt < 1:211107/149‑1 (MQ=255) cacCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCAc < 5:211612/149‑1 (MQ=255) ccgtccgtAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCAcc > 7:243802/1‑147 (MQ=255) gtAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAATAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGtt < 3:49519/149‑1 (MQ=255) cGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGat < 5:134488/149‑1 (MQ=255) cgcgAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCAc > 3:188553/1‑149 (MQ=255) gcgAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCAc > 2:229734/1‑148 (MQ=255) gCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGaa < 2:140740/102‑1 (MQ=255) gCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGaa > 1:140740/1‑102 (MQ=255) catcatCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt < 1:17797/78‑1 (MQ=255) catcatCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt > 2:17797/1‑78 (MQ=255) ctctTCGATATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGAAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGGTTGATATCACCATCAATGGTGACTTTTGCCg < 5:214043/149‑1 (MQ=255) ctTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGtt > 3:163036/1‑149 (MQ=255) agaTATGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCTGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg < 4:29193/129‑1 (MQ=255) agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg < 4:32346/129‑1 (MQ=255) agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg > 3:32346/1‑129 (MQ=255) agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCACCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg > 3:29193/1‑129 (MQ=255) cGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCTGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAaccacc < 8:243802/149‑1 (MQ=255) acatCAATCACATCGCGCGGACCTTTCCCCGGGGATTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGCTTTCACCAAGCAGGTTGATATCACCATCAATGGTGACTTTAGCCGTGCCCGCGa < 8:181689/130‑1 (MQ=255) aCTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGa > 7:181689/1‑133 (MQ=255) | AATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACC > NC_000913/2123026‑2123291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |