Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,419,801 | 0 | T | G | 66.7% | 5.2 / 15.9 | 18 | L41V (TTG→GTG) | cycA | D‑serine/D‑alanine/glycine transporter |
Reads supporting (aligned to +/- strand): ref base T (1/5); new base G (12/0); total (13/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTT > NZ_CP009273/4419712‑4419885 | atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgt < 1:287266/90‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt < 2:14564/90‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt < 1:254221/90‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt < 1:544837/90‑1 (MQ=255) tCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGgtt < 1:530495/90‑1 (MQ=255) cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTgg > 1:380996/1‑90 (MQ=255) cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCTTTGTACGGGGGTGTTTTTgg > 2:14979/1‑90 (MQ=255) cGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCATTGGTACGGGGGTGGTTATgg > 1:534816/1‑90 (MQ=255) gAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCt > 1:148843/1‑90 (MQ=255) gAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGGACGGGGGTGTTTTTTGGGGCt > 1:509340/1‑90 (MQ=255) aGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCCTTTGTACGGGGGTGTTTTTTGGGGCTggga > 2:546095/1‑88 (MQ=255) gTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGGCTGGCaa > 1:374855/1‑90 (MQ=255) tCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaaa < 1:546095/90‑1 (MQ=255) tCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCaaa < 2:495494/90‑1 (MQ=255) cGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGGGCCCTTGGTACGGGGGTGTTTTTTGGGTCTGGCaaaa > 2:540123/1‑90 (MQ=255) gcAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTTGGGGCTGGCAAAACGATTAGcc > 2:47290/1‑90 (MQ=255) ccGACATATTCAGCTTATTGCCATTGGCGGGGGCATTGGTACGGGGGGGTTTTTGGGGTCTTGCAAAAAGATTTGCCCTGCCGGGGCGTc > 2:301623/1‑90 (MQ=255) cGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCg < 2:61846/90‑1 (MQ=255) ttCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCAttt < 1:29329/90‑1 (MQ=255) ttGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTtata > 2:163879/1‑90 (MQ=255) ttGCCATTGGCGGTGCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTTGTTtata > 1:121293/1‑90 (MQ=255) gCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATg < 1:158645/90‑1 (MQ=255) cGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCtttt > 1:29573/1‑90 (MQ=255) | ATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGATTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATATGATCATTGGTTTTATGCTCTTTT > NZ_CP009273/4419712‑4419885 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |