| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,028,160 | G→C | 16.4% | intergenic (‑143/+28) | dgcQ ← / ← yodC | putative diguanylate cyclase DgcQ/protein YodC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,028,160 | 0 | G | C | 16.4% | 44.9 / 3.2 | 23 | intergenic (‑143/+28) | dgcQ/yodC | putative diguanylate cyclase DgcQ/protein YodC |
| Reads supporting (aligned to +/- strand): ref base G (10/9); new base C (2/2); total (12/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01 | |||||||||||
ACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAG > NC_000913/2028025‑2028302 | acgcgcAACCTGGCTGTGTAGGATTCGGGCATAAACGGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGc < 1:427895/138‑1 (MQ=255) tGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGCTAAGACGCGACTATCATTCATTGCAGTCATTAGCACACCAGGCCGGATGTGATCAGGCCCTATGGTaa > 2:83113/1‑143 (MQ=255) ttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTCGCACACCAGGCCGGATGTGATCCGGCCTGATGTTCATGAtc > 1:111915/1‑142 (MQ=255) ttATGATTCTGGCATAACCTGGCGACAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAcc < 2:270916/143‑1 (MQ=255) aTTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCCTTAGCACCCCAGGCCGGATGTGATCCGGCCTGCTGATAATGATccccctcc > 2:230582/1‑139 (MQ=255) cgggCATAGCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca < 2:204882/141‑1 (MQ=255) gCATAACCTGGCGAGAGCGTACCGCAAAACAGCGATAAGTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt < 2:471871/143‑1 (MQ=255) ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCTGAAGAACGActcgcct > 1:546056/1‑143 (MQ=255) cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTGCTTCCGCCGAACGAcccccctc > 2:428773/1‑143 (MQ=255) aaaCAGCGATAATTTTATAATGGTGGCATGATTAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc < 1:525571/143‑1 (MQ=255) cAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa < 1:488608/143‑1 (MQ=255) cAGCGATAATTTTATAATGGTCGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCCTTCTTCCGCCGAACGACCCCCCTCCCCCgcg > 2:579903/1‑141 (MQ=255) aGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa > 1:367091/1‑143 (MQ=255) actggtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGACACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa < 2:301906/141‑1 (MQ=255) aTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt < 1:510970/143‑1 (MQ=255) gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAg > 1:479325/1‑143 (MQ=255) tGAGAAGACGCTCCTATCATGCGTTGCAGGCATTAGCGCACCCGGCTGGATGTCATCCGGCCTCATCCTAACGATCCCCTTTCACACTTCTTCCGCAGAACGTCTCCCCTCCCCCGTAACAAGCTCGGTTTCATGAAAAGCtt < 2:561805/143‑1 (MQ=255) aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATATCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc > 1:345375/1‑143 (MQ=255) ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa > 1:68981/1‑143 (MQ=255) acCAGGCCGGGTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTCCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCACCGtcattaaa < 2:268819/143‑6 (MQ=255) tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACACGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCACCGACACTcccc > 2:652688/2‑140 (MQ=255) tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAAc > 2:298659/2‑143 (MQ=255) ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc < 2:367091/143‑1 (MQ=255) gATCCTGCCTCATTTCAATGATACCCTTCCACACATCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAg < 1:136371/143‑1 (MQ=255) | ACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAG > NC_000913/2028025‑2028302 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |