| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,075,148 | 1 | . | T | 27.9% | 75.3 / 51.9 | 43 | coding (327/1605 nt) | ydhS | conserved hypothetical protein with FAD/NAD(P)‑binding domain |
| Reads supporting (aligned to +/- strand): ref base . (0/31); new base T (0/12); total (0/43) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGG > minE/1075093‑1075194 | tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGTCGAATATTTCCGCGATCAA‑TTTTTa < 1:121895/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:418804/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:238358/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:238700/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:257555/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1022814/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:340480/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:340768/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:351425/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:264926/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:468622/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:510532/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:61793/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:628618/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:737258/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:793625/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:886762/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:990701/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1024156/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1180819/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1036113/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1158495/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1149476/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1059421/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1074701/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1076456/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1093816/62‑1 (MQ=255)tGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1114297/62‑1 (MQ=255) gCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:914246/61‑1 (MQ=255) aGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:498930/50‑1 (MQ=255) ttACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTa < 1:1162180/45‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:444274/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:425565/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:1147012/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:678957/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:687458/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:722176/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:1062467/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:1171542/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:125967/62‑1 (MQ=255) tatTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTAGCAGTgg < 1:816863/62‑1 (MQ=255) atTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:149443/61‑1 (MQ=255) atTTCCGCGATCAATTTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTgg < 1:366360/61‑1 (MQ=255) | TGCACGATCGTCAGTTTTTACCGCGAATTCTGCTGGGCGAATATTTCCGCGATCAA‑TTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGG > minE/1075093‑1075194 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |