Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 510,080 | 0 | A | C | 55.0% | 19.9 / 28.3 | 21 | V79G (GTG→GGG) | ybbJ | NfeD family protein |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/11); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGC > NZ_CP009273/509996‑510166 | aGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgc > 1:474483/1‑90 (MQ=255) aGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgc > 2:222746/1‑90 (MQ=255) cGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgccg > 2:742757/1‑90 (MQ=255) gaTCCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCCCCCGCCGCGAc < 2:49061/90‑1 (MQ=255) cccaccaCAAACCGTCGCCAAACCGGCTCCTCCCCGCGTTGTTTTAAAGGACTATCGCTGTCTTTTTGTTCCCGCGCCCGCCGCGACAAc < 2:368954/89‑1 (MQ=255) cccAGCAAAAACCGTCGGCCAATCAGCTGCCCCCCGCGCTGTTTTAAATGACTATCGCTGTGTTTTTGTTCGCGCCCCCGCCGCGACAAc < 2:340608/89‑1 (MQ=255) cccAGCAAAAACCGCCGCCCAATCCGCTGCCCCCCGCGCTGTTTTAAATCACTACCGCTGTGTTTTTTTTCGCCCCCCCGCCGCGACAAc < 1:570711/89‑1 (MQ=255) cAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCa > 2:674345/1‑90 (MQ=255) cAGCAAAAAACGTCGGCCAATCAGCTGCTCCCCGCGCTGTTTTAAATGACTATCGTTGTTTTTTTGTTCGCCCCCCCGCCGCGACAACCa < 2:1071854/90‑1 (MQ=255) cccAATCGCCTGCGCCCCGCGTTGTTTAAAATGACTACCGCTGTTTTTTTGTCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCa < 1:222746/89‑1 (MQ=255) tCAGCTGCGCCCCGCGTTGTTTAAAATGCCTTCCGCTGTTTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAggcgg < 1:151008/90‑1 (MQ=255) gctccCCGCGCTGTTTAAAATGACTTCCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca < 1:108583/90‑1 (MQ=255) gccccccGCGTTGTTTAAAATGCCTTTCGTTGTTTTTTTTTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca < 1:1017325/86‑1 (MQ=255) gcTGTTTAAAGTGACTCCCGCTGTTTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGa < 1:775086/90‑1 (MQ=255) aCTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAc > 2:1051459/1‑90 (MQ=255) aCTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAc > 1:1039290/1‑90 (MQ=255) cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAcc > 2:199165/1‑90 (MQ=255) ttgtttttGTCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCa < 1:742757/86‑1 (MQ=255) gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCa > 1:931228/1‑90 (MQ=255) tttGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc > 1:29978/1‑90 (MQ=255) cgcCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGc < 1:199165/90‑1 (MQ=255) | AGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGC > NZ_CP009273/509996‑510166 |
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: ‑ |