Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,177,828 | 1 | . | A | 63.0% | 67.4 / 31.1 | 27 | intergenic (‑3/‑371) | ribB/ubiK | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/ubiquinone biosynthesis accessory factor UbiK |
Reads supporting (aligned to +/- strand): ref base . (5/5); new base A (4/13); total (9/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTC > NZ_CP009273/3177745‑3177913 | gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 2:154/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 1:250970/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 1:270255/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 2:358795/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 1:484118/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 1:520418/1‑90 (MQ=255) gcgTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaa > 2:30584/1‑90 (MQ=255) ttCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCa > 2:222090/1‑90 (MQ=255) ttCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCa > 1:575804/1‑90 (MQ=255) aCGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACt < 2:205030/90‑1 (MQ=255) gcagcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTaaa > 2:2627/1‑90 (MQ=255) agcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAAt > 2:404205/1‑90 (MQ=255) gcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAAt > 2:79210/1‑90 (MQ=255) gcgcagcCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAAtt > 2:274095/1‑90 (MQ=255) ccAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGt < 1:418259/90‑1 (MQ=255) aGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTa > 2:450130/1‑90 (MQ=255) tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaaaacc < 1:373772/70‑4 (MQ=255) tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑aaaaaaaacc > 2:373772/1‑67 (MQ=255) tGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCa < 2:575804/90‑1 (MQ=255) gCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCa < 2:114858/90‑1 (MQ=255) acacGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAggg < 2:142439/90‑1 (MQ=255) acacGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAggg < 1:228902/90‑1 (MQ=255) acacGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGc < 1:129005/90‑1 (MQ=255) cacGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 1:561190/90‑1 (MQ=255) acGTTCGAAAGGCGTGCCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 2:100361/90‑1 (MQ=255) acGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTGTGGTTACCAGAATCAGGGCa < 2:548237/90‑1 (MQ=255) acGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 2:250970/90‑1 (MQ=255) acGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 2:280033/90‑1 (MQ=255) acGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 2:443469/90‑1 (MQ=255) gCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGt < 1:364078/90‑1 (MQ=255) gTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAggg > 2:171775/1‑75 (MQ=255) gTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAggg < 1:171775/75‑1 (MQ=255) aaGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa < 2:172134/70‑1 (MQ=255) aaGAGGAAAGTAGCGTCTGATTCATGGTAAAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCa > 1:172134/1‑70 (MQ=255) agGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGc < 1:154/90‑1 (MQ=255) aTGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTc < 2:484118/90‑1 (MQ=255) | GCGTCCTTCACGCAGCGCAGCCAGTGCATTTTCAACACGTTCGAAAGGCGTACCAAAAGAGGAAAGTAGCGTCTGATTCATGGT‑AAAAAAACCTCACTAAAATTATGGTTACCAGAATCAGGGCAGTCTTAGGAGTGGCGGCATATAGCCAAAATAACGTGAGCGGGTC > NZ_CP009273/3177745‑3177913 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |