Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,362,264 | C→G | 11.5% | Q48H (CAG→CAC) | tpiA ← | triose phosphate isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,362,264 | 0 | C | G | 11.5% | 150.8 / 11.7 | 52 | Q48H (CAG→CAC) | tpiA | triose phosphate isomerase |
Reads supporting (aligned to +/- strand): ref base C (23/23); new base G (2/4); total (25/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01 |
CTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGC > NC_002947/5362117‑5362409 | cTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCtg > 1:821843/1‑150 (MQ=255) tGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCAGGTGATtgatg < 2:713793/150‑1 (MQ=255) ggCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCAGGTGATtgatgg < 2:990855/150‑2 (MQ=255) ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAAc > 1:796492/1‑150 (MQ=255) ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAAc > 2:343637/1‑150 (MQ=255) ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAAc > 2:423408/1‑150 (MQ=255) ccAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAAc > 1:1059861/1‑150 (MQ=255) cAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCAGGTGATtgatggc < 1:1075602/148‑3 (MQ=255) cAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCAGGTGATtgatggc > 2:1075602/1‑146 (MQ=255) aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCAGGTGATtgatggcctg < 2:370130/150‑6 (MQ=255) aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACaa > 2:676116/1‑150 (MQ=255) aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACaa > 1:261612/1‑150 (MQ=255) tCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAgg > 1:115536/1‑150 (MQ=255) cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGCGGTCAGGGCGGCGCCAGTGATTTGGTTATGTGCGCGGAGCAGGAA‑CACCTGATTGATGAACAAGGCCGGTgg > 1:569784/1‑150 (MQ=255) cAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTgg < 1:343637/150‑1 (MQ=255) ttGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAAt > 1:607432/1‑150 (MQ=255) aCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc < 2:138792/150‑1 (MQ=255) cTGNTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGa > 2:703012/1‑150 (MQ=255) aCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAg > 2:532984/1‑123 (MQ=255) aCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAg < 1:532984/123‑1 (MQ=255) gctgctCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGc > 2:1002458/1‑150 (MQ=255) tctggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGc > 2:1084750/4‑150 (MQ=255) gCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCg > 1:949999/1‑150 (MQ=255) aCTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCgaga < 1:703012/150‑1 (MQ=255) aCTGCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGGGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCgaga < 1:140124/150‑1 (MQ=255) gcgcccTGTGCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATTAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg < 2:359303/147‑1 (MQ=255) gcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg < 1:423408/150‑1 (MQ=255) cTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACATTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATCCGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcag < 1:868820/150‑1 (MQ=255) gggTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagctcagc < 2:821843/150‑1 (MQ=255) agctgAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCa < 1:1030865/95‑1 (MQ=255) agcagAATTATGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCa > 2:1030865/1‑95 (MQ=255) gcagAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAgcg > 2:807178/1‑150 (MQ=255) gAATTCTGTCCGCCGACAGTGATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg < 2:417621/150‑1 (MQ=255) aaTTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGACTGATGAACAAGGCTGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCggg < 2:607432/150‑1 (MQ=255) ttCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg > 2:448573/1‑147 (MQ=255) ttCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCgg < 1:448573/147‑1 (MQ=255) cgcCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCa < 2:115536/150‑1 (MQ=255) aGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa > 1:594248/1‑150 (MQ=255) aGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa < 1:1002252/150‑1 (MQ=255) aGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa < 1:925247/150‑1 (MQ=255) aGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa < 1:676116/150‑1 (MQ=255) gTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa > 1:463835/1‑86 (MQ=255) gTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa < 2:463835/86‑1 (MQ=255) tAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTcc > 2:39603/1‑147 (MQ=255) tAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTcc < 1:39603/147‑1 (MQ=255) tttCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTAc < 1:1002458/150‑1 (MQ=255) cTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTCCCAg < 2:885681/150‑1 (MQ=255) tttACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg < 1:1035261/100‑1 (MQ=255) tttACCGGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg > 2:1035261/1‑100 (MQ=255) gttcATCAATCAGGTGATTGATGGACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGc > 1:659919/4‑150 (MQ=255) tcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGc > 2:771935/1‑150 (MQ=255) tcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCGACCATAGGGCGACGCATGc < 2:949999/150‑1 (MQ=255) | CTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGC‑CAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGC > NC_002947/5362117‑5362409 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |