Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,734,269 | 0 | T | G | 56.2% | 20.4 / 43.7 | 32 | V402G (GTG→GGG) | avtA | valine‑‑pyruvate transaminase |
Reads supporting (aligned to +/- strand): ref base T (1/13); new base G (18/0); total (19/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.03e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACA > NZ_CP009273/3734186‑3734353 | aGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGAt > 2:699198/1‑90 (MQ=255) gggCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGAtt > 1:55508/1‑90 (MQ=255) ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGc > 1:154996/1‑90 (MQ=255) ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGc > 2:491662/1‑90 (MQ=255) tAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGa < 2:253447/90‑1 (MQ=255) aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag < 2:492320/90‑1 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGCGGgagagaga > 1:475494/1‑88 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGAGTCTTGTGGGAGAGATa > 2:118008/1‑90 (MQ=255) cATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGGCTgg > 2:398771/1‑90 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt < 1:280443/84‑1 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAGAAAAAGCCTGGGCTg > 1:478812/1‑90 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt > 2:280443/1‑84 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGGTg > 2:147191/1‑90 (MQ=255) aTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGa < 2:303535/90‑1 (MQ=255) tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAg < 2:535852/90‑1 (MQ=255) tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTTGCGGAAGAGATAGAAAAAGACTGGGCTGaaaa > 1:448939/1‑89 (MQ=255) tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAAATTCTTGCGGAAAAGAAAGAAAAAACCCGGGGTGAAAg > 2:549708/1‑90 (MQ=255) gcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTTTTGCGGGGGAGATAGAAAGAGCCTGGGCTGAAAGTCACt > 1:455115/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 1:180120/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTTTGGCGGGAGAGATAGAAAAAGCCTGGGGTGAAAGTCACTaa > 1:640723/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAGAAAAAACCTGGGCTGAAAGGCACGaa > 1:475348/1‑90 (MQ=255) gTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGCGGGAGAGATAGAAAGGGCCTGGGGGGAAAGTCAGTAAAAAATTCCt > 2:410260/1‑90 (MQ=255) gagCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGCGGAAGAGATAGAAAAAGCCCGGGCTGAAAGGCACTAAAAAATGCCTGATgcg > 1:639735/1‑90 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc < 2:317776/90‑1 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc < 1:548202/90‑1 (MQ=255) gagaAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGc < 1:549708/90‑1 (MQ=255) gagaAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGc < 2:639735/90‑1 (MQ=255) gagaAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGc < 2:85107/90‑1 (MQ=255) gagaAAATTGAGGCGGGGGGGAAGATTCTGGGGGAAGAGATAGAAAGAGCCTGGGCTGAAAGGCACTAAAAATTGCCTGATGCGCTACGc > 1:4626/1‑90 (MQ=255) aaTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc < 2:180120/90‑1 (MQ=255) aTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCa < 2:263756/90‑1 (MQ=255) gggggTGAAGATTCTGGCTGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa < 1:223186/90‑1 (MQ=255) | AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACA > NZ_CP009273/3734186‑3734353 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |