Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,712,188 | T→C | D26G (GAT→GGT) | pdxH ← | pyridoxamine 5'‑phosphate oxidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,712,188 | 0 | T | C | 90.0% | 64.5 / ‑5.1 | 31 | D26G (GAT→GGT) | pdxH | pyridoxamine 5'‑phosphate oxidase |
Reads supporting (aligned to +/- strand): ref base T (1/0); major base C (15/12); minor base A (3/0); total (19/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01 |
GATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCG > NZ_CP009273/1712061‑1712305 | gATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGt < 2:626155/139‑1 (MQ=255) tGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTaa > 1:93585/1‑139 (MQ=255) cgcTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAAcc < 1:861876/139‑1 (MQ=255) cgcTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGAACGCGGCG‑GCGTaaac > 1:745641/1‑137 (MQ=255) aTAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCt < 1:707370/139‑1 (MQ=255) gACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGATAATGGATCGGCGGAAGGATCGCGGCGAGGGCAACCGCCCTTTggcag > 1:488729/1‑136 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGACCGCGGGG‑GCGGAACCCGCCTTTGGTGTAtt > 1:565272/1‑139 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGACCGCGGCG‑GCGGAACCCCCCCTTTGGGTAtt > 2:510199/1‑139 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACcgcg < 1:538838/113‑1 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACcgcg > 2:538838/1‑113 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTAtt < 2:684099/139‑1 (MQ=255) ccATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGGCACACCCCCTTGGTGGAtt > 2:565811/1‑139 (MQ=255) agacagaTCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCa < 1:540579/133‑1 (MQ=255) aGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGACCGCGGCG‑GCGGAAACCGCCCTTTGTGTATTTTCGGCGCAGATGcgcg > 2:820779/1‑139 (MQ=255) cGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GGGTAAACCCCCTTTTGTGTATTCACGGCGCAGATGcgcggg > 1:250265/1‑137 (MQ=255) cGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAAACCGCCTTTTGTGGATTCACGGCGCAGATGCGCGAt > 1:424262/1‑139 (MQ=255) cGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGGAAACCGCCCTTTGGGTATTTACGGCGCAGATgctgtct > 2:540579/1‑134 (MQ=255) accaccATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCttt < 1:181665/110‑1 (MQ=255) accaccATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCttt > 2:181665/1‑110 (MQ=255) accaTCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGACCGCGGCG‑GCGGAAACCGCCTTTGGGGGATTTACGGCGCAGCTGCGCGATTtgct > 1:208594/1‑139 (MQ=255) gggTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGAACGCGGCG‑GCGTAACCCCCCTTTTGTGTATTTACGGCGCCGAAGGGCGGTTTGCTGCAAATcgtcgt > 2:220520/1‑139 (MQ=255) gCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGGAGACCGCGGCG‑GCGTAACCCCCCCTTTGGGTATTTCCGGCGCAGATGCCCGGTTTGCTGCAATTTGTCGGTAt > 2:266515/1‑136 (MQ=255) acaAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCtt < 2:754220/138‑1 (MQ=255) cttgggagagCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCATTATCAGACATGGTTTTCTTTACGGa < 2:488729/134‑1 (MQ=255) ccTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGa < 2:622034/139‑1 (MQ=255) ccTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGa < 2:93585/139‑1 (MQ=255) agagCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCa < 1:211547/139‑1 (MQ=255) ccAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCGCCTTTTGTGTTTTCACGGCGCCGCTGCGCGATTTTCTGCAATTCTTCGTTTTCCGACCTGGTTTTCTTTTCGGATTTTCAGTgg > 1:427389/1‑139 (MQ=255) cAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGACCGCGGCG‑GCGTAACCCCCCCTTTGGGTATTTACGGCGCCGAAGCGCGGTTTGCTGCAAATCTTCGTTATTAAACATGGTTTTTTTTTTGGGATGTCAGTggg > 2:71429/1‑139 (MQ=255) cGTTCAAAAAGGGTTAATGGATCGGCGGGGAGAAAGCGGCG‑GCGGAAACCCCCTTTTGTGGATTCACCGCGCAGATGCGCGGTTTTGTGCAAATCGTCGTTTATAGACATGGGTTTTTTTATGGGATGTTAGTgggg > 2:827874/1‑136 (MQ=255) gTTAATGGATCGGCGGGAAGACCGCGGGG‑GGGTAACCCCCCCTTGGTGTATTTACGGCGCAGATGCGCGATTTGCTGCAAATCGTCCTTTTTAGACATGGTTTTTTTTACGGGTATTCAGTGGGGGACGCCATtgtggg > 1:411992/1‑137 (MQ=255) | GATGCGCTGATAAGGCTGACCATGTTCATCCACGGTAGCGACCACCATCGCGGTAGGGTCCGCCAGTTTGGCTTCACAAGCCTGAGAGAGCCAACGTTCAAAAAGGGTTAATGGATCGGCGGGAAGATCGCGGCG‑GCGTAACCCGCCTTTGGTGTATTCACGGCGCAGATGCGCGATTTGCTGCAATTCGTCGTTATCAGACATGGTTTTCTTTACGGATTGTCAGTGGGTGACGCTATTGTGCG > NZ_CP009273/1712061‑1712305 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |