| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 949,854 | 0 | A | C | 60.0% | 10.9 / 16.2 | 21 | G23G (GGT→GGG) | focA | formate transporter FocA |
| Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/12); total (6/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.22e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949777‑949942 | aaTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCg > 1:573476/1‑90 (MQ=255) aTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCgg > 1:1140322/1‑90 (MQ=255) gCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGGTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGcc > 2:683112/1‑90 (MQ=255) ttaagaTAAAACCCCCGGGGGTAACCGCAAGATAAAAAGTTTTAAGGGGTTTTTTCGTTGCTTTATAGCCCCCCGCCTCTTCGGCCACtt < 1:1263090/87‑1 (MQ=255) ataaaaaCCCCGGCGGAAAACGCCAGATAGAAACTCTTACGCGGATGTTTGGTTGCTTTATAGACTCCCGCCTCTTCGGCCACTTTGGcc < 1:975239/87‑1 (MQ=255) cgggggAAACCGCAAGAAAAAAAGTCTTAGGGGGTTGTTTCGTTTTTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGcagcag < 2:239949/90‑1 (MQ=255) cggcggTAACCGCCAGATAGAAAGTCTTAAGCGGATTTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGcagcag < 1:91991/90‑1 (MQ=255) ggcggTAATCGCAAGAAAAAAAGTCTTAAGGGGTTGTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg < 1:613514/90‑1 (MQ=255) gcggTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCTTCTTCGGCCACTTTGGCCATTGCAGCAGGa < 1:288981/90‑1 (MQ=255) gcggTAACCGCCAGATAGAAAGTCTTAAGGGGTTTTTTCGTTGTTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGa < 1:683112/90‑1 (MQ=255) aGAAAGTCTAAAGCGGATGTTTCGTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAg < 1:244152/90‑1 (MQ=255) cTTAAGGGGTTTTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTc < 2:348661/90‑1 (MQ=255) ttAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCa > 1:348661/1‑90 (MQ=255) tAAGCGGATGTTTCGTTGCTTTATAGACCCCCGCCTCTTAGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAg < 1:13679/90‑1 (MQ=255) ggATTTTTGTTTGTTTTAAAAACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTc < 2:1083450/90‑1 (MQ=255) ttGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc < 1:1074303/90‑1 (MQ=255) gCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc > 1:1226507/1‑90 (MQ=255) ttttAAAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGACAGCTTTCACACTAActctct < 2:475089/89‑1 (MQ=255) ttttAAACACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctct < 2:75103/89‑1 (MQ=255) tttATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGCCAGCTTTCACACTAActctctc > 1:833498/1‑90 (MQ=255) ttAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctctct < 2:56038/90‑1 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 1:6243/1‑90 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 2:828896/1‑90 (MQ=255) | AATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949777‑949942 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |