Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 687,303 | 0 | A | G | 50.0% | 43.8 / 13.5 | 46 | intergenic (‑63/+27) | corC/ybeY | CNNM family magnesium/cobalt transport protein CorC/rRNA maturation RNase YbeY |
Reads supporting (aligned to +/- strand): ref base A (5/18); new base G (23/0); total (28/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.39e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTT > NZ_CP009273/687219‑687388 | tgtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCt > 1:763359/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg > 1:650411/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg > 2:486747/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTg > 2:362883/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg > 1:451860/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg > 1:1133335/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTg > 2:173571/1‑90 (MQ=255) gtgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTg > 2:1173914/1‑90 (MQ=255) gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc < 1:688242/90‑1 (MQ=255) gtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGcc < 1:310008/90‑1 (MQ=255) aaTTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGGCTGcccgg > 1:482675/1‑87 (MQ=255) tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCCGCCAGTCaa > 2:562434/1‑90 (MQ=255) tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTTCCCGCCaa > 1:654362/1‑90 (MQ=255) tGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCaa < 2:1252682/90‑1 (MQ=255) tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTTCCCGTCaaaa > 1:214243/1‑89 (MQ=255) tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGAGGGCTGCCAGTCAAAt > 2:286252/1‑90 (MQ=255) tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGGTGCCCGTCcaaa > 2:654965/1‑89 (MQ=255) cgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCt < 2:1156451/90‑1 (MQ=255) cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg < 2:844916/90‑1 (MQ=255) aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGGCATTCAAAGCCCTGGcga > 1:364367/1‑88 (MQ=255) cgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCCGGCAAAAGCCCGGCAAagtattt > 1:170620/1‑89 (MQ=255) ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTtctc < 1:603584/90‑1 (MQ=255) tcttATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGGCTGCCAGTCAAATGCCTGGCAAACTATTCTTTCTCGGc < 1:1115543/90‑1 (MQ=255) tctcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGc < 1:1297868/90‑1 (MQ=255) tcATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCaa > 2:1252497/1‑90 (MQ=255) ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGGGTCTGCCAGTCAAAAAACTGGCAAATTATTTTTTTTCGGCAAAGTAc > 1:115311/1‑90 (MQ=255) ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCTGCAATGTAc > 1:1508530/1‑90 (MQ=255) ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTGCCAGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCAATGTAc > 2:1297868/1‑90 (MQ=255) ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc < 2:823162/90‑1 (MQ=255) ttagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc < 2:1424023/90‑1 (MQ=255) gttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTCTTTCTCGGCGATGTACggg > 1:1315137/1‑89 (MQ=255) taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATc < 1:963313/90‑1 (MQ=255) aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAAATATTTTTTTTCGGCAAAGTACGGATcc > 2:1457470/1‑90 (MQ=255) aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCAGTCAAATGCCCGGCAAAATATTTTTTTTCGGCAAAGTACGGGTcc > 2:1581332/1‑90 (MQ=255) aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGGGCGGCTTCCCGACAAAAACCCGGCAAATTATTTTTTTTCGGCCATGTACTGAttc > 1:1383955/1‑88 (MQ=255) aaaaCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCAt < 1:1252497/90‑1 (MQ=255) aaCGCCGTCGTTAATCACCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTTTCGGCGAAGTACGGATCCCCCTAg > 1:1032546/1‑90 (MQ=255) cGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGcc < 2:763359/90‑1 (MQ=255) gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGccc < 1:212668/90‑1 (MQ=255) ccGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCa < 2:1383588/90‑1 (MQ=255) gtcgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCaga < 2:1148522/90‑1 (MQ=255) cgtTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagc < 1:916231/90‑1 (MQ=255) cACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa < 1:223726/90‑1 (MQ=255) aCCAACGGCGGGGGCGTCTGCCCGTCAAAAGCCCGGCCAATTATTTTTTTTCCGCAATGTACGGGTCCTCATAGCCCAGAGCAAACATaa > 1:1091586/1‑90 (MQ=255) aaCGGCGGGGGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTTTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct > 1:1424357/1‑90 (MQ=255) ggggGCGGCTGCCCGTCAAAAGCCTGGCAAATTATTCTTTTTCGGCAAAGTACGGGTCCTCATAACCCCGAGCAAGCATAATCTCTGttt > 1:1532825/1‑90 (MQ=255) | TGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTT > NZ_CP009273/687219‑687388 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |