Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,171,107 | 0 | C | A | 46.7% | 97.3 / 32.4 | 60 | E55* (GAG→TAG) | nudF | ADP‑ribose diphosphatase |
Reads supporting (aligned to +/- strand): ref base C (29/3); new base A (0/28); total (29/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.36e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCC > NZ_CP009273/3171025‑3171191 | gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt > 2:1052965/1‑90 (MQ=255) gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt > 1:983144/1‑90 (MQ=255) gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc > 1:790017/1‑90 (MQ=255) gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc > 2:1898441/1‑90 (MQ=255) aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc > 2:962206/1‑90 (MQ=255) aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGACGTGACCGCGCTCAAAAATTTCccgc > 1:478912/1‑90 (MQ=255) aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGACGTGACCGCGCTCAAAAATTTCccgc > 1:261168/1‑90 (MQ=255) cGAATCTGCTCAATCAGCACACCTTCACCCCGCCCGGGGCAAAGGGG‑TAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCcgccgc < 1:1501867/90‑1 (MQ=255) cGAATCTGCTAAATCAGCACAACTCCACCCCGCACTGGGTCAAGGGG‑TAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTCCCcgccgc < 1:307083/90‑1 (MQ=255) aTCTGCTCAATCAGCACAACTTCATCACGCACGGGGTAAAGGGG‑TAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCAcc < 2:2118238/90‑1 (MQ=255) aTCTGCTCAATCAGCACAACTTCACCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCAcc < 2:1993912/90‑1 (MQ=255) gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt > 2:810341/1‑90 (MQ=255) tCAATCAGCACACCTTCATCCCGCACGGGGCAAAGGGG‑TAGCAAGACTGGGGCGTGGCGGCGCTAAAAAATTTCCCGCCGCACCTCATGa < 2:49176/90‑1 (MQ=255) aGCACAACTTCACCACGCACGGGGCCAAGGGG‑TAGAAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCAtt < 1:510481/90‑1 (MQ=255) cacaACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTg > 1:2014039/1‑90 (MQ=255) caACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc < 1:961725/90‑1 (MQ=255) caACTCCACCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc < 1:757852/90‑1 (MQ=255) aCTTCATCCCGAACGGGGTAAAAGGG‑AAGAAAGATGGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCg < 2:689612/90‑1 (MQ=255) cTTCATCCCGCACTGGGCAAAGGGG‑TAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt < 1:1514768/90‑1 (MQ=255) cTTCACCCCGCACTGGGTAAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt < 2:684869/90‑1 (MQ=255) cTTCACCACGCACGGGGTCAAGGGG‑TAGCAAGCCTGCGGGGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGt < 1:800722/90‑1 (MQ=255) ttCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt > 1:233345/1‑90 (MQ=255) tcatcaCGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg > 1:1578106/1‑90 (MQ=255) tcatcaCGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCCCCGCACCTCATGACTCATTTGCCCGTTg < 1:1451163/90‑1 (MQ=255) gcGCACTGGGCAAAAGGG‑AAGCAGGTCGGGGGGGGGCCCGCGTAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg < 1:1687332/89‑1 (MQ=255) ccGCCCTGGGTAAAAGGG‑TAGCAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg < 2:1182056/89‑1 (MQ=255) aCGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg > 1:120544/1‑90 (MQ=255) gCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc > 2:1017529/1‑90 (MQ=255) gCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc > 1:1343672/1‑90 (MQ=255) ccACTGGGTCAAGGGG‑AACAAAGACCGCGGCGTGCCCGCGCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc < 1:1160383/89‑1 (MQ=255) cACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg > 1:384154/1‑90 (MQ=255) cACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg > 2:800722/1‑90 (MQ=255) ccTGGGTCAAGGGG‑AAGCAAGCCGGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 1:1398068/89‑1 (MQ=255) aCTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 1:545768/90‑1 (MQ=255) aCTGGGTCAAAGGG‑AAGAAAGACTGCGACGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:261168/90‑1 (MQ=255) aCTGGGTAAAAGGG‑TAGCAAGACTGCGACGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa < 2:478912/90‑1 (MQ=255) gTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg > 1:287539/1‑90 (MQ=255) taaaaGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa < 2:1721213/88‑1 (MQ=255) cAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGac > 2:1687332/1‑89 (MQ=255) cAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa > 1:354992/1‑90 (MQ=255) aaagggg‑TAGCAAGCCTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa < 2:695306/86‑1 (MQ=255) aagggg‑TAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:2156584/87‑1 (MQ=255) aaaGGG‑TAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:1792756/90‑1 (MQ=255) aaaGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:384154/90‑1 (MQ=255) aaaGGG‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:911057/90‑1 (MQ=255) aaaGGG‑TAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:1301006/90‑1 (MQ=255) aaGGGTTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 2:287539/90‑1 (MQ=255) aaGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt > 2:1060291/1‑90 (MQ=255) aGGG‑TAGCAAGACTGCGGCGTGCCCGCGAAAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc < 1:1795654/90‑1 (MQ=255) ggg‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt > 2:2129208/1‑90 (MQ=255) ggg‑TAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt < 2:1613849/90‑1 (MQ=255) g‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat > 2:997198/1‑90 (MQ=255) tAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata > 1:1225020/1‑90 (MQ=255) tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa > 2:200923/1‑90 (MQ=255) tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa > 1:271497/1‑90 (MQ=255) gcggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt > 2:1174792/1‑90 (MQ=255) gTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa < 1:166661/69‑1 (MQ=255) gTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa > 2:166661/1‑69 (MQ=255) gcgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc > 1:1558999/1‑90 (MQ=255) gcgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc > 1:1948761/1‑90 (MQ=255) | GCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGG‑TAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCC > NZ_CP009273/3171025‑3171191 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |