Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 949,854 | 0 | A | C | 55.0% | 7.1 / 17.9 | 20 | G23G (GGT→GGG) | focA | formate transporter FocA |
Reads supporting (aligned to +/- strand): ref base A (7/2); new base C (0/11); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.64e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949765‑949942 | atatAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGaca > 2:810812/1‑90 (MQ=255) gaagaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGACTTAAGCGGATGTTTTGGTGGTTTATAGACACCCGc > 1:560499/1‑90 (MQ=255) aaagaCAAATGCGATTGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGTTTTATAGACCCCCGc < 1:387502/89‑1 (MQ=255) agaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCct > 2:77341/1‑90 (MQ=255) agaCGAATGCGATTGAGATGAAAACACCGGCGGGAATCGGCAGATAGAAAGGCTTAAGCGGATGTTTCGTTGCTTTATAGACAACCGCct > 2:672021/1‑90 (MQ=255) aTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCgg > 2:659766/1‑90 (MQ=255) aTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCgg > 1:861225/1‑90 (MQ=255) tgggATTGAGATAAAACCCCCGGGGTTAACCCCCAGATAGAAAGTTTTAGGGGGTTTTTTGTTTGTTTTATACCCCCCCGCCTCTTCGGc < 1:460003/87‑1 (MQ=255) tGCGATTGAGAAAAAACCCCCGGGGGTAACCGCCAGATAGAAAGTCTTAAGCGGTTTTTTCGTTGTTTTATAGACCCCCGCCTCTTCGGc < 1:277895/90‑1 (MQ=255) gCGATGGAGATAAACACCCGGGGGGTAACCCCCAAAAAAAAATTCTAAGGCGGTTGTTTCGTTGTTTTAAAGACCCCCGCCTCTTCGGcc < 1:794106/90‑1 (MQ=255) cggcggTAACCCCCAAAAAGAAAGTTTTAAGGGGTTTTTTCGTTGTTTTAAACCCCCCCGCCTCTTCGGCCACTTTGGCCATTGcagcag < 1:939424/90‑1 (MQ=255) gATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCaa > 2:710121/1‑90 (MQ=255) taaaGCGGTTGTTTCGTTGTTTTATAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCa < 2:205585/88‑1 (MQ=255) ggATGTTTTTTTGTTTTAGAGCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTc < 1:828056/90‑1 (MQ=255) atttttCCTTCCTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcac < 2:42887/87‑1 (MQ=255) ttttCTTTGTTTTAAAAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCacac < 1:797711/89‑1 (MQ=255) tCGTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTaa < 2:861225/90‑1 (MQ=255) ttGTTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc < 1:77341/90‑1 (MQ=255) tGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcac > 2:177768/1‑81 (MQ=255) tGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcac < 1:177768/81‑1 (MQ=255) tttATACCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctctc < 2:678395/90‑1 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 1:261161/1‑90 (MQ=255) cacCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTaa > 1:25718/1‑90 (MQ=255) | ATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAA > NZ_CP009273/949765‑949942 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |