Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 32,653 | 0 | T | C | 47.4% | 41.7 / 13.9 | 38 | S613P (TCC→CCC) | carB | carbamoyl‑phosphate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base T (17/3); new base C (0/18); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.02e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC > NC_000913/32508‑32802 | gCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGa < 2:607104/150‑1 (MQ=255) gcCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACtt > 1:440248/1‑151 (MQ=255) gAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCg > 2:322267/1‑149 (MQ=255) cGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAAct > 1:383649/1‑151 (MQ=255) gttCAGGTTATAAAATTCGACTACTGTGGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAACCCATTATGGTTAACTGTACCCC‑GGAAC‑CCGTCCCCCCCGACTACGACACTTCCGCCCGCCTCTACTTCGAGCCGGTAAc < 2:579051/147‑1 (MQ=255) ggccAGGGTATCGAATTCGACTCCTTTTGCGTACCCGCCTCGCTGGCCCTGCGCGAAGCCGTTTACAAAACCTTTAGGTTTACTTGTACCCC‑GAAAC‑CCGTCTCCCCCGACTCCGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAAct < 2:214036/147‑1 (MQ=255) gTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAActct > 1:69653/1‑151 (MQ=255) gggTATGAAATTCGACTACTTTTGCGTCCCCGCCTCGCTGGCGCTGCGCGAAGACGGTTACAAAACCATTATGTTAACCTGTACCCC‑GAAAC‑CCGTCCCCCCCGCCTACAACACTCCCGACCGCCTCTACTTCGAGCCGGTAActct < 1:67673/147‑1 (MQ=255) gggTATCGATTTCGCTTTCTTTTGGGTTCCCGCCTGGCCGGCGCCGGCCGAGGGGGGTTCCCAACCCTTATTGGTAAACTGACACCA‑GAAAC‑CCGCCCCCACCGACTACCGCATTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGa < 2:315623/150‑1 (MQ=255) atttCGACTACTGTTGCTTCCCCGCCTCGCTGGCGCTGCGCAAAGACGTTTACGAACCCATTATGTTTAACTGTACCCC‑GAAAC‑CCTTCTCCACCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGAtgtg < 2:477882/147‑1 (MQ=255) aTTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTg > 1:320137/1‑151 (MQ=255) cGACTACTGTTCCGTACCCGCCTCGCTGGCGCTGCGCGAGGCCGTTTACAAACCCTTTATGTTTAACTGTACCCC‑GGAAC‑CGGTCTCCCCCGACTACGACCCTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGa < 2:381006/150‑1 (MQ=255) tGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGc > 2:160225/1‑151 (MQ=255) gCGTACACGCCGCGCGGGCGCTGCGCGAAGGCGGTTACGAACCCATTATGGTAAACTGTAACCC‑GAAAC‑CCGTCTCCCCCGACTAAGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGt < 1:167964/149‑1 (MQ=255) gTACCCCCCCCGCGGGCGCTGCGCAAAGACGTTTACAAACCCTTTATGTTTAACTGTACCCC‑GAAAC‑CCGTCCCCCCCGACTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg < 2:320170/151‑1 (MQ=255) cccGCCTCGCGGGCGCTGCGCGAAGCGGTTTACAAACCCATTATGTTTAACTGT‑ACCC‑GGAAACCCGTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTAt < 1:542587/144‑1 (MQ=255) cccGCCCCGGTGGCGCTGCGCGAAGGGGGTTACAAAACCATTTTGTTTAACTGT‑ACCC‑GGAAACCCGTCTCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCg < 1:372515/146‑1 (MQ=255) acGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCgaga > 1:365767/1‑150 (MQ=255) gCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGc > 2:409825/1‑151 (MQ=255) ggCGCGGCGGGAAGA‑GGTTGAAAAACCCTTTGGGTTAAATGAAACCC‑GGAAACCCGTCCCCCCCGATAACGACACTCCCAACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGaaa < 1:52986/148‑1 (MQ=255) cgAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTc > 2:192826/1‑151 (MQ=255) aaGACGGTTACGAAACCATTATGTTTACCTGTACCCC‑GAAAA‑CCGTCCCCCCCGCCTACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCg < 2:157079/147‑1 (MQ=255) ttACAAAACCTTTATGGTTAACTGTACCCC‑GAAAC‑CCTTCCCCCCCGACTACGACCCTCCGGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATc < 1:173341/139‑1 (MQ=255) tACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcg > 1:319192/1‑151 (MQ=255) aacccATTAGGTTAAACTGT‑ACCC‑GGAAACCCGTCCCCACCGACCACGACACTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCa < 1:670761/148‑1 (MQ=255) aaaCCATTATGTTTACCTGTACCCC‑GGAAC‑CCGTCCCCCCCGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggc < 1:192805/146‑1 (MQ=255) aaaCCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTc > 1:617495/1‑150 (MQ=255) aCCATTATGGTTAACTGTAACCC‑GTAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTAcggcgg > 1:192408/1‑146 (MQ=255) aTTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGAccc > 1:653560/1‑151 (MQ=255) tGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt > 2:532027/1‑151 (MQ=255) gTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCt > 1:93852/1‑149 (MQ=255) ctgaaCCCGGGAAC‑CCGCCCCCCCCGACAACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGaa < 1:6394/144‑1 (MQ=255) cgttcccccccGCCTACGCCCCTCCCGCCCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCt < 1:271360/139‑1 (MQ=255) cccccccGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgct < 2:447304/146‑1 (MQ=255) cccccGACTACGACCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAg < 1:221002/144‑1 (MQ=255) aCCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAgctgct > 2:629563/1‑150 (MQ=255) cggACAACGCCCCTCCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGctgctg < 1:578397/148‑1 (MQ=255) cgacaccgccctttCCGCCCGCCTCTTCTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGc < 2:520439/139‑1 (MQ=255) acTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGGAGCTGCTGGCGTACCGGt > 1:117989/1‑149 (MQ=255) ttCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATc > 2:568454/1‑151 (MQ=255) | GCTTGGCGGCGGCCCGAACCGTATCGGTCAGGGTATCGAATTCGACTACTGTTGCGTACACGCCTCGCTGGCGCTGCGCGAAGACGGTTACGAAACCATTATGGTTAACTGTAACCC‑GGAAA‑CCGTCTCCACCGACTACGACACTTCCGACCGCCTCTACTTCGAGCCGGTAACTCTGGAAGATGTGCTGGAAATCGTGCGTATCGAGAAGCCGAAAGGCGTTATCGTCCAGTACGGCGGTCAGACCCCGCTGAAACTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATC > NC_000913/32508‑32802 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |