Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,691,100 | 0 | T | G | 40.5% | 50.1 / 11.2 | 37 | V428G (GTG→GGG) | bcsE | cellulose biosynthesis c‑di‑GMP‑binding protein BcsE |
Reads supporting (aligned to +/- strand): ref base T (4/18); new base G (14/1); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.71e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAAC > NZ_CP009273/3691011‑3691188 | ggTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGt > 2:824437/1‑90 (MQ=255) gggTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGc < 1:828185/86‑1 (MQ=255) gggTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGtt < 1:276753/90‑1 (MQ=255) gggTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGc > 2:828185/1‑86 (MQ=255) gggTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGtt > 1:474143/1‑90 (MQ=255) tCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTtctct > 1:160544/1‑90 (MQ=255) tCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCATACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTtttct > 2:359324/1‑90 (MQ=255) aaCAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTc < 2:767691/90‑1 (MQ=255) ccTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTaa < 2:474143/90‑1 (MQ=255) tGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTTTTCTCATTTTTTCGGGTTAACg > 2:266340/1‑90 (MQ=255) cTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAAACGGCTGGGGGTGGTTTTCTCCTTTTTTCGGATTAACGGTCTg > 2:676432/1‑90 (MQ=255) gtgtCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGTTGGGGCTGTTTTTCTCATTCTGTCGGGTTAACGATCTggg > 1:926051/1‑89 (MQ=255) tCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAc < 1:824437/90‑1 (MQ=255) tCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATAc < 2:640553/90‑1 (MQ=255) cTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGt < 2:526047/90‑1 (MQ=255) cTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGt < 2:781571/90‑1 (MQ=255) cTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGt < 1:90289/90‑1 (MQ=255) cTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGt < 1:1039636/90‑1 (MQ=255) cTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGt < 1:676432/90‑1 (MQ=255) aaCCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTTTTCTCCTTTTGTCGGGTTAACGATCTGGAAACCGCGGTg > 2:90289/1‑90 (MQ=255) tACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCa > 1:24740/1‑90 (MQ=255) aCCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTtctc < 1:1012383/48‑1 (MQ=255) aCCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCat > 1:348653/1‑90 (MQ=255) aCCGGCGATATCATGACCATTGGCGGTAATCGGCTGGGGCTGTTtctc > 2:1012383/1‑48 (MQ=255) ggCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTTTGTCGGGTTAACGATCTGGGTACCGCGTTGAATCATAtt > 1:847142/1‑90 (MQ=255) tCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCAt < 1:266340/90‑1 (MQ=255) tCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCAt < 2:999860/90‑1 (MQ=255) aTGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTg < 1:445440/90‑1 (MQ=255) aTGACCATTGGCGGTAATCGGCTGGGGCTGTTTTTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGGATCATATTTTTCCAttt > 2:199155/1‑89 (MQ=255) gACCATTGGCGGTAATCGGCTGGGGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGcc > 1:1035880/1‑90 (MQ=255) aTTGGCGGTAATCGGCTGGGGCTGTTTTTCTCATTTTTTCGGGTTTACGATCTGGATACCGCGGTGAATCATATTTTTCCCTTTCCTACt > 2:82597/1‑90 (MQ=255) ttGGCGGTAATCGGCTGGGGCTGTTTTTCTCATTTTTTCGGGTTAACGGTCTGGATACCGCGGTGAATCATTTTTTTCCCTTGCCTACTg > 2:624730/1‑90 (MQ=255) gCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCg < 2:248069/90‑1 (MQ=255) tAATCGGCTGGGGCTGTTTTTCTCCTTTTTTCGGATTAAAGGTCTGGGGACCGCGTTGAATCATTTTTTTCCCTTTCCTACTGGGGACAt > 2:961790/1‑90 (MQ=255) ggCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTtct > 2:660562/1‑90 (MQ=255) tGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCaa < 1:97228/90‑1 (MQ=255) ggTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCaaa < 2:1035880/90‑1 (MQ=255) gTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAAc < 1:242901/90‑1 (MQ=255) | GGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTGGTGCTGTTTCTCTCATTCTGTCGGATTAACGATCTGGATACCGCGTTGAATCATATTTTCCCATTGCCTACTGGCGACATTTTCTCAAAC > NZ_CP009273/3691011‑3691188 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |