Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,549,566 | 0 | T | C | 54.5% | 9.6 / 27.2 | 23 | S215G (AGT→GGT) | rtcA | RNA 3'‑terminal phosphate cyclase |
Reads supporting (aligned to +/- strand): ref base T (5/5); new base C (0/12); total (5/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.57e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NZ_CP009273/3549434‑3549702 | ctcTGCCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGGAAAAC‑CCCCCCGc < 1:422051/137‑1 (MQ=255) cccGACGACAAAAAGGCGTTCGGTGATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCAGGTTATAAATATTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCa < 2:31512/138‑1 (MQ=255) cgacgaAAAAAAAGGGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACg < 1:605764/139‑1 (MQ=255) gacgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGc > 1:438263/1‑139 (MQ=255) aaaaaaagGCTTTCGGTGATTTTTTCACTTGCGACTACAAGCGAGGCGGTATCCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTTTTCATGCAGGAAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAg < 2:362778/137‑1 (MQ=255) ggcttgcGGTAATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGCCGCGCGGAAGTTTATAAATTTTCTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg < 1:208175/133‑1 (MQ=255) attttttCACTTCCGACTCAAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGAAGTTTATGAATATTCTGTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac < 2:564871/136‑1 (MQ=255) atatTTTCATTTCCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca < 1:2799/138‑1 (MQ=255) ttttCTCTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTCCATGCAGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 1:263293/138‑1 (MQ=255) ggcgggTTTTCCCCGCCCCGTGGCCGCCCGGCAGGTTATGAATTTTCTTTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 2:597597/135‑1 (MQ=255) agtcGGTATTCCCCGCCCCCTGGCCGCGCGGGAGTTTATGAATCTTCTGTAAAGCCAGGGAAAAA‑CGCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt < 1:224907/136‑1 (MQ=255) gTATTCCCCGGCCCCTGGCCGCGCGGCAGGTTATAAATATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca < 1:344939/139‑1 (MQ=255) gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca < 1:603268/139‑1 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 1:358848/100‑1 (MQ=255) ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa > 2:358848/1‑100 (MQ=255) cgcgcgGCAGTTTATAAATATTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc < 1:194681/139‑1 (MQ=255) aaaTATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt < 2:209687/138‑1 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca < 1:322633/99‑1 (MQ=255) tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca > 2:322633/1‑99 (MQ=255) gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa < 1:229543/62‑1 (MQ=255) gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa > 2:229543/1‑62 (MQ=255) gAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga > 1:171096/1‑139 (MQ=255) a‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc < 1:16275/139‑1 (MQ=255) | CACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC > NZ_CP009273/3549434‑3549702 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |