Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,977,855 | (G)6→7 | coding (331/1080 nt) | rffT → | 4‑acetamido‑4,6‑dideoxy‑D‑galactose transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,977,849 | 1 | . | G | 90.9% | 66.5 / 0.8 | 22 | coding (325/1080 nt) | rffT | 4‑acetamido‑4,6‑dideoxy‑D‑galactose transferase |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (9/11); total (11/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 |
CGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑GGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTATTTGCCACCCGCGGCGATTTGAGCTTTTTTGCCAAAACGCA > NC_000913/3977712‑3977982 | cGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑gg < 1:765850/140‑1 (MQ=255) cGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑gg < 2:120138/140‑1 (MQ=255) cTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑gggg < 2:21883/140‑1 (MQ=255) gCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑ggggggg < 1:429269/140‑2 (MQ=255) cGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑gggggggc < 1:171515/140‑3 (MQ=255) aaGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑gggggggcaga < 2:846267/140‑6 (MQ=255) ggTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTg > 1:755470/1‑140 (MQ=255) cATCGCGAAAGCAAAAGCTAACCGTCAGCTGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTAt < 1:285809/140‑1 (MQ=255) gCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCa < 2:493547/140‑1 (MQ=255) cTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCtt < 1:762679/140‑1 (MQ=255) cttcCACGGTCGGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACATGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCa < 2:617152/140‑1 (MQ=255) cttcCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCa < 2:879216/140‑1 (MQ=255) cACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑GGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACg > 2:373083/1‑140 (MQ=255) aGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCt > 2:541991/1‑140 (MQ=255) aTCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCa > 1:4063/1‑140 (MQ=255) tCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCaa > 1:832442/1‑140 (MQ=255) cccACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCaaa < 2:827938/140‑1 (MQ=255) cacTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑GGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGa > 2:797716/1‑140 (MQ=255) acTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAG‑CCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAg < 1:522380/140‑1 (MQ=255) ctggctCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCt < 1:541991/140‑1 (MQ=255) ctggctCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCt < 2:755470/140‑1 (MQ=255) ggctCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCtgt > 1:641517/1‑140 (MQ=255) gctCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCtgtg < 1:364196/140‑1 (MQ=255) tctGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTa > 2:525753/1‑140 (MQ=255) tgctgAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTAtt > 1:519226/1‑140 (MQ=255) ctgAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTATTTg > 1:97117/1‑140 (MQ=255) ccagccagTTTTTCTGGCATATCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTATTTGCCACCCGCGGCGATTTGAg < 2:4063/140‑1 (MQ=255) atatCTGGGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTATTTGCCACCCGCGGCGATTTGAGCTTTTTTGCCAAAACGCa > 1:246965/1‑140 (MQ=255) | CGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTGGCATATCT‑GGGGGGCAGACCTGTATGAGCTTTCCAGTGGCTTGAGATATAAGCTTTTTTACCCACTACGTCGCCTGGCGCAAAAGCGAGTCGGCTGTGTATTTGCCACCCGCGGCGATTTGAGCTTTTTTGCCAAAACGCA > NC_000913/3977712‑3977982 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |