Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,852,466 | 0 | T | G | 46.5% | 34.6 / 31.6 | 72 | V672G (GTG→GGG) | mutS | DNA mismatch repair protein MutS |
Reads supporting (aligned to +/- strand): ref base T (23/15); new base G (0/33); total (23/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.65e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852378‑2852555 | caAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTAtggtg < 2:2885975/90‑1 (MQ=255) caAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTAtggtg < 2:2808923/90‑1 (MQ=255) aaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGgaga < 1:1025875/90‑1 (MQ=255) aaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGgaga < 1:616169/90‑1 (MQ=255) tctcgatcgCTTTTTTCCCCGGGTGGGCGCGGAAGATGCCCTGGCGCCCGGCCC‑CTCACCTTTTAGGGGGGAGATGACTGAAACCGCCaa < 2:691249/88‑1 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 1:2250201/1‑90 (MQ=255) atcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAat > 1:1983392/1‑90 (MQ=255) gatcgCTTCTTTCCCCGGGTGGGCGGGGAAGATGCCCGGGCGTCGGGGCG‑CTCACCTTTTATGGTGGAGATGACTAAAACCGCCAATAtt < 2:874513/90‑1 (MQ=255) atcgCATTTTTCCCCTCGTAGGCGGGGCAGACGACCTGGCGTCCGGGCC‑CTCACCTTTTATGGGGGAGATGACTAAAACCGCCAATAttt < 1:698125/90‑1 (MQ=255) atcgCATCTTTACCCGGGTGGGCGCGGCAGATGCCCTGGCGTCCGGGCC‑CTAAACTTTTAGGGTGGAGATGACTAAAACCGCCAATAttt < 1:1993178/90‑1 (MQ=255) gCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTAc > 1:222158/1‑90 (MQ=255) tCTTTACCCTCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATa > 2:283922/1‑90 (MQ=255) tttCCCCGGGTGGGCGCGGCAGATGCCGTGGCCTCCGGGCC‑CACAACCTTTATGGTGGAGATGACTAAACCCGCCAATATTTTACATAAc < 2:236596/90‑1 (MQ=255) tttCCCCGCGTGGCCGCGGCAGATGCCCTGGCGCCCGGGCG‑CCCACCCTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAAc < 1:2286576/90‑1 (MQ=255) tccccGGGTGGGGGGGGCAGATGCCCGGGCGCCCGGGCC‑CTCAACCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 2:2275551/88‑1 (MQ=255) tccccGGGTGGGCGCGGCAGATGACCTGGCGCCCGGGCG‑CTACACTTTTATGGTGGAGATGACTGAACCCGCCAATATTTTACATAACGc < 2:1316428/88‑1 (MQ=255) tccccGCGTGGGCGCGGCAGATGCCCTGGCGCCGGGGCG‑CTCACCTTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 1:2417612/88‑1 (MQ=255) tccccGCGTAGGCGCGGCAGATGCCCCGGCGCCGGGGCC‑CCCACCCTTTATGGTGGAGATGACAAAAACCGCCAATATTTTACATAACGc < 2:2109658/88‑1 (MQ=255) tccccGCGGGGGCGCGGCAGATGCCCGGGCCCCCGGGCC‑CTAACCTTTTGGGGGGGAGATGACTAAACCCGCCAATATTTTACATAACga < 2:3254076/88‑2 (MQ=255) tccccGCGGAGGCGCGGCAGATGCCCTGGCGTCGGGGCG‑CTCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 1:2711631/88‑1 (MQ=255) tccccGCGGAGGCGCGGCAGATGCCCGGGCGCCCGGGCC‑CCCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGc < 1:2418057/88‑1 (MQ=255) tACCCGCGTAGGCGCGGCAGATGCCCTGGCGCCCGGGCC‑CTCACCTTTTAGGGTGGAGATGATAAAAACCGCCAATATTTTACATAACGc < 1:1710396/90‑1 (MQ=255) tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACga > 1:3254076/1‑89 (MQ=255) tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc > 2:1128192/1‑90 (MQ=255) cccGCGTGGGCGGGGCAGATGCCCGGGCGCCCGGCCC‑CCCACCTTTTTGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCa < 2:1079791/90‑1 (MQ=255) ccGCGTGGGCGCGGCAGATGCCCGGGCGCCCGGGCC‑CCCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc < 1:2785571/90‑1 (MQ=255) ccGCGTCGGGGCGGCAGAGGCCCGGGCGTCGGGGCC‑CTAACCTTTTATGGGGGAAATGACTAAAACCGCCAATATTTTACATAACGCCAc < 1:654613/90‑1 (MQ=255) ccGCGGGGGCGCGGAAAAGGCCCGGGCGCCGGGGCC‑CTCACCTTTTAGGGGGGAGATGCCTGAAACCGCCAATATTTTACATAACGCCAc < 1:2356580/90‑1 (MQ=255) cgcgTAGGCGCGGCAGATGCCCTGGCGCCGGGGCG‑CTCCCCCTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAcc < 1:316622/90‑1 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 2:2030064/1‑90 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 1:3139304/1‑90 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 1:866091/1‑90 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 1:397068/1‑90 (MQ=255) cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGa > 1:812164/1‑90 (MQ=255) gtgggCGCGGCAGATGCCCTGGCGCCGGGGCG‑CTCAACCTTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGaa < 1:2367513/87‑1 (MQ=255) gTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGaa > 1:332012/1‑90 (MQ=255) gggggcgGAAGATGCCGGGGCGGCCGGGCC‑CCCAACCTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATa < 2:1301856/86‑1 (MQ=255) gggCGGGCCAGTTGACC‑GGGGGCCGGGCGCCTAAACTTTTTGGGGGGAGATGACGAAAACCGCCAATATTTTACATAACGCCACCGAATa < 1:657822/89‑1 (MQ=255) gggCGCGGCAGATCCCCGGGCGCCCGGGCG‑CCCACCTTTTAGGGGGGAGATGACTGAACCCGCCAATATTTTACATAACGCCACCGAATa < 2:1248375/89‑1 (MQ=255) gggCGCGGAAGATGCCCTGGCGTCGGGGCG‑CTCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATa < 1:150394/89‑1 (MQ=255) aGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATa > 1:515884/1‑90 (MQ=255) ggCGCGGCAGATGCCCTGGCGCCCGGGCC‑CTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:2967790/90‑1 (MQ=255) ggCGCGGCAGATGCCCTGGCGCCCGGGCC‑CCCACCCTTTATGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:701177/90‑1 (MQ=255) ggCGCGGCAGATGCCCGGGCGCCGGGGCG‑CTCAACCTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:1427286/90‑1 (MQ=255) ggCGCGGCAGATGCCCGGGCGCCGGGGCG‑CCCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:2990367/90‑1 (MQ=255) ggCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCACCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:3126296/90‑1 (MQ=255) ggCGCGGCAGATGACCTGGCGCCGGGGCG‑CTCACCTTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:955247/90‑1 (MQ=255) gcgcGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACa > 2:433542/1‑90 (MQ=255) gcgGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGt > 1:235285/1‑90 (MQ=255) cgGCAGATGCCCTGGCGTCGGGCCC‑CTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 1:61553/90‑1 (MQ=255) ggCAGATGCCCGGGCGTCCGGGCC‑CTCCCCCTTTAGGGGGGGGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:225393/90‑1 (MQ=255) ggCAGATGCCCGGGCGCCCGGGCC‑CCCACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 1:3284448/90‑1 (MQ=255) ggCAGATGCCCGGGCGCCCGGGCC‑CCCAACTTTTAGGGGGGAGATGATTAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:1229569/90‑1 (MQ=255) ggCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt > 2:923745/1‑90 (MQ=255) ggCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt > 2:1580675/1‑90 (MQ=255) ggCAGAGGACCTGGCGTCCGGGCC‑CACATCCTTTAGGGTGGAGATGATTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:1049971/90‑1 (MQ=255) gaaGATGCCCGGGCGCCGGGGCC‑CTAACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTg < 1:1548820/88‑1 (MQ=255) cAGATGCCCGGGCGCCCGGGCC‑CCCACCCTTTAGGGGGGAGATGACGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTgg < 1:43391/90‑1 (MQ=255) cccTGGCGCCGGGCCG‑CCCAACCTTTATGGGGGAGATGACTTAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa < 1:2865431/89‑1 (MQ=255) ccTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAAt > 2:1349611/1‑90 (MQ=255) cTGGCGCCGGGGCG‑CTCACCTTTTAGGGGGGAGATGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATg < 1:3065737/90‑1 (MQ=255) tGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATgg > 1:959736/1‑90 (MQ=255) tGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATgg > 2:1353679/1‑90 (MQ=255) ggCGTCCGGGCG‑CTCACCCTTTGGGGGGGAGATGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGa < 2:397068/90‑1 (MQ=255) ggCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGa > 2:492802/1‑90 (MQ=255) ggCGTCCGGGCA‑CCAAACTTTTAGGGGGGAGTTGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGa < 2:1868289/90‑1 (MQ=255) cGTCCGGGCG‑CTAACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATg < 2:3212642/90‑1 (MQ=255) g‑cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGc > 2:1764566/1‑90 (MQ=255) g‑cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGc > 2:1573502/1‑90 (MQ=255) cAACTTTTTGGGTGGAGATGCCGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg < 1:2290314/90‑1 (MQ=255) acttttAGGGGGGAGATGACTAAAACCGCCAATA‑TTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa < 1:170058/87‑1 (MQ=255) ccTTTAGGGTGGAGATGACAAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa < 2:2798853/90‑1 (MQ=255) ttATGGCGGACATCACTCAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt < 2:1869690/90‑1 (MQ=255) ttAGGGTGGAGATGATTAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt < 1:3286298/90‑1 (MQ=255) gtggAGATGACAGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAcc < 2:2567015/90‑1 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:2405369/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:1451874/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:2323191/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:2671772/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:2386452/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:2336858/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:1997483/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:1576875/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:1491884/1‑90 (MQ=255) tggAGATCACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:360305/1‑90 (MQ=255) | CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCG‑CTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852378‑2852555 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |