Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,630,918 | G→A | 13.7% | D114N (GAC→AAC) | lon‑II → | DNA‑binding ATP‑dependent protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,630,918 | 0 | G | A | 13.7% | 65.7 / 2.8 | 29 | D114N (GAC→AAC) | lon‑II | DNA‑binding ATP‑dependent protease |
Reads supporting (aligned to +/- strand): ref base G (13/12); new base A (2/2); total (15/14) | |||||||||||
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.89e-01 |
CGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCT > NC_002947/2630779‑2631060 | cGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGa < 2:605308/150‑1 (MQ=255) cAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa > 1:477472/1‑150 (MQ=255) acctgctgAAGCTGCCTGCTGACACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTGCGAAGTTGAAGGGCACATCCGTGCCGATGTCTCCCTGATCGACCAAACCGATACTGCCGAGCGCGAGTCGGCAg > 2:345680/3‑150 (MQ=255) gctgAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGtt < 2:477472/150‑1 (MQ=255) tgAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCg < 2:263649/150‑1 (MQ=255) tGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGcgcg < 2:608304/129‑1 (MQ=255) tGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGcgcg > 1:608304/1‑129 (MQ=255) cttgacggTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa < 2:415474/118‑1 (MQ=255) cttgacggTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGaa > 1:415474/7‑124 (MQ=255) gTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAg < 2:138793/150‑1 (MQ=255) cAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGtt < 2:352373/150‑1 (MQ=255) tGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCAACGAAACCGATGCTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTGGAACAGt > 1:111104/1‑150 (MQ=255) cGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCt < 2:472073/150‑1 (MQ=255) agcagcGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTgg < 1:193441/150‑1 (MQ=255) gcctgtcCGTAGTTGAAGGGCAGGTCCTTGCCGAAGTCTCCCTGATCAACGAAACCGATACTGCCGGGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAg < 2:111104/145‑1 (MQ=255) ttCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTg > 1:236842/1‑150 (MQ=255) tctcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgc < 1:56468/150‑1 (MQ=255) tcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgcgc > 1:314306/1‑150 (MQ=255) tcCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGgcgc > 2:482015/1‑150 (MQ=255) tcCCTGATCAACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGaa < 1:598093/75‑1 (MQ=255) tcCCTGATCAACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGaa > 2:598093/1‑75 (MQ=255) ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTTCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCAAGCCGAAGTCCTGTCTTTGCTGAACAGCATCGACGAGCGAGGGCGCCt > 1:299317/1‑150 (MQ=255) ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTTAACAGCATCGACGAGCCAGGGCGCCt > 2:494644/1‑150 (MQ=255) ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCt > 1:610636/1‑150 (MQ=255) ccTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCCTCGCTGAACAGCATCGACGAGCCAGGGCGCCt > 1:412249/1‑150 (MQ=255) cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCa > 2:206599/1‑102 (MQ=255) cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCa < 1:206599/102‑1 (MQ=255) cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCt < 2:243345/122‑1 (MQ=255) cTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCt > 1:243345/1‑122 (MQ=255) | CGGTACCGTTGCCACCGTGCTGCAACTGCTGAAGCTGCCTGATGGCACCGTCAAGGTGCTGGTCGAGGGCGAGCAGCGGGGTGCGGTCGAGCGCTTCTCCGAAGTTGAAGGGCACATCCGTGCCGAAGTCTCCCTGATCGACGAAACCGATACTGCCGAGCGCGAGTCGGAAGTGTTCGTGCGCACGCTGCTGTCGCAGTTCGAACAGTACGTTCAGCTGGGCAAGAAGGTGCCAGCCGAAGTGCTGTCTTCGCTGAACAGCATCGACGAGCCAGGGCGCCT > NC_002947/2630779‑2631060 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |