Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,595,044 | 0 | T | A | 48.6% | 26.5 / 19.2 | 35 | Y146F (TAC→TTC) | livM | branched‑chain amino acid ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base T (14/4); new base A (0/17); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCACC > NC_000913/3594900‑3595181 | aGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAAc > 1:181635/1‑151 (MQ=255) cgcgCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTAccc > 2:39172/1‑151 (MQ=255) agcaCCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGCGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCccagcaccagc > 2:321308/1‑150 (MQ=255) ccccgAAAACCGGGGGGGGAGCCCGCCCCCCGCGCCATTAAACCCACAAAC‑CGG‑CGGGAGGTCCAAGAGGCCCAGGCCGTAATAGTGATTGAGCAGCACAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACcagcag < 1:141182/146‑1 (MQ=255) aCCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGAc > 2:463678/1‑151 (MQ=255) gAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGAcc > 2:318385/1‑148 (MQ=255) aaaaCCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGAGATTGAGCAGCGCAAAAGTGTAAGCGCCGCTGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGa > 1:211034/1‑150 (MQ=255) gAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGAcc > 2:173775/1‑150 (MQ=255) aagaaacccggccgccgtTGCCA‑TAAACCCGAAATT‑GGCAAGGAAGGTCA‑AAAACCCCAATCCGAAAGAGTGATTGAGCAGCGAAAAGGTGTAGGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACcaac < 2:401515/140‑1 (MQ=255) gaacgCCCGCCGCCGCTCCCA‑TAAACCAACAAACGGGC‑GGGCGGGTCA‑AAAAGCCCAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACcaacaa < 2:293719/146‑1 (MQ=255) ccccgccgcTGCCATTAACCCAGAAACC‑GGC‑GGGAGGGTCA‑AAAGGCCCAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTc < 2:43260/149‑1 (MQ=255) gcTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCg > 1:198333/1‑150 (MQ=255) ttAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCg > 2:132910/1‑149 (MQ=255) aacc‑GGC‑AGGCAGGTCC‑AGACCCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGAt < 1:463559/141‑1 (MQ=255) ggc‑gggcgGGCCC‑AGAAGCCAAAGCCGTAATGGTGATTGAGCGGCGAAAAAGTGTAAGCGCCGATGGCGAAAAACCCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGt < 1:490545/145‑1 (MQ=255) ggccaggTCC‑AGACCCCCAACCCGAAATAGGGATGGCGCAGGGCAAAATTGTAACCCCGAAGGGAGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATgg < 1:173715/148‑1 (MQ=255) aggcaggTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATgg > 2:127182/1‑151 (MQ=255) gccgggTCA‑AGACGCCCAGGCGGTAATAGTGATGGAGCAGCGCAAAAGTGTAGGCGCGGATGGGATAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATg < 2:77653/145‑1 (MQ=255) caggTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGTCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTc < 1:360543/150‑1 (MQ=255) gggcaagaaaGCCAAAGCTGAAAAAGTGTTTGAGCAGCGCAAAAGTGAAAGCGCGGATGGAGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAg < 1:282618/143‑1 (MQ=255) ggccc‑AGACCCCCAAGCCGTAATGTTGATAGAGAAGCGCAAAGGGGTAAGCGCCGTGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTc < 1:313654/145‑1 (MQ=255) cccAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGAAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTgg < 2:54918/145‑1 (MQ=255) caaaGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATc < 2:228846/149‑1 (MQ=255) aaGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATcc > 2:351003/1‑150 (MQ=255) ccGTAAGAGAGATTGGGCGGCGCAAAAGTGGAAGCGCCGGGGGCAAAAAAACCTCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCg < 2:179596/151‑1 (MQ=255) gTGATGGAGCAGCGAAAAAGGGAAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGt < 1:268915/151‑1 (MQ=255) ttGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGaa > 1:181072/1‑150 (MQ=255) ttGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTg < 2:181695/148‑1 (MQ=255) accatcgaaaaaGTAAAACCGCCGATGGCATAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGaa < 1:48336/139‑1 (MQ=255) gcagcGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAAcc < 2:181132/149‑1 (MQ=255) cgcAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAAc > 1:399762/1‑151 (MQ=255) aaaaGTGTAGGCGCGGGGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACgg < 1:319298/150‑1 (MQ=255) aaGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGgcc > 2:185093/1‑150 (MQ=255) tggaaGCGCCGATGGCGAAAAACCCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGgccac < 1:533604/146‑1 (MQ=255) ccGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCAcc < 2:474054/147‑1 (MQ=255) | AGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATC‑GGC‑AGGCAGGTCC‑AGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCACC > NC_000913/3594900‑3595181 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |