Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 473,861 | 0 | G | T | 73.3% | 36.9 / 57.6 | 60 | G71G (GGC→GGA) | pdeB | cyclic‑guanylate‑specific phosphodiesterase PdeB |
Reads supporting (aligned to +/- strand): ref base G (12/4); new base T (30/14); total (42/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-02 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGT > NZ_CP009273/473774‑473947 | gttATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTcc > 2:821587/1‑90 (MQ=255) gttATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTcc > 1:361691/1‑90 (MQ=255) tATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCtt > 2:748187/1‑90 (MQ=255) tATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCtt > 1:455648/1‑90 (MQ=255) aTCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTg < 1:748187/90‑1 (MQ=255) gATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCa > 2:919225/1‑90 (MQ=255) gATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCa > 2:197259/1‑90 (MQ=255) gATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCa > 2:59613/1‑90 (MQ=255) gATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCa > 1:364393/1‑90 (MQ=255) taAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCtt > 2:696885/1‑90 (MQ=255) ccACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCa < 2:516535/90‑1 (MQ=255) aCTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGc > 2:51089/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCt > 2:440012/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCt > 2:325301/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCt > 2:699465/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCt > 1:932476/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCt > 2:145218/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCt > 1:72499/1‑90 (MQ=255) cTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCt > 1:630460/1‑90 (MQ=255) cTTCATGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGGTCCCTGCCATCTTTCCAGCt > 2:85636/1‑90 (MQ=255) atatAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAgc < 1:999142/90‑1 (MQ=255) ataGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAgcgc > 2:678203/1‑90 (MQ=255) ataGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAgcgc > 2:441811/1‑90 (MQ=255) ataGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAgcgc > 1:697433/1‑90 (MQ=255) aTAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCttt > 2:518635/1‑90 (MQ=255) aaCTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTcc > 1:142954/1‑90 (MQ=255) aaCTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTcc > 1:351795/1‑90 (MQ=255) aaCTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTcc > 1:558597/1‑90 (MQ=255) tGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCtt < 2:85947/90‑1 (MQ=255) tGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCtt < 2:76505/90‑1 (MQ=255) aaGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTtgtg > 2:875451/1‑90 (MQ=255) tACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGc < 1:696885/90‑1 (MQ=255) tACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGc < 2:136812/90‑1 (MQ=255) aCCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGcc < 1:105689/90‑1 (MQ=255) cccGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGATCCTGCAGCGCATCTTTCCCTTGTGTCGCCa > 2:887840/1‑90 (MQ=255) cccGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGATCCTGCAGCGCATCTTTCCCTTGTGTCGCCa > 2:459245/1‑90 (MQ=255) cccGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCa > 1:533111/1‑90 (MQ=255) cccGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCa > 2:847299/1‑90 (MQ=255) ccGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCAc < 1:192170/90‑1 (MQ=255) gCATTTCCATGAGATGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTACTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTAt > 1:871837/1‑90 (MQ=255) cATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTAtt > 2:302132/1‑90 (MQ=255) aTGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGc < 2:1011463/74‑1 (MQ=255) aTGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGc > 1:1011463/1‑74 (MQ=255) agaTGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATa > 1:793995/1‑90 (MQ=255) agaTGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATa > 1:809344/1‑90 (MQ=255) gaTGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAg < 2:48854/90‑1 (MQ=255) gaTGGGCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGGGTCGCCACCTTATTGGCTCTAATaa > 1:979316/1‑89 (MQ=255) gggCTTCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCCAATAGCGa > 2:503903/1‑90 (MQ=255) tCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCATCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGAcgc > 1:392031/1‑88 (MQ=255) tCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCATCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGAcgc < 2:392031/88‑1 (MQ=255) tCGCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCgg < 1:197259/90‑1 (MQ=255) tCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCgg > 2:478000/1‑90 (MQ=255) gCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTtg < 1:431780/54‑1 (MQ=255) gCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTtg > 2:431780/1‑54 (MQ=255) gCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGgga < 1:1059478/90‑1 (MQ=255) gCTACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGgga < 1:821587/90‑1 (MQ=255) tACAGGCAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGggagg < 2:735674/90‑1 (MQ=255) cAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGaa > 1:546239/1‑90 (MQ=255) cAGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAgg < 2:793995/90‑1 (MQ=255) aGCTCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCACTTGTGTCGCCACCTTATTGGCTAGAATAGCGACGCGGGGGGAATAGGt > 1:848455/1‑90 (MQ=255) | GTTATCGATATAAGCCACTTCCTGAATATAGCGATAACTGTAAGATACCCGACGCATTTCCATGAGATGGGCTTCGCTACAGGCAGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGT > NZ_CP009273/473774‑473947 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |