| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,532,979 | 0 | G | A | 40.6% | 40.4 / 12.6 | 32 | G191S (GGC→AGC) | cysK | cysteine synthase A, O‑acetylserine sulfhydrolase A subunit |
| Reads supporting (aligned to +/- strand): ref base G (7/12); new base A (13/0); total (20/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.46e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTCAGCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAAC > NC_000913/2532838‑2533122 | ttCAGCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCAACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCCGGCTGGCGTCCGCCGCCCCATTAAAAG‑AACCcaa > 2:189338/1‑150 (MQ=255) tCAGCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGA‑CGGCGGCTGGGGCCAGCCCCTCCATTAAAAG‑AACCaaa > 1:391123/1‑149 (MQ=255) gCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCaa < 2:507348/151‑1 (MQ=255) gCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑AGCTGGCTGGGGGCAGCCGCTACATTAAAAG‑CACCaaa > 2:46134/1‑146 (MQ=255) aaTCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTATAGCAG‑CAGGTGGCCGCCGCCACACTAAAAG‑CACCaaa > 1:190324/1‑144 (MQ=255) aTCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCCGGTTGGGGTCCGCCGCGACATTTAAAG‑CACCaaa > 2:577041/1‑143 (MQ=255) ccGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑Acc < 1:186857/150‑1 (MQ=255) cGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑AGCCGGCTGGGGTCAGCCCCCCCATTAAAAG‑GACCCAAGGCaa > 1:47829/1‑145 (MQ=255) cGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGA‑AGCCGACTTGGGTCCGCCGCCCCATTTAAGG‑GCCCAAAAGCaaaaac > 1:96998/1‑145 (MQ=255) cccTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTAtc < 2:201291/151‑1 (MQ=255) aaTTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGA‑AACCGGATGCGGGCCCGCGCTACATTTAAGGACACCAAAGGAAAG‑Accagctttttt > 1:504275/1‑137 (MQ=255) tCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTTACTGGCGGCAGCCGCCCCATTAAAAA‑GACCAAAAGCAAG‑ACCCAGCTTATTTTTTTCGc > 2:203796/1‑150 (MQ=255) aaaaGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑AGCTGACTGGGGGCAACCGCCACATTAAAAG‑CACCAAAAGCAAA‑AACCAACTTTTCTCTTTCGcc > 2:294549/1‑146 (MQ=255) aaaGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCCGACTGGGGGCAGCCCCTCCCTTTAAAG‑CACCAAAAGCAAA‑ACCGAACTTTTCTCCTTCGcc > 1:247802/1‑145 (MQ=255) aGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑AGCTGAATGGGGGCCGCCGCCACATTAAAAG‑CACCAAAAGCAAA‑AACCATCTTATCTTTGTTGCCGTTGAGcc > 2:36407/1‑151 (MQ=255) ggAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CCCTAGATGGGGGCAGGCCCTCCAATAAAAG‑CACCAAAAGGAAGAACCGACCTT‑TTTCTTTCGCCGGTTAGGCAACCGCATCTcc > 2:65033/1‑150 (MQ=255) gagaTATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAg < 1:274066/151‑1 (MQ=255) gagaTATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAg < 2:398007/151‑1 (MQ=255) atatGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTa < 2:533079/151‑1 (MQ=255) gATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAgg < 1:85928/151‑1 (MQ=255) ccGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑AGCTGACTGGGGGCCGCCGCCACATTTAAAA‑AACCAAAAGCAAA‑ACCGAACTTATTTTTTTCGCCCTTGAGCCAACCCATTCTCCAGTTATTGcccccggg > 2:512856/1‑144 (MQ=255) gTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAgg < 1:521937/151‑1 (MQ=255) gTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAgg < 2:574043/151‑1 (MQ=255) tatttattGCTGGCGTTGGGACTGGCGGTA‑AGCTGACTGGCGTCCGCCGCGACATTAAAGG‑CAACAAAGGCAAG‑ACCGACATTTTCTCCGTTGCCGGTTAGCCAACCGATTTTCCCGTTATCGCCCAGGCGCTGGCCGGTGgagaga > 1:51459/1‑147 (MQ=255) tttattGCTGGCGTTGGGACTGGCGGTA‑CGCCGACTGGCGGCAGCCGCTACATTAAAAG‑CATCAAAAGCAAA‑ACCGATCTTTTCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGCTATCGCCCAGGCGCCGGCAGGGGAAGAGATTaaa > 2:358674/1‑150 (MQ=255) cGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTTTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGt > 2:230008/1‑151 (MQ=255) gaa‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGGAt > 1:305211/3‑151 (MQ=255) cGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTgg > 2:284215/1‑151 (MQ=255) ccGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCg < 2:423274/151‑1 (MQ=255) aCATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTaa < 1:515830/151‑1 (MQ=255) aCATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTa < 1:6526/150‑1 (MQ=255) cATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATACCGGCTAAc > 1:307580/1‑151 (MQ=255) | TTCAGCAATCCGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTA‑CGCTGACTGGCGTCAGCCGCTACATTAAAGG‑CACCAAAGGCAAG‑ACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCGCTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAAC > NC_000913/2532838‑2533122 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |