| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,634,341 | 0 | C | T | 47.8% | 11.9 / 11.2 | 23 | F40F (TTC→TTT) | dtpB | dipeptide/tripeptide permease DtpB |
| Reads supporting (aligned to +/- strand): ref base C (2/10); new base T (11/0); total (13/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.73e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAAC > NZ_CP009273/3634208‑3634465 | gAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTttttttttt > 2:413149/1‑136 (MQ=255) ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTtaa > 1:366696/1‑138 (MQ=255) tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt < 1:413149/139‑1 (MQ=255) aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTCTTCGTTTAACAGCTTGGGtttt > 1:511967/1‑137 (MQ=255) tACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAGCTTGGGTTTTCg > 2:87563/1‑139 (MQ=255) caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAACTTGGATTCTCGCaa > 2:37248/1‑137 (MQ=255) acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGc < 1:87563/139‑1 (MQ=255) cacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAACAACTTGGGGTCTCCCAagag > 2:91413/1‑137 (MQ=255) acCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAg < 1:37248/139‑1 (MQ=255) ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCTTTAAAAAGCTTGGGTTTTCGCAAGAGCAggggttt > 2:56427/1‑134 (MQ=255) gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAACTTGGGTTTTCCCAAGAGGAGGGTTTTTTCACt > 1:234407/1‑139 (MQ=255) gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTTGGTTCTCGCAAGAGCAGGCTTTTGTCAcc > 1:122672/1‑138 (MQ=255) tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGATTCTCGCAAGAGCAGGGTTTTGTTACttt > 2:265027/1‑139 (MQ=255) cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCTTTAAAAAGAGTTGGTTTTCGCCAAAGCAGGGCTTTTTGACttttt > 1:97859/1‑138 (MQ=255) agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTTAACAGCCTTGGTTCTCGCAAGAGCAGGGTTTTTTCCCtttttg > 2:413239/1‑137 (MQ=255) gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtg < 1:470300/138‑1 (MQ=255) gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtg < 1:91413/138‑1 (MQ=255) aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcg < 1:265027/139‑1 (MQ=255) cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt < 1:56427/138‑1 (MQ=255) ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATgg < 2:350440/135‑1 (MQ=255) ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGc < 1:36742/136‑1 (MQ=255) cGTACTGGCGGTTTTTTTCGTTAAACAGCTTGGGTTTTCGCCAGAGCAGGCCTTTTTCCCTTTTTGTGCTTTTTCTTCGCTGGTCTATGGCCTCATTTCCCTTTGTGGCTTTTTTTGCGACCACCCTCTGGGGGCCcaa > 1:315326/1‑139 (MQ=255) cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGAACCaaa < 2:174209/139‑1 (MQ=255) gTACTGGCGGTTTTTTTCGTTAAAAAACTTGGGTTCTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGCTTTTTCTGCGCTGGTCTATGGCCTCATTTTCCTTTGCGGCTTTTTTCGCGGCCACCCCGTGGGGGCCCAAc > 1:325915/1‑139 (MQ=255) | GAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAAC > NZ_CP009273/3634208‑3634465 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |