Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,012,918 | A→G | D351G (GAC→GGC) | fliK → | flagellar hook‑length control protein FliK |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,012,918 | 0 | A | G | 98.1% | 165.5 / ‑6.1 | 52 | D351G (GAC→GGC) | fliK | flagellar hook‑length control protein FliK |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (20/31); total (20/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGA > NZ_CP009273/2012838‑2013001 | cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATAc > 1:961545/1‑90 (MQ=255) cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATAc < 2:704293/90‑1 (MQ=255) tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGt > 1:1059056/1‑90 (MQ=255) tCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGt > 1:209323/1‑90 (MQ=255) cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCgg < 1:1353643/88‑1 (MQ=255) cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCgg > 2:1353643/1‑88 (MQ=255) cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGtt > 2:162641/1‑90 (MQ=255) cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGtt > 2:411905/1‑90 (MQ=255) cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc < 1:429000/90‑1 (MQ=255) cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc < 2:1241386/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAACCCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:1340089/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:1488642/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:1399540/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:1211725/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 2:386357/90‑1 (MQ=255) ggCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtctc < 1:387162/90‑1 (MQ=255) gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTa > 2:374910/1‑90 (MQ=255) cccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAg < 2:767954/90‑1 (MQ=255) cccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAg < 2:143701/90‑1 (MQ=255) ccAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAgg < 1:1520068/90‑1 (MQ=255) cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAggg < 1:888017/90‑1 (MQ=255) acagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCgtgt > 2:1170180/1‑90 (MQ=255) gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc > 1:990556/1‑90 (MQ=255) gcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAAc < 2:990556/90‑1 (MQ=255) cGCACCGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGCTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc < 1:1350494/90‑1 (MQ=255) cGCACAGCAATCCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc < 1:1033647/90‑1 (MQ=255) cGCACAGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAgc < 1:297890/90‑1 (MQ=255) caGCAAACCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGCGTAACAGGCAACAgcggcg > 2:801223/1‑90 (MQ=255) aaCCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACaa > 2:1450343/1‑60 (MQ=255) aaCCATGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACaa < 1:1450343/60‑1 (MQ=255) tGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCg > 2:1489702/1‑60 (MQ=255) tGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCg < 1:1489702/60‑1 (MQ=255) tGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt > 2:442238/1‑90 (MQ=255) tGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt > 1:422520/1‑90 (MQ=255) tGAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt > 1:1242185/1‑90 (MQ=255) gAACCTCTGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc < 1:993527/90‑1 (MQ=255) ctctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCt < 1:1122448/90‑1 (MQ=255) tctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc > 1:1207621/1‑40 (MQ=255) tctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGtc < 2:1207621/40‑1 (MQ=38) tctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTa < 2:248911/90‑1 (MQ=255) ctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa > 1:937875/1‑90 (MQ=255) ctGGCGGGGGAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTaa > 2:975522/1‑90 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 2:1392591/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 1:162641/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 2:231633/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 1:411905/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 1:1181211/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 2:920619/90‑1 (MQ=255) gggggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 1:1018756/90‑1 (MQ=255) gggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCaga > 1:470951/1‑90 (MQ=255) gggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa > 2:1332525/1‑88 (MQ=255) gggAAGGCGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCa < 1:1332525/88‑1 (MQ=255) | CTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTCCCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGA > NZ_CP009273/2012838‑2013001 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |