Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,621,984 | 0 | A | C | 70.0% | 1.4 / 13.4 | 20 | V105G (GTC→GGC) | pdeF | cyclic‑guanylate‑specific phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/14); total (6/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.58e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT > NZ_CP009273/2621895‑2622072 | gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGa > 2:655070/1‑90 (MQ=255) gaaaaCAAAATCTGCCAGAATATACGCGGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTaa < 2:221795/88‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAAcc < 2:244696/90‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGACTAAAGCCCCTGTAAcc < 1:659492/90‑1 (MQ=255) aaaCACAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAAcc < 2:74510/90‑1 (MQ=255) aaaCAAAATCTCCCAAAATATACGCTGGAAAATTAAACGGGTATCCCCATGTGAAACGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAAcc < 1:371451/90‑1 (MQ=255) cacaATCTGCAAAAATATACGCTGGGAAATAAAACGGGTATCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGcc < 1:105581/90‑1 (MQ=255) acaATCTGCAAAAATATCCGGGGGAAAATTAAACGGTTACCCCCATGTAAACGGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGccc < 2:739290/90‑1 (MQ=255) aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCa < 1:256924/90‑1 (MQ=255) gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACaga > 1:49221/1‑90 (MQ=255) cAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAgaga < 1:340748/90‑1 (MQ=255) ataCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgata > 2:258346/1‑90 (MQ=255) cgcgGGAAAATAAACCGGGTATCCCCATGTAAAACTTTTTTCCGCGGGGAAAAAAAGCCCTTATAACCGCCCCAACAGAGaatgataatg < 2:731500/90‑1 (MQ=255) ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATgataatgataat > 1:237515/1‑90 (MQ=255) gTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt > 1:225677/1‑90 (MQ=255) tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc > 1:540211/1‑90 (MQ=255) cgcgGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc < 1:655070/90‑1 (MQ=255) ccGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 2:171985/90‑1 (MQ=255) gggggAATAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc < 1:192567/89‑1 (MQ=255) gggggAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc < 2:38910/89‑1 (MQ=255) gACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATAAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt > 1:370536/1‑90 (MQ=255) | GCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTT > NZ_CP009273/2621895‑2622072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |