Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,126,661 | 0 | T | G | 49.3% | 56.4 / 38.7 | 71 | H219P (CAC→CCC) | yghS | dTMP kinase |
Reads supporting (aligned to +/- strand): ref base T (1/35); new base G (35/0); total (36/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.67e-19 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGT > NZ_CP009273/3126573‑3126751 | cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 1:2993080/90‑1 (MQ=255) cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 1:2655321/90‑1 (MQ=255) cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 1:2529698/90‑1 (MQ=255) ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑c < 2:2761272/90‑1 (MQ=255) ttAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑c < 2:1160786/90‑1 (MQ=255) aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGcc < 1:1182064/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATcc < 1:956047/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATcc < 2:2385670/90‑1 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATcc > 1:2385670/1‑90 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGCGccctcc > 2:1343545/1‑87 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCTcc > 1:2451978/1‑90 (MQ=255) tcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTcc > 2:2641503/1‑90 (MQ=255) ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑TGCCCCTCCa > 2:1611124/1‑90 (MQ=255) ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCATccc > 1:3154195/1‑89 (MQ=255) ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCATCCa > 2:1696389/1‑90 (MQ=255) cccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAg < 2:581847/90‑1 (MQ=255) cTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCTcccggt > 1:3073227/1‑86 (MQ=255) aTCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCa > 1:306214/1‑90 (MQ=255) tCCATTTTTTTGATCAGAAAAGCGCGGGGTGATATGCCCCTATTGACGCTTGCTGAGTTTCATGGGCGGGGTGT‑CGCCCCTCCAGTTCaa < 2:2008243/90‑1 (MQ=255) cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAgg < 2:2043095/90‑1 (MQ=255) cATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAgg < 1:1444461/90‑1 (MQ=255) ttttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCCCCCGTTCaaaggt > 2:2473128/1‑86 (MQ=255) tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCTCCAGATGAAGGagc > 1:1349907/1‑88 (MQ=255) ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCAAGGGTCt > 1:1887680/1‑90 (MQ=255) tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCtt < 1:2578926/90‑1 (MQ=255) tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑AGCCCCTCCCGATCCAGGGTCtt > 1:1457546/1‑90 (MQ=255) gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCTCCCGTTCAAGGGTCtttt > 2:1775499/1‑89 (MQ=255) gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCCCCCGGTCAAGGGACtttt > 2:2261763/1‑89 (MQ=255) gTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCCCCAGTTCAAGGAGCtttt > 2:1003248/1‑89 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcg < 1:532407/90‑1 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcg < 1:1464286/90‑1 (MQ=255) gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcg < 1:1602094/90‑1 (MQ=255) aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcgc < 1:419597/90‑1 (MQ=255) aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcgc < 1:1067026/90‑1 (MQ=255) aTCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTgcgc < 1:1135280/90‑1 (MQ=255) tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGcc < 1:2870180/90‑1 (MQ=255) tCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGcc < 2:1527516/90‑1 (MQ=255) ggAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCAtt < 2:1468228/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGGAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGAATg > 1:1954362/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg < 2:699912/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg < 1:1003248/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATg < 2:502999/90‑1 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCTCCCGGTTAAGGAGCTTTGCGCCATTGAATg > 2:1915276/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCAAGGGTCTTTTGGCCCTTGAATg > 2:1067026/1‑90 (MQ=255) gcgcGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGGTCAAGGATCTTTTCGCCCTTGAATg > 2:2336734/1‑90 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg > 2:2501607/1‑50 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 1:2501607/50‑1 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGAATGTc > 1:2630548/1‑90 (MQ=255) gcggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCAAGGATCTTTTCGCCCTTGGATGTc > 1:1234296/1‑90 (MQ=255) cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg > 2:299009/1‑49 (MQ=255) cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtg < 1:299009/49‑1 (MQ=255) cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCTCCCGGTCCAGGATCTTTTGGCCATTGAATGGCa > 2:2971498/1‑90 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑GGCCCCTCCAGGTCAAGGATGTTTGGGGCATTGAATGTCAGt > 1:1897310/1‑90 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCCCCAGTTTAAGGATCTTTGCGCCCTTGGATGTCAGt > 1:1193046/1‑90 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑TGCCCCTCCAGTTCAAGGATCTTTTGGCCATTGAATTTTAGt > 1:699912/1‑90 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑GGCCCCTCCCGTTCGAGGATCTTTTCGGCATTTAATGtttgt > 1:2131600/1‑86 (MQ=255) gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCAGTTTAAGGATCTTTTTGCCCTTGAATGTCAgg > 1:2790575/1‑89 (MQ=255) tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc < 1:2463155/90‑1 (MQ=255) tGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgc < 1:2641503/90‑1 (MQ=255) gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGG‑CGCCCCCCCCGGTCAAGGAACTTTGCGCCATTGAATGTGAGTTgtggggga > 2:1312344/1‑86 (MQ=255) tAGTGACGCTTGCAGAATTTCATCGGCGGGGGGT‑CGCCCCTCCAGTTCAAGGGTCTTTTCGCCCTTGAATGTCAGTTGCGGCGGAACggg > 1:2710670/1‑89 (MQ=255) aCGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAttt < 1:1696389/90‑1 (MQ=255) cTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGAtttttt < 2:2710670/90‑1 (MQ=255) ttGCAGAATTTCATCGGCGGGGGGT‑CGCCCCCCCCGGTCAAGGGGCTTTTCGCCCTTTTATGTCAGTTGCGGGGGGACGGCGGttttttt > 2:671013/1‑89 (MQ=255) tGCAGAATTTCATCGGCGGGGGGT‑CGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGGAACGGCGAtttttttc > 1:3164133/1‑88 (MQ=255) tGCAGAATTTCATCGGCGGGGGGG‑CGCCCCTCCCGTTCAAGGATTTTTGTGCCCTTTAATGTCAGTTGGGGGGGGACGGCGGtttttttc > 2:1306234/1‑88 (MQ=255) tGCAGAATTTCATCGGCGGGGGGG‑CGCCCATCCAGCTCAAGGAGCTTTGCGGCCTTGAATGGCAGGTTCGGCGGGACTGCGGtttttttc > 2:2883760/1‑88 (MQ=255) gAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTa < 1:206774/90‑1 (MQ=255) aTTTCATCGGCGGGGGGT‑CGCCCCTCCAGTTCAAGGATCTTTTCGCCCTTTAATGTCAGTTGCGGGGGAACGGCGGTTTTTTTCTTTaaa > 1:2402047/1‑89 (MQ=255) aTTTCATCGGCGGGGGGG‑CGCCCCTCCCGGTCAAAGATCTTTTCGCCCTTTAATGTCAGGTGGGGGGTAAAGGGGGTTTTTTTCTGTaaa > 1:2752185/1‑89 (MQ=255) ttCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGc < 1:1657700/90‑1 (MQ=255) ttCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGc < 2:529936/90‑1 (MQ=255) aTCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTga < 1:1487113/90‑2 (MQ=255) cggcggGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCa < 2:2776214/90‑1 (MQ=255) cggcggGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCa < 2:286211/90‑1 (MQ=255) cggGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGtt < 1:2155808/90‑1 (MQ=255) gggTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg < 2:1496586/90‑1 (MQ=255) gggTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtgg < 1:2971498/90‑1 (MQ=255) tgt‑CGCCCATCCAGTTCAAGG‑TCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTggtggt < 1:2832509/90‑1 (MQ=255) | CGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGT‑CGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGT > NZ_CP009273/3126573‑3126751 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |