Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 669,219 | 0 | G | C | 56.2% | 4.3 / 15.8 | 16 | H341D (CAC→GAC) | leuS | leucine‑‑tRNA ligase |
Reads supporting (aligned to +/- strand): ref base G (4/3); new base C (0/9); total (4/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.92e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGAGCCGTCAGCTGCCAGGATAACCGGTTT‑GATGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAAT > NZ_CP009273/669140‑669304 | tCAGACCCGTCAGCGGCCAGGATAACCGGTTT‑GATGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTAcc < 2:292243/90‑1 (MQ=255) agagCCGTCAGCTGCCAGGATAACCGGTTT‑GATGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGc > 1:527947/1‑90 (MQ=255) gATACCGGGTTG‑GTTGTTCGGCCCTTTTTAAGGGGAAACCTGGTGGTCGCGCGGGTCGGCCCCCGGTACCGCCATAACTGCGCCCGTGcc < 1:744380/90‑1 (MQ=255) gataCCGGTTTGGATGTTCGGGCCTTTTTTAGGGGCAACCTCGTGGCCGCGCTGGTCGCCCCCGGGTACCGCCATAACTGCGCCCGTGcc < 1:472175/87‑1 (MQ=255) tt‑GATGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATc < 2:527947/90‑1 (MQ=255) aTGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAAt > 1:337523/1‑90 (MQ=255) tGTTCAGGCCTTTTTTAGAGGCAAACCCGTAGCCGCGCTGGTCGTCCCCCGGTCCCGCCATAACTGCGCCCGTGCCGTACTCCATCAATa < 1:214379/90‑1 (MQ=255) cAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGaa > 2:722020/1‑90 (MQ=255) aGTCGCCCTGGTCGTCCCCGGGTCCCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATttctt < 1:89912/90‑1 (MQ=255) gccgcgCTGGTCGCCCCCGGTTCCCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATttcttc < 1:722020/88‑1 (MQ=255) gccgcgCTGGTCGCCCCCCGTTCCCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATttcttc < 1:390333/88‑1 (MQ=255) cgcgCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTAc > 1:716280/1‑47 (MQ=255) cgcgCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTAc < 2:716280/47‑1 (MQ=255) gcgcGGGTCGCCCCCGGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGcc < 2:337523/90‑1 (MQ=255) gccttGGTCGTCCCCGGTTCCCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGcc < 1:546413/86‑1 (MQ=255) tCGTCCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAAt < 2:82661/90‑1 (MQ=255) | TCAGAGCCGTCAGCTGCCAGGATAACCGGTTT‑GATGTTCAGGCCGTATTTAGAGGCAAACTCGTAGTCGCGCTGGTCGTGCCCCGGTACCGCCATAACTGCGCCCGTGCCGTACTCCATCAATACGAAGTTTGCTGCCCAAACGGGAATTTCTTCGCCCGTTAAT > NZ_CP009273/669140‑669304 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |