Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 949,854 | 0 | A | C | 57.1% | 7.8 / 11.7 | 14 | G23G (GGT→GGG) | focA | formate transporter FocA |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/8); total (6/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949772‑949942 | agaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGGTTTATAGACACCCGCct > 2:463342/1‑90 (MQ=255) agaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCct > 1:666442/1‑90 (MQ=255) agaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCct > 2:559258/1‑90 (MQ=255) ttGAGATAAAACCCCCGGCGGAAACCCCCAAAAAAAAATTTTTAAGCGGTTTTTTCTTTTTTTTAAAGACCCCCGCCTCTTCGGCCACtt < 1:350564/90‑1 (MQ=255) gagaTGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTg > 1:579397/1‑90 (MQ=255) gaTGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGc > 1:432806/1‑90 (MQ=255) cccGGCGGTAATCGCCAGAAAGAAATTCTTAAGCGGTTTTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagc < 1:382425/89‑1 (MQ=255) cccGGCGGAAACCGCCAGAAAAAAAGTCTTAAGCGGTTGTTTCTTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagc < 1:645156/89‑1 (MQ=255) accGCCAGAAAAAAATTCTAAAGGGGTTTTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaa < 2:10626/88‑1 (MQ=255) aTCGCCAGATAAAAAGTCTTAAGGGGATTTTTCTTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaa < 1:607144/90‑1 (MQ=255) ttAAGGGGTTTTTTCGTTGTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCa < 1:444491/90‑1 (MQ=255) ttAAGGGGTTGTTTCTTTTTTTTAAACACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCa < 2:666442/90‑1 (MQ=255) ttGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc < 2:579397/90‑1 (MQ=255) ttATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctctct > 1:85206/1‑90 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGcagcag > 2:374987/1‑36 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 2:603753/1‑90 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 2:603874/1‑90 (MQ=255) | AGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949772‑949942 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |