Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,432,473 | 0 | A | G | 47.2% | 26.4 / 14.6 | 36 | S857P (TCT→CCT) | nifJ | pyruvate:ferredoxin (flavodoxin) oxidoreductase |
Reads supporting (aligned to +/- strand): ref base A (16/3); new base G (0/17); total (16/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAA > NZ_CP009273/1432388‑1432559 | ccAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGc < 2:1250201/90‑1 (MQ=255) aatttCGGCATTATCTTAAAAAAGAGAGTTCGCCCAGGCCGCCCCCCGCCCGTTGGCATGGGTGGTATACGGTGTAGGGGGCAGGTTAcc < 1:768603/87‑1 (MQ=255) ggCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATa > 2:424777/1‑90 (MQ=255) aTCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTga > 2:170932/1‑90 (MQ=255) cAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGGGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAAc > 2:860819/1‑90 (MQ=255) aaaTAGAGATTTCCCCCATGCCGGCCCACGACCGTTGGCATGGGTGGTATCCGGTGTAGGGGGCAGGTTACCGCCATAAATTGAAGAACa < 2:1260000/90‑1 (MQ=255) gagagTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAg > 2:523662/1‑90 (MQ=255) agagTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGt > 2:1111325/1‑90 (MQ=255) gTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGc > 2:972075/1‑90 (MQ=255) gTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGc > 1:1259999/1‑90 (MQ=255) cggggccgccCCCCGCCCTTGGGCTTGGGTGGTATAGGGTGTGGGGGGCGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGc < 1:113306/86‑1 (MQ=255) cggggccgccCCACGCCCGTGGGAAGGGGTGTTATCGGGTGAAGGGGGGAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGc < 1:1470789/86‑1 (MQ=255) tGCCGCCCCCCGCCCTTTGGCATGGGGGGTATACGGTGAAGGGGGAGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGAt < 2:578497/90‑1 (MQ=255) tGCCGCCCCACGCCCGTGGGCATGGGTGGTATAGGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGAt < 2:709904/90‑1 (MQ=255) gccgccCCACGCCGGTGGGCATGGGGGTTATGGGGTGGGGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATc < 1:795252/90‑1 (MQ=255) ccACGACCGTGGGCATGGGTGGTATAGGGTGTAGGGGGCGGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATc < 2:1643652/90‑1 (MQ=255) cGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATcctg > 2:1470790/1‑88 (MQ=255) ccggTTGGAATGGGTGGTATGGGGTGTAGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc < 1:645829/87‑1 (MQ=255) cccGTGGGCATGGGGGTTATCGGGTGAGGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc < 1:1111324/89‑1 (MQ=255) aCCGTGGGCAGGGGTGATATGGGGTGAAGGGGGGAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTc < 1:424777/90‑1 (MQ=255) gTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGcc > 2:1047998/1‑90 (MQ=255) tGGCATCGGTGGTATCCGGTGTAGAGGGAAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCAt < 1:366678/90‑1 (MQ=255) tGGCACCGGTGTTACACGTTTTAGGGGGCAGGTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCAt < 2:119252/90‑1 (MQ=255) gggtggtATCGGGTGTAGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCt < 2:1428177/89‑1 (MQ=255) ggtggtATCGGGTGAGGGGGGAAGTTCCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTg < 1:357209/90‑1 (MQ=255) gtggtATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGa > 1:1205700/1‑90 (MQ=255) gtATGCGGTGAAGGGGGAGGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTc < 2:1272539/90‑1 (MQ=255) tggggTGTAGGGGGGGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGt < 2:827064/87‑1 (MQ=255) gtgtAGGGGGCGGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt < 2:1205701/90‑1 (MQ=255) gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACGGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt > 1:1846685/1‑90 (MQ=255) gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCTTAGAGCTGAGTCAGTAAtt > 1:210951/1‑90 (MQ=255) gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTaact > 1:1428176/1‑88 (MQ=255) gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt > 2:1413703/1‑90 (MQ=255) gtgtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAtt > 2:1449913/1‑90 (MQ=255) ggtAGGGGGCAGTTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAAttt < 1:523662/89‑1 (MQ=255) gtAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTa > 1:619671/1‑90 (MQ=255) gagagGGCAGGTTACCGCCATAAATTAACGAACAGCAAGTGGGGGTAGAGATGAACATCCGGCCGCCATAGAGATGAGCCAGTAATTTaa > 1:170060/2‑90 (MQ=255) | CCAAGGCCAAATTCGGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAA > NZ_CP009273/1432388‑1432559 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |