Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,608,161 | 0 | A | C | 53.3% | 13.6 / 20.4 | 15 | V267G (GTA→GGA) | sad | succinate‑semialdehyde dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (0/8); total (5/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCC > NZ_CP009273/1608026‑1608292 | cGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCt < 1:616088/139‑1 (MQ=255) agagCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGa > 2:565010/1‑138 (MQ=255) ttCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTc > 2:576497/1‑139 (MQ=255) cccATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCg > 1:24765/1‑139 (MQ=255) cATTTTAAAGGCTCCCGCAGCTGCCAAAAAACGTTCGGTAAA‑TGCCAAAGCAATTCCCTCTTCGATAATAAAGCTTTTTGCCGCTGCACATCCCTGCCCGGTATTCTGATAACGTCCGGCTCCCGCCGCTTTCACCGcc < 1:563768/139‑1 (MQ=255) cGCAGCTGCCAAAAAACGTTCGGAAAA‑TGCCGAAGCAATTCCCTCTTCGATAAAAAAGCTTTTTGCCGCTGCACATCCCTGTCGGGTTTTCTGATAACGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCa < 1:565010/139‑1 (MQ=255) ccAAAAACCTTTCGGAAAA‑TGCCAAAGCATTCCCCTTTTCAATAAAAAGGCTTTTTGCCGCTGCACATCCCTGCCCGGTATTCTGATACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAg < 2:362892/139‑1 (MQ=255) gTTCGGTAAA‑GGCCAAAGAAATCCCCTCTTCAATAAAAAAGCTTTTTCCCGCTGCACATCCCTGCCGGTTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATaaa < 1:11844/139‑1 (MQ=255) ttCGGTAAATTGCCAACCCAATCCCCTCTGCAAAAAAAAAGCGTTGTGCCGCGCCCCATCCCTGCCCGGTATTCTGATAACGCCCGGCTACCGCCGCTTCCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATaaa < 2:32710/139‑1 (MQ=255) aGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTa < 1:417867/55‑1 (MQ=255) aGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTa > 2:417867/1‑55 (MQ=255) aGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGcc > 1:214010/1‑138 (MQ=255) ctttttGCCGCTGCCCATCCCGGCCGGGTATTCTGATAACGCCCGGCTACCGCCGCTTTCCCCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTg < 2:202066/137‑1 (MQ=255) ctttttCCCGCTCCACATCCCTGCCCGTTATTCTAATACCTTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTg < 2:24765/137‑1 (MQ=255) cacaTCCCTGTCCGTTTTTCTGATAACTTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc < 1:183209/137‑1 (MQ=255) | CGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAA‑TGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCC > NZ_CP009273/1608026‑1608292 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |