Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,611,723 | 0 | G | A | 41.7% | 22.9 / 11.2 | 24 | K196K (AAG→AAA) | nikE | nickel import ATP‑binding protein NikE |
Reads supporting (aligned to +/- strand): ref base G (3/11); new base A (10/0); total (13/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAG > NZ_CP009273/3611588‑3611849 | gTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTTGCCCGCTGAAAAAACTa > 1:66477/1‑139 (MQ=255) aGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaac < 1:388358/139‑1 (MQ=255) aGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaac < 2:123347/139‑1 (MQ=255) cggcCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAATCTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGGGTTATGCCTCCTGTAAAAACCacaacaacc > 2:402960/1‑137 (MQ=255) gcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCCCCGCCTGCCTGTTCATCACCCa < 1:140958/139‑1 (MQ=255) gcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCa < 1:56693/139‑1 (MQ=255) gTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGcc < 1:16095/139‑1 (MQ=255) actactGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCCTTCGCCCGCTTAAAAAACTACAACAACAACTTGGCACCGCCCGCCTGTTCATCACCCCCGACTTAAGCCCGGTGGAAAGc > 2:244790/1‑139 (MQ=255) ttCTGGATGACGCTGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCCCGACTTACGCCTGGTAGAACGCTTTTGCCa < 2:163249/139‑1 (MQ=255) ttCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCa < 1:236505/139‑1 (MQ=255) tCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAACTACAACAACAGTTTGGGACCGCCCGCCCGTTCCTTACCCCCCACTTATGCCTGGTAGAAAGCTTTTGcccg > 2:387662/1‑137 (MQ=255) tCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGGCATTTGCCTCCTGAAAAAAATACAACAACAGATTGGCGCCGCCCGCCTGTTTCTCACCCCCCACCTACCCCCGGGAGGACGCCTTTGcccg > 2:654568/1‑137 (MQ=255) gagacagTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa < 1:579851/133‑1 (MQ=255) aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa < 2:43650/138‑1 (MQ=255) aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTCGCCTGCTGAAAAAACTACAACAAAAATTTTGCACCGCCCGCCTGTTCATCACCCACCACCTACGCCCGGTAGAAAGCTTTTGCCAGCgggggat > 1:60893/1‑135 (MQ=255) aTCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGAc < 1:654568/138‑1 (MQ=255) aTCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAAATAAAACAACAGTTTGGGACCGCCCGCCCGTTCATTACCCACGACCTACGCCCGGTAGAACGCTTTTTCCCGCGGGGAAATGTTTTGGGCAACGGAc > 2:560259/1‑138 (MQ=255) cGTGTTACAGGCGGGGGTCATTCGCCTGCTGAAAAAACTACAAAAAAAGTTTGGCACCGCCTGCCTGTTCATCACCCCCGCCTTACGCCCGGGAGAACGCTTTTGCCAGCGGGGAATGGGTATGGACAACAGACaaaac > 1:128592/1‑137 (MQ=255) tgtTACAGGCGGGTGTCATTTTCCTGCTGAAAAAACTACAACAACAGTTTGGCACCGCCCGCCTGTTCATCACCCACGACTTACGCCCGGTAGAACGCTTTTTCCAGCGGGGAATGGTTATGGACAAAGGACAAAtcgt > 1:140024/1‑139 (MQ=255) tgtTACAGGCGGGTGTCATTCGCCCGCTGAAAAAGCTACAACCAAAGATTGGCACCGCCCGCCTGTTCATCAACCCCGACTTACGCCTGGTAGAACGCTTTTTCCAGCGGGGAATGGTTATGGACAACGGACAAAacgg > 2:638169/1‑138 (MQ=255) ggCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAAcc < 2:60893/139‑1 (MQ=255) ggTGTTATTCTCCCGCTGGAAAAACAACAAAAAAAATTTGGGACCGCCCGCCTGTTCATCACCCACGACTTTTGCCCGGTAGAAAGCTTTTTCCACCGGGTAATGGGTTTGGACAACGGACAAAACGGCGAAACCCAgg > 2:496637/1‑139 (MQ=255) cGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGtggtgg > 1:633832/1‑136 (MQ=255) ctgctgAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGgagag > 1:379834/1‑138 (MQ=255) | GTTAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAG > NZ_CP009273/3611588‑3611849 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |