| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,532,988 | 0 | G | A | 40.5% | 53.7 / 18.9 | 43 | G194S (GGC→AGC) | cysK | cysteine synthase A, O‑acetylserine sulfhydrolase A subunit |
| Reads supporting (aligned to +/- strand): ref base G (12/13); new base A (17/0); total (30/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.55e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGC > NC_000913/2532842‑2533135 | gCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCaa < 2:388973/151‑1 (MQ=255) cAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAg < 2:212234/151‑1 (MQ=255) aaTCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGTCTGGCGGGAGCCGCTC‑C‑ATTTAAAGGGCCAAAAGCa > 2:474324/1‑148 (MQ=255) tGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTAAGCTTACTTGCGTCAGCCGCGC‑C‑ATTAAAGGAACCAAAAGCAAGAACGAttttttttc > 1:281955/1‑145 (MQ=255) gAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACTGGCGGCAGCCGCCA‑C‑CTTAAAAGGACCAAAAGAAAAACCCAACtttt > 1:316042/1‑144 (MQ=255) aTTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACCCTGACTGGCGGCAGCCGCCA‑C‑ATTAAAAGCACCCAAAGGAAAACCGAACtttt > 1:77784/1‑141 (MQ=255) aTTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACCCCGACTGGGGGCCGCCGCCA‑CAATTAAAGCACCCAAAGGCAAAACGGacttttttctttt > 1:714596/1‑138 (MQ=255) gAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGGCTGGCGGCCGCCGCCACC‑ATTAAAGAC‑CCAAAAGGAAGACCGAACTTATTTCTGTCGCCggt > 2:298067/1‑148 (MQ=255) gAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGCCTGGCGCCAGCCCCTC‑C‑ATTAAAAGCACCAAAGGGAAAAACGAGCTTTTCTCTTTTGcc > 2:108314/1‑147 (MQ=255) aaaaGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGGGGCAGCCGCCA‑CAATTAAAGG‑ACCAAAAGCAAAAACGATCTTTTCTCTTTCGcc > 1:756094/1‑146 (MQ=255) aaGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTTGGCGGAATGGTGGCAGTCGCTC‑C‑ATTAAAAGGACCCACTGGAAGAACGAGATTTTCTTTGTctgcgt > 2:117841/1‑141 (MQ=255) accaccGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCCA‑C‑ATTAAAAGCACCAAAAGCAAAAACCAACTTTTTTCTTTCGCCGTTGAggc > 1:574615/1‑147 (MQ=255) ccaccGGTCCGGGGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCa < 1:559338/149‑1 (MQ=255) acaccGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACCGGCGTGAGCCGCCA‑C‑CTTTAAGGGACCAAAAGGAAGACCGAACTTTTCTCTGTCGcccgttg > 2:256003/2‑140 (MQ=255) caccGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAAcc < 1:566130/151‑1 (MQ=255) gTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACACAGACTGGCGTCCGACGCTA‑C‑ATTTAAAGCACCAAAAGCAAGACCGATCTTTTCTTTTTCGCCCTTTAGCCAAcc > 2:799059/1‑146 (MQ=255) ccGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGAAAGCCGGCTGGCGTCTGCCGCCA‑C‑ACTTAAGAGACCAAAAGCGAGAACGATCTTTTCTACTTTGCCGCTGGGGCAACCGAta > 2:474834/1‑147 (MQ=255) ggAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTTACTGGCGGCAGCCGCCA‑C‑ATTTAAAGCACCCAAAGCAAGAACCGCCTTTTCTCTTTTGCCCGTGAGCCCACCCGTTCTcc > 1:56835/1‑150 (MQ=255) ggAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTTACTGGCGGCAGCCCCTC‑C‑ATTTAAAGCACCACAAGCAAAACCCAGCTTTTTTCTTTCGCCCTTGAGCCAAcc > 2:196352/1‑142 (MQ=255) gagaTATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGCACGCTTAATGGGGTGAGGCGGCA‑A‑ATTTAAAGGACCAAAAACAAGACCCATCTTATCCCTTTCGCCGGTGAGCCCAccccttc > 1:401172/1‑141 (MQ=255) gATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGGCTGGCGGGCGCCGCTC‑C‑ATTAAAAGCAACCAAAGCAAAAACGATCTTTTTTCTTTTGCCGCTGAGCCCACCGATTTTCCCGTTATTGccc > 1:16924/1‑148 (MQ=255) aCCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGCACCCCGATTGGCGTCTGCCCCTA‑C‑ATTAAAAACACCAAAAGGAAAAACGCTCTTATTTCTTTCGCCCGTGGGGCAAACGAATCTCCTCTTATCccccc > 2:728196/1‑145 (MQ=255) cAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTg < 1:815738/151‑1 (MQ=255) ggTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTATGCTGACCTGCGTGTTCCGCTT‑C‑AATAAAAGAACCAAAGGGAAGAGCGCACTTGTCTTTGTTGCCCTTGTGCCCAACCCTTCTCCAATTTTAGCCCATGGCCCTGGC‑GGTGa > 1:828937/1‑150 (MQ=255) ggTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGACACCGGCTGGGGGTAGCCGCCA‑C‑AATAAAAGAACCCAAAGGAAGAACGATCTTATCTTTTTCCCCGTTTTGCCCACCCATTTTCCAGACATCGCCCC‑GGCGcgggcggggg > 2:3958/1‑147 (MQ=255) tGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGAt < 2:710866/150‑1 (MQ=255) gTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACTGGCGGCAGCCGCCA‑C‑ATTAAAGGCACCAAAAGCAAGACCGATCTTTTCTCTTTCGCCGTTGAGCCCACCGATTCTCCCGTTATTGCCCA‑GGCGCCGGCAGGGGGAGAGAtt > 1:696207/1‑150 (MQ=255) atttattGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTaa < 2:16059/150‑1 (MQ=255) tattGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAAcc < 1:245886/150‑1 (MQ=255) gACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAAtt < 1:7094/151‑1 (MQ=255) ggCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCa < 1:83359/149‑1 (MQ=255) gCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGgc > 1:626633/1‑151 (MQ=255) cTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTgg > 2:629952/1‑151 (MQ=255) ggCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGtt > 1:637193/1‑151 (MQ=255) ggCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGt < 1:323756/150‑1 (MQ=255) gTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGtttt > 2:736554/1‑150 (MQ=255) gCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGc < 1:415153/151‑1 (MQ=255) gCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGc > 1:151084/1‑151 (MQ=255) tA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCt > 2:630565/1‑150 (MQ=255) a‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTa > 1:240346/1‑150 (MQ=255) aCCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCa < 1:581308/151‑1 (MQ=255) cAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCa > 1:658539/1‑149 (MQ=255) aaaGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGc > 1:945264/1‑151 (MQ=255) | GCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTA‑C‑ATTAAAGGCACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCA‑GGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAAGC > NC_000913/2532842‑2533135 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |