Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 726,200 | 0 | A | C | 54.2% | 18.3 / 25.3 | 24 | T388P (ACC→CCC) | rhsC | rhs element protein RhsC |
Reads supporting (aligned to +/- strand): ref base A (7/4); new base C (0/13); total (7/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACG > NZ_CP009273/726111‑726284 | gTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatca > 1:593896/1‑90 (MQ=38) ttACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac < 1:354817/90‑1 (MQ=40) tACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCatcaccatcacc > 1:180121/1‑90 (MQ=38) gCGAGGGGCGGGTGACAGAACAGCAAACCCC‑GGCAGGTTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCACCACCCCCGACAGcc < 2:593896/90‑1 (MQ=11) cGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCt > 2:370954/1‑90 (MQ=38) ggcgggTGAAAGAACAGCAAACCCC‑GGCAGGTTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCCCCATCACCGACAGCCTGAAcc < 1:306985/90‑1 (MQ=14) gaaaGAACAGCAAACCCC‑GGCGGGCTAAATTTACACGTACCAGTATGAGAAGGCCCGCATCCCCATCCCCGACAGCCTGAACCGCCGTGa < 2:555970/87‑1 (MQ=21) gACAGAACAGCTAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCCCCGACAGCCTGAACCGCCGTGa < 1:671656/90‑1 (MQ=255) aaCAAAACACCCAACCCC‑GGCAGGCTAAATTAACCCGTATAAGTATAAGAGGCCCCGCACCCCCACCCCCGACGGCCTGAACCGCCGTGa < 2:486821/89‑1 (MQ=255) acagaGCTAACCCC‑GGGGGGCTTAATTTACACGTATCAGTATGAGAAGGCCCGCACCCCCATCACCGACAGCCTGAACCGCCGTGAAGTc < 1:271567/86‑1 (MQ=14) caaaCCCGGGCGGCCTTAAGTAACACGTATCAGTAGGAGAAGGCCCGCACCCCCACCCCCGACAGCCTAACCCGCCGTGAAGTCCTGcac < 2:618496/89‑1 (MQ=255) cccc‑GGCGGGTTTAAGTTACACGTATAAGTAGGAGAAGGCCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc < 2:372676/89‑1 (MQ=255) cccc‑GGCAGGCTAAAGTTACCCGTATCAGTATGAGAAGGCCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc < 2:193695/89‑1 (MQ=14) cccc‑GGAAGGTTTAAGTTACACGTACCAGAAGGAGAAGGCCCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc < 1:482229/89‑1 (MQ=14) ttAAGTTACTCGTATCAGTAGGAGAGGGCCCGCACCCCCACCCCCGACGGCCGAACCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGc < 2:46693/90‑1 (MQ=255) ttAAGTTACACGTATAAGTATGAAAAGGCCCGCATCCCCACCACCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGc < 1:717923/90‑1 (MQ=255) acacGTATCAGGATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGa > 1:788691/1‑90 (MQ=255) aTCAGTATGAGAAGGCCCGCACCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg < 1:472450/90‑1 (MQ=17) tGAGAAGGCCCCCACCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaa < 1:370954/90‑1 (MQ=17) agaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACg > 2:287570/1‑56 (MQ=255) agaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACg < 1:287570/56‑1 (MQ=255) agaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACg > 2:378270/1‑56 (MQ=255) agaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACg < 1:378270/56‑1 (MQ=255) aaGGCCCCCACCCCCACCCCCGACAGCCTGACCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAg < 2:179355/90‑1 (MQ=17) ggACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGa > 1:683458/1‑90 (MQ=255) gACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGaa > 1:245461/1‑90 (MQ=255) gCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg > 1:223780/1‑90 (MQ=255) | GTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACG > NZ_CP009273/726111‑726284 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |