| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 276,800 | 0 | A | C | 56.2% | 37.0 / 40.0 | 32 | Y299D (TAT→GAT) | yagA | IS481‑like element ISEc19 family transposase |
| Reads supporting (aligned to +/- strand): ref base A (10/4); new base C (0/18); total (10/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.55e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGA > NZ_CP009273/276721‑276874 | tCCATTTTCCTGACCATCCCCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGccc < 2:46638/90‑1 (MQ=255) aCTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCGGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGt > 1:439430/1‑90 (MQ=255) tttCCTGCCCATCCCCCCTTCATCTTTTTCCGGGGGCTTTGTTTTGCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTAc < 2:24024/90‑1 (MQ=255) tttCCTGCCCATCCCCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTAc < 1:91609/90‑1 (MQ=255) tttCCTGCCCACCCCCCCTTCACCGTTTCCGGGGGGCGTTGTGTTCCCGCTGTCCTGCCGCGCTGCCGGCTGATCCCGCGAGCCGGGTAc < 2:214476/90‑1 (MQ=255) tttCCTGACCATCCCCCCTTCATCTTATCCGGGGGGGGTTGTGTTGCCGCTGTACTGCCGCGCTGCCGGATGATACCGCGACCCGGGTAc < 2:214782/90‑1 (MQ=255) gACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCat > 1:253193/1‑90 (MQ=255) tCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCa > 1:190052/1‑90 (MQ=255) tCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCa > 1:245190/1‑90 (MQ=255) tCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCa > 1:412408/1‑90 (MQ=255) aCCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAgc > 1:300030/1‑90 (MQ=255) ctccATCGTATCCGGGGGGCGTTGTTTTGCCGCTGTCCTGCCGCGCTGGGGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCt < 2:320551/87‑1 (MQ=255) gggggggGTTGTGTTGCCGCTGTCCTCCCGCGCGGGGGGCTGATCCCGCGGGCCGGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTc < 2:439430/89‑1 (MQ=255) ggggggCGTTGTGTTGCCGCTGTCCTGCCGCGCTGCCGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTc < 1:491703/89‑1 (MQ=255) cGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGGGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTc < 1:99730/90‑1 (MQ=255) ggggggTTGTGTTCCCGCTGTCCTGCCGCGCGGGCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 2:460970/90‑1 (MQ=255) ggggCTTTGTGTTGCCGTTGTCCTGCCGCGCTGCCGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:202397/90‑1 (MQ=255) ggggCTTTGTGTTGCCGCTGTCCGGCCGCGCTGGCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 2:55130/90‑1 (MQ=255) ggggCTTTGTGTTGCCGCTGTACTGCCGCGCTGGCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:227664/90‑1 (MQ=255) ggggCGTTGTTTGCCCCCTGTCCCGCCGCGCTGCGGGCTGATCCCGCGAGCCCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 2:152608/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:129322/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTCCGGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:28258/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTCCCGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:218614/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa > 1:460970/1‑90 (MQ=255) ggggCGTTGTGTTCCCGCTTTCCTCCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:355398/90‑1 (MQ=255) ggttttGTGTGCCCGCTGTCCCGCCGCGCTGTGGTTGGATCCCGCGGCCCCGGACCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAgg < 2:407039/86‑1 (MQ=255) gtTGCCTCTGTACTCCCGCGCTGCCGGCTGATCCCGCGCGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGAc < 1:478965/90‑1 (MQ=255) tGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa > 2:424033/1‑78 (MQ=255) tGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:424033/78‑1 (MQ=255) gCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGt > 2:246565/1‑90 (MQ=255) tgtgCTGCCGCGCGGCGGGCTGATCCCGCGGGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGcc < 1:310999/90‑1 (MQ=255) gcgcTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGa > 2:110080/1‑90 (MQ=255) | TCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGA > NZ_CP009273/276721‑276874 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |