Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 5,159 | 0 | G | A | 54.5% | 5.0 / 27.6 | 22 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (9/1); new base A (0/12); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA > NZ_CP009273/5074‑5243 | tttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa > 1:1224376/1‑90 (MQ=255) aaaaaaaaGGAAAGGAAAAAAGAAAAAATTCAAAAAAAAGGGGGAAACTTAGAAATAAGGTTGGGGGAAAAAAACACCCGCCGTTCTCAt < 2:1224376/87‑1 (MQ=255) aaTGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAATTTAGAGATTGGGATTGGGGAGAATAACAACCGCCGTTCTCATCGAGTa < 2:477906/90‑1 (MQ=255) aaTGAAAGGAAAAAAGAAAAAATTCAAAAAAAAGGGGTTAACTAAGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTa < 1:289989/90‑1 (MQ=255) agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc > 2:728296/1‑90 (MQ=255) gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc > 2:781609/1‑90 (MQ=255) aataatTGTGGTAATTAAGAATTTAGGATGGGGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc < 1:717181/90‑1 (MQ=255) aaaaTGGGTTAACTTAGAGATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:1193471/89‑1 (MQ=255) aaaaTGCGTTAACTAAGAGATTAGGATGGGGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:521874/89‑1 (MQ=255) aaaaTGCGGAAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:347689/89‑1 (MQ=255) aaTGGGGTAACTTAAAAATTAGGATGGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg < 1:948011/90‑1 (MQ=255) gttAACTTAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:615141/88‑1 (MQ=255) ggTAATTTAGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:467948/90‑1 (MQ=255) ggTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:401995/90‑1 (MQ=255) gatAATTTAGAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:1186711/88‑1 (MQ=255) gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt > 2:1119551/1‑90 (MQ=255) aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAATGGGCAATGATAAAAGGAGTAACCTGTGa > 2:931451/1‑90 (MQ=255) aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa > 2:167006/1‑90 (MQ=255) aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa > 1:652831/1‑90 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa > 1:904460/1‑69 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:997777/1‑90 (MQ=255) gCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 2:904460/69‑1 (MQ=255) | TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA > NZ_CP009273/5074‑5243 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |