Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,436,636 | 0 | A | C | 45.6% | 78.5 / 22.5 | 57 | V274G (GTA→GGA) | secY | preprotein translocase subunit SecY |
Reads supporting (aligned to +/- strand): ref base A (25/6); new base C (0/26); total (25/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.52e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA > NZ_CP009273/3436502‑3436773 | acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc > 1:251101/1‑139 (MQ=255) acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc > 2:691538/1‑139 (MQ=255) acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc > 2:314006/1‑139 (MQ=255) acacATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc > 1:435693/1‑139 (MQ=255) aaaGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCa > 2:326491/1‑138 (MQ=255) aGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCatat > 1:599355/1‑139 (MQ=255) gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc > 1:551540/1‑139 (MQ=255) gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCAc > 2:372084/1‑139 (MQ=255) cccAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACttt > 1:408478/1‑138 (MQ=255) acagCTGAAAATACAGCAAAATTGTTGTCAGCCAGTCCCAA‑CCAGTCCCCCCCCCGACCCAGGCCGCGTTGGCCGCCGGGAACAAAATAATACTGGAAGCAAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAg < 2:903521/136‑1 (MQ=255) ggaaaaTACAGCGAATTTGTTGACAGCC‑GTTCCAACCCAGTCCCCCCCCCGACCCATGCCGCGATGGCCGCCGGAAACAAAATAATACTGGAGGCAAAGATTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTa < 2:245961/136‑1 (MQ=255) cccAAATACGGCAAAATGGTTGTCAGCCAGTTCCAC‑CAATTACCCCCCCCGAACCATGCCGCGTTGGTCGCCGGGAACAAAATAATATTGGAAGGGAATATTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTa < 2:371090/137‑1 (MQ=255) cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTCCCGCCCCCGAACCATGACGCGATGGTCGCGGGAAACAAAATAATACGGGAAGCGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTa < 1:64191/137‑1 (MQ=255) cAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTaaa > 1:106865/1‑139 (MQ=255) cAAATACAGCAAAATTGTTTTCAGCCGGTTCCAC‑CCATTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATCCTGGAGGCGAATATTGCCGTGTTTCCCCCCGCCATATTCACTTTCAGCGGTa < 1:363019/137‑1 (MQ=255) aaaaTCCAGCAAAATTTTTGTCAGCCAGTTCCAC‑CCAGTCCCCCCCCCGAACCATGACGCGATGTCCGCCGGGAACAGAAAAATACTGAAAGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTa < 1:355043/136‑1 (MQ=255) aaaaTACAGAGAATTGGTTGTCAGCCAGTCCAAA‑CAATTCCCCCCCCCGACCCATGCCGGGGTGGTCGCCGGAAAAAAAATAATACGGGAAGGGAAGATTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTa < 2:951230/136‑1 (MQ=255) aaTACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtg > 1:572673/1‑139 (MQ=255) cAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg > 1:371090/1‑139 (MQ=255) gcaataTTGTTGCCAGCCAGTTCCAC‑CCACTCCCCCCCCCGAACAATGCGGCGATGGCCGCGGGGAAAAAAATAATTCTGGAAGCAAAGTTTGCCGGAATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 2:106865/134‑1 (MQ=255) gcaaaaTTGTTGTCAGCCAGTCCCAC‑CCAGTCCCGCCCCCGAACCATGACGCGAGGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:327465/136‑1 (MQ=255) gcaaaaTTGTTGTCAGCCAGTCCCAA‑CCAGTACCCCCCCCGAACCATGACGCGATGGTCGCCGGGAACAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:691538/136‑1 (MQ=255) gcaaaaTTGTTGCCAGCCAGTTCCAA‑CCAGTCCCCCCCCCAAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:858609/136‑1 (MQ=255) gCGAATTTGTTGTCACCCATTCCCAC‑CCAGTCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAGGCGAAGATTGCCGGGTTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:778260/139‑1 (MQ=255) gCGAAATTGTTGTCATCCAGTCCCAA‑CCAGTCCCCCCCCCAAACCATGACGCGTTGGTCGCCGGGAACAGAAAAATACGGGAAGGGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 2:251101/139‑1 (MQ=255) gCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgctgt > 1:903521/1‑135 (MQ=255) gCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct > 1:296078/1‑139 (MQ=255) gCGAAATTGTTGCCAGCCAGTTCCAA‑CCAGTCCCGCCCCCGAACCATGACGCGATGGGCGCCGGGAACAGAATAATCCTGGAAGCAAATAGTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:922161/139‑1 (MQ=255) cccaaaTTGTTTACACCCATTACCAA‑CCATTCCCGCCCCCAACCCATGACGCGATGGTCGCCGGAAACAGAAAAATCCTGAAAGCGAAGTTTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:540358/136‑1 (MQ=255) ttttttgcCAGCCATTCCCAA‑CCAGTACCGCCCCCAACCCATGCCGCGATGGCCGCGGGGAACAAAAAAATACTGAAAGCGAAGATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg < 1:741767/135‑1 (MQ=255) atggttgcCACCCAGTTCCAC‑CCAGACCCCCCCCCAACCCAGGCCGCGATGGTCGCGGGAAACAAAATAATACTGAAAGCGAAGATTGCGGTGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg < 2:600442/136‑1 (MQ=255) tgttgcAAGCCA‑TTCCAACCCAGTACCGCCCCCGAACCAGGACGCGATGGCCGCCGGGAAAAAAATAATACGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc < 1:656712/138‑1 (MQ=255) aGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa > 1:83780/1‑139 (MQ=255) cacccAGTCCCGCCCCCGAACCAGGCCGCGATGTTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACAcgacga < 2:269371/135‑1 (MQ=255) gacccATGACGCGATGGCCGCCGGGAACAAAATAATCCGGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg < 2:551540/136‑1 (MQ=255) gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg < 1:219826/139‑1 (MQ=255) gAACCATGACGCGATGGTCGCCGGAAACAGAATAATACTGGAAGCAAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg < 1:914106/139‑1 (MQ=255) tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt > 1:825882/1‑138 (MQ=255) tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt > 2:101812/1‑138 (MQ=255) tGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt > 2:483037/1‑138 (MQ=255) gACGCGAGGGCCGCCGGGAACAAAATAATCCTGAAAGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa < 1:314006/139‑1 (MQ=255) gcgATGGTCGCGGGGAACAAAATAATCCTGAAAGCAAAGTTGGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACca < 2:408478/139‑1 (MQ=255) gcgcGCCGGAAACAAAATAATACGGAAAGCGAAGAGTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg < 2:52224/136‑1 (MQ=255) aataatACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc > 1:912856/1‑139 (MQ=255) taCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaa > 2:575707/1‑138 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct > 2:33922/1‑59 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct < 1:33922/59‑1 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa > 2:767902/1‑139 (MQ=255) gAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGAAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGAGTAGTTTACCACAATGCGGCGTTGACCACGCTAAACAAAt > 2:194262/1‑138 (MQ=255) aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg < 1:167360/96‑1 (MQ=255) aaGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg > 2:167360/1‑96 (MQ=255) aGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac > 2:796176/1‑138 (MQ=255) gAAGATTGCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca < 2:104271/139‑1 (MQ=255) ttGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa < 1:576940/139‑1 (MQ=255) ggggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc < 1:50061/138‑1 (MQ=255) cGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc < 1:418162/139‑1 (MQ=255) aTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg < 2:296078/139‑1 (MQ=255) aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCa < 1:575707/138‑1 (MQ=255) | ACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCA > NZ_CP009273/3436502‑3436773 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |