| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,740,909 | 0 | T | G | 37.8% | 48.2 / 13.3 | 37 | V402G (GTG→GGG) | avtA | valine‑‑pyruvate aminotransferase |
| Reads supporting (aligned to +/- strand): ref base T (8/15); new base G (14/0); total (22/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.00e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTT > NC_000913/3740778‑3741046 | ccTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGAtt < 2:49868/139‑1 (MQ=255)ccTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGAtt > 1:49868/1‑139 (MQ=255) cACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGaa < 2:843255/140‑1 (MQ=255) gcggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAagag > 2:71347/1‑140 (MQ=255) gcggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAagag > 1:949269/1‑140 (MQ=255) cggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAgaga < 1:148834/140‑1 (MQ=255) cgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAg > 1:175385/1‑140 (MQ=255) gTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTg < 2:949269/140‑1 (MQ=255) ggCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAAAAAGAGCCTGGGCTGaa > 2:807004/1‑140 (MQ=255) ggCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGaa > 1:758804/1‑140 (MQ=255) cAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTg > 2:648451/1‑140 (MQ=255) ggCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCt > 1:877626/1‑140 (MQ=255) gCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTg < 2:265736/140‑1 (MQ=255) aTAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcg < 2:47908/140‑1 (MQ=255) ccGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACg < 2:758804/140‑1 (MQ=255) ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTa < 1:648451/140‑1 (MQ=255) ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTa > 1:91157/1‑140 (MQ=255) ccGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc < 1:377869/140‑1 (MQ=255) cATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAgg < 2:877626/140‑1 (MQ=255) tGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAttt > 1:656409/1‑140 (MQ=255) gcgcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCa < 1:718316/140‑1 (MQ=255) gcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAat > 2:301012/1‑140 (MQ=255) tGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAatata < 1:71347/140‑1 (MQ=255) tGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAatata > 1:377495/1‑140 (MQ=255) tGAACTACGTACCAGAGCCCGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAatatc > 1:104270/1‑139 (MQ=255) cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTGTCAGGCCTACATGATTTCTGCAATATATTGAAtt > 1:540426/1‑140 (MQ=255) tACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTg > 1:212917/1‑140 (MQ=255) cAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGgaa > 1:229130/1‑140 (MQ=255) gagCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGgaaga > 1:319807/1‑140 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGgaagaa > 1:944356/1‑140 (MQ=255) gagaAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGt > 1:287899/1‑140 (MQ=255) gAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGAt < 1:807004/140‑1 (MQ=255) aGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATa > 1:208047/1‑140 (MQ=255) ggCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATaa < 1:839620/140‑1 (MQ=255) ggCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATaa > 1:947946/1‑140 (MQ=255) gCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAg > 2:871398/1‑140 (MQ=255) ggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGtt < 2:319807/140‑1 (MQ=255) | CCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATAAGGCGTT > NC_000913/3740778‑3741046 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |