Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,171,107 | 0 | C | A | 55.2% | 15.4 / 42.1 | 29 | E55* (GAG→TAG) | nudF | ADP‑ribose diphosphatase |
Reads supporting (aligned to +/- strand): ref base C (11/2); new base A (0/16); total (11/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGC > NZ_CP009273/3171028‑3171190 | gCACTCCAAATCTGCCCAATCAGCACATCTTAAACAC‑CACTGGGTTCAAAGGGTAGCAAGACGGGGGGGTGCCCGCGC‑AAAAAATTTccc < 1:167252/90‑1 (MQ=255) caaaTCTGCTCAACCACCACAACTTCATCCCCCACGGGG‑TCAAGGGGTAGCAAGACTGGGGGGTGACCGCCCTAAAAATTTCCCcgccgc < 1:831261/88‑1 (MQ=255) cGAATCTGCTCAATCACCACACCTTCATCACGCCCGGGG‑TAAAGGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCcgccgc < 1:725057/90‑1 (MQ=255) aaTCTGCTCAACCCGCCCAACTTCACCCCGCACTGGG‑TAAAAGGGCAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCAc < 2:445844/90‑1 (MQ=255) gCTCAATCAGCACAACTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt > 2:831261/1‑90 (MQ=255) gCTCAATCAGCACAACTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt > 2:725057/1‑90 (MQ=255) gCTCAATCAGCACAACTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt > 1:69374/1‑90 (MQ=255) caccTTCACCACGCACGGGG‑TAAAAGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc < 2:71129/87‑1 (MQ=255) cTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt < 2:700441/90‑1 (MQ=255) tcccGCACGGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt < 2:661419/87‑1 (MQ=255) ccgggg‑TCAAAGGGTAGCAAGACGGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:538947/87‑1 (MQ=255) aCTGGG‑TCAAAGGGAAGCAAGACTGCGGCGTGCCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:489323/90‑1 (MQ=255) aCTGGG‑TAAAAGGGAAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:183367/90‑1 (MQ=255) aaaGGGTAGCAAGACTGGGGCGTGCCCGCCCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:81605/90‑1 (MQ=255) aaaGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:103538/90‑1 (MQ=255) aaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa > 1:59176/1‑89 (MQ=255) aaGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:59176/89‑1 (MQ=255) aGGGTAGCAAGACTGCGGCGTGCCCGCTCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc < 1:337004/90‑1 (MQ=255) gTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat > 1:616657/1‑90 (MQ=255) gTAGCAAGACTGCGGCGTGACCGCCCTCAAAAATTTCCCGCCGCACCTCATGGCTCATTTGCCCGTTGAATGGACGATGGCGAAATCtat > 1:595066/1‑90 (MQ=255) aGCAAGCCTGGGGCGTGCCCGCCCAAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa < 1:514885/90‑1 (MQ=255) aGCAAGCCGGCGGGGGGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa < 2:595066/90‑1 (MQ=255) aCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt > 2:526722/1‑90 (MQ=255) aCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt > 2:692935/1‑90 (MQ=255) aCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt > 2:239138/1‑90 (MQ=255) cTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa > 1:725622/1‑90 (MQ=255) tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa > 2:823869/1‑90 (MQ=255) gcggcgTGCCCGCCAAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt < 1:532842/90‑1 (MQ=255) cgcgCCAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGc < 2:69374/90‑1 (MQ=255) | GCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGG‑TCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGC > NZ_CP009273/3171028‑3171190 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |