Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,734,269 | 0 | T | G | 55.3% | 30.5 / 47.7 | 38 | V402G (GTG→GGG) | avtA | valine‑‑pyruvate transaminase |
Reads supporting (aligned to +/- strand): ref base T (6/11); new base G (21/0); total (27/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT > NZ_CP009273/3734186‑3734358 | aGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGgagaa > 1:3735/1‑89 (MQ=255) aGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGaa > 1:90694/1‑89 (MQ=255) gggCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGAtt > 2:244975/1‑90 (MQ=255) ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTggg > 1:295988/1‑89 (MQ=255) ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGc > 2:462162/1‑90 (MQ=255) ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGc > 1:65587/1‑90 (MQ=255) tAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGa < 2:267610/90‑1 (MQ=255) aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAag > 1:13414/1‑90 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATa > 2:37355/1‑90 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATa > 2:293447/1‑90 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATa > 2:94427/1‑90 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATa > 2:725377/1‑90 (MQ=255) gCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAag > 2:385696/1‑90 (MQ=255) cGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt < 2:146533/90‑1 (MQ=255) cGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt < 1:244975/90‑1 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTg > 2:686571/1‑90 (MQ=255) aaTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAACCTGGGCTg > 1:1588/1‑90 (MQ=255) aTGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGGGGGAGAGAAAGAAAGAGCCTGGGCTGAAAGt > 1:410853/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 2:186860/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 2:344801/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 2:274910/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 2:301399/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 2:323491/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGGCACTaa > 2:752364/1‑90 (MQ=255) aCCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa < 1:483662/90‑1 (MQ=255) aCCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa < 1:86896/90‑1 (MQ=255) aCCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa < 1:186860/90‑1 (MQ=255) ccAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGGTGAAAGTCACTAAAAAATGCCTGAt > 1:112571/1‑90 (MQ=255) gagCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGCTGAAAGGCACTAAAAATTGCCTGATgcg > 1:556200/1‑90 (MQ=255) gCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCt > 2:548948/1‑90 (MQ=255) cGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTAc < 1:344801/90‑1 (MQ=255) ggAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACg > 1:327288/1‑90 (MQ=255) ggAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACg > 2:587675/1‑90 (MQ=255) gagaAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGc < 1:546207/90‑1 (MQ=255) aaTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc < 1:329618/90‑1 (MQ=255) aTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCa < 1:548948/90‑1 (MQ=255) ggggggGAAGATTCTGGCGGGAGAGAAAGAAAGAGCCTGGGGTGAAAAACACTAAAAATTGCCTGATGCGCTACGCTTATCAAGCCTTCa > 2:662838/1‑90 (MQ=25) gggggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa > 1:815866/1‑66 (MQ=255) gggggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa < 2:815866/66‑1 (MQ=255) tGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAtt < 2:327288/90‑1 (MQ=255) tGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAtt < 1:385696/90‑1 (MQ=255) | AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT > NZ_CP009273/3734186‑3734358 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |