Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,561,000 | 0 | A | C | 55.2% | 27.3 / 38.2 | 30 | V131G (GTG→GGG) | yjiR | PLP‑dependent aminotransferase family protein |
Reads supporting (aligned to +/- strand): ref base A (11/2); new base C (0/16); total (11/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTA > NZ_CP009273/4560914‑4561089 | gCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCGccc < 1:899637/90‑1 (MQ=255) gCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCAccc < 2:260093/90‑1 (MQ=255) aTCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCgg > 1:233639/1‑90 (MQ=255) gggCATTCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCt > 1:1062612/1‑90 (MQ=255) gggCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCt > 2:398891/1‑90 (MQ=255) aTCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGc > 1:227674/1‑90 (MQ=255) tctcGCAATCCTCGCTGCCCGGCTAACAAGTCATACCCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGc < 2:904174/90‑1 (MQ=255) tctcGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCACCCGGCTTAGCTCCCGCCACAGc < 1:421523/90‑1 (MQ=255) tctcGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGc > 1:818520/1‑90 (MQ=255) cAATCCTCGCTGACCGGCAAACAAGCCAAAACCGAAAAGGGTTTGCAGTTTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGttt < 1:815313/90‑1 (MQ=255) tACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAg > 1:989015/1‑90 (MQ=255) cccGGCTAACAAGTCAAACCCGAGAACGGTTTGCAAATTATGCTGCCCCCCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGgcgc < 2:1062612/89‑1 (MQ=255) aCCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGgcgc > 2:933969/1‑90 (MQ=255) cGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCtt < 2:333505/90‑1 (MQ=255) cGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCtt < 2:224803/90‑1 (MQ=255) aCAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACAt < 1:478372/90‑1 (MQ=255) tCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGgt < 1:326744/90‑1 (MQ=255) cccGAGAACGTTTTGCAGATTATGCTGCCCCCCCCGGTTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCt < 2:703922/89‑1 (MQ=255) gagaACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCtttt > 2:519722/1‑90 (MQ=255) aCGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCtt < 2:370353/90‑1 (MQ=255) aCGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCtt < 2:370352/90‑1 (MQ=255) ttttGCAGATTATGCTGCCCCCCCCGGTTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAAc < 2:744918/89‑1 (MQ=255) tGCAGTTTATGCTGCACCCCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGa < 1:547687/90‑1 (MQ=255) tGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGa < 1:869788/90‑1 (MQ=255) tGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGa < 2:351483/90‑1 (MQ=255) cAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAAc < 1:999172/90‑1 (MQ=255) aTTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAt > 2:914042/1‑90 (MQ=255) tgctgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGaa > 1:319975/1‑90 (MQ=255) gctgcCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAAc < 2:256987/90‑1 (MQ=255) tgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTg > 2:144639/1‑90 (MQ=255) caccacCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCtt < 2:865876/47‑1 (MQ=255) caccacCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCtt > 1:865876/1‑47 (MQ=255) accCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTa > 2:447845/1‑90 (MQ=255) accCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTa > 1:765984/1‑90 (MQ=255) accCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGCGCTTTTGCTTACCGGAACACTGGAACTGTCGCTa > 2:654735/1‑90 (MQ=255) | GCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTA > NZ_CP009273/4560914‑4561089 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |