| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,197,867 | 0 | T | G | 42.4% | 47.2 / 13.0 | 33 | A47A (GCA→GCC) | cirA | ferric iron‑catecholate outer membrane transporter |
| Reads supporting (aligned to +/- strand): ref base T (4/15); new base G (14/0); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.06e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACA > REL606/2197795‑2197942 | gAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCa < 2:503529/76‑1 (MQ=255) ccTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAg < 1:696159/76‑1 (MQ=255) tCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGAtt > 1:707528/1‑76 (MQ=255) ttCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATtttttt > 2:1147219/1‑76 (MQ=255) tGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTTAAATTTTTTTTCa > 2:875483/1‑76 (MQ=255) tGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCa > 1:1120080/1‑76 (MQ=255) ccGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGa < 2:1298020/76‑1 (MQ=255) nggTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTAAGATTTTTTTCCCCCGaa > 1:272141/2‑76 (MQ=255) cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa < 2:204210/76‑1 (MQ=255) cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa < 1:609436/76‑1 (MQ=255) cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa < 1:1120809/76‑1 (MQ=255) ggTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAg < 1:741296/76‑1 (MQ=255) tttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCCCGGCAGaag > 2:40399/1‑74 (MQ=255) cGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAg < 2:459062/76‑1 (MQ=255) cGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTAAAATTTTTTTTCACAGGAAATGCAg > 1:1293601/1‑76 (MQ=255) cGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAGATTTTTTTCCACGGAAGATGCAg > 2:1067904/1‑76 (MQ=255) gCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGGATCTTTAAAATTTTGTTTTCCGGAAGATGCAgg > 1:236852/1‑75 (MQ=255) cTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTg < 1:1047024/76‑1 (MQ=255) cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCACGGAAGATGCAGTGACa > 2:1053692/1‑76 (MQ=255) tCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCCCGGAAGATGCAGTGACAAACa > 2:1392809/1‑76 (MQ=255) tCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCg < 1:1340335/76‑1 (MQ=255) ccTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAGATTTTTTTTCACGGAAGATGCAGTGACAACCATCGt > 2:75831/1‑76 (MQ=255) cTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGtt < 2:150309/76‑1 (MQ=255) gggTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTc < 2:549215/76‑1 (MQ=255) aaTGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACGGAAGATGCAGTGACAAACAACGTTTCGccc > 2:902314/1‑75 (MQ=255) aaTGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTCCACGGAAGATGCAGTGACAACCATCGTTTTGccc > 1:1183445/1‑75 (MQ=255) tGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCcatc < 1:48856/76‑1 (MQ=255) gACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGAAGATGCAGTGACAAACCTCGTTTCGCcctca > 2:1380578/1‑76 (MQ=255) gATACTGGCAGGTGCGTCTTTAAAATTTTTTTTCTCGGAAGATGCAGAGACAACCATCGTTTTGCCATCATCATCg > 2:510284/1‑76 (MQ=255) gATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCg < 1:635911/76‑1 (MQ=255) gATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCg < 2:744146/76‑1 (MQ=255) aCTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGAAGATGCAGGGACAACCATCGTTTTTCCCTCATCATCGAcc > 2:47537/1‑76 (MQ=255) cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCg > 1:1434191/1‑76 (MQ=255) gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTa > 2:1330722/1‑76 (MQ=255) tGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAAcaca > 2:1084123/1‑76 (MQ=255) tGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAAcaca < 1:560327/76‑1 (MQ=255) tGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAAcaca > 2:577826/1‑76 (MQ=255) | GAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAACACA > REL606/2197795‑2197942 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |