| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,004,454 | 0 | T | G | 52.0% | 32.5 / 40.6 | 25 | T82P (ACC→CCC) | yigI | 4HBT thioesterase family protein |
| Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (13/0); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.92e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-06 | |||||||||||
CGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GATTTTTAAAGGCCAG > NC_000913/4004304‑4004574 | cggcgacTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGt < 2:420388/151‑1 (MQ=255) ggcgacTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGtt < 1:714019/151‑1 (MQ=255) ttATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTa < 2:398669/150‑1 (MQ=255) ccTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTaaagggatttccc > 1:622296/1‑140 (MQ=255) gCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGc < 1:673236/151‑1 (MQ=255) aCGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCACCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTATCGCGCTTCTCACGcac > 1:594253/1‑151 (MQ=255) aaCAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGAGTTCCCCCGCCCAcc > 1:74006/1‑150 (MQ=255) aCAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAg < 2:174453/151‑1 (MQ=255) aGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGAACTTCCACCCACACCAGa > 2:91493/1‑150 (MQ=255) gCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACcg < 2:12195/151‑1 (MQ=255) tagtagCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTTAGGTGATTCCCCCACAAAccccaac > 2:274721/1‑140 (MQ=255) cAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGTCGGCGACATc < 2:260924/151‑1 (MQ=255) aaCGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGACTTCCCACCCACACCAAACCGGCGGCGGCATcc > 1:446766/1‑146 (MQ=255) ctgGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGACTTCCCCCCCCCCCCCCACCGGCGGGGGGCGCCCGCGCCGACG‑Cgacgac > 1:656671/1‑148 (MQ=255) cgcAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAAGGACTTTCCCCCCCCACCAAAACCGCGGGGGCAACCAACGCCCCCG‑CCAAGACcccccccc > 2:553044/1‑149 (MQ=255) aaCGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCtt < 2:401467/151‑1 (MQ=255) aCGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGACTTCCCCCCCCCCCCACACCGGCGGCGGCACCCACCGCCCACG‑CCCAGACACCCCCCCGCAAAAAGCttt > 1:686310/1‑151 (MQ=255) ccccATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCccac < 2:564060/150‑1 (MQ=255) aTCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTTAGGGACTTCCCCCCCCCACCCCACCGGGGGCGGGATCCCACGCCCGCG‑CAATGAACCCCCCCCGCCAAAAACTTTTCG‑cccccgttccccccc > 1:137705/1‑136 (MQ=255) gcgATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTAAGGGACTTCCCCCCCACACCACACCGGCGGCGGCAACCCGCGCCCGCG‑CAAATACCCCCCCCCGCCAAAAACTTTGCG‑GCCCGCTGGccaccacc > 1:505955/1‑144 (MQ=255) ccGCTGGCTTAGTTCGTCTTCCCTGATGGTTGCGTGGCGGGTTAAGGTACTGCCCACGCGCACCAGACCGGCGGCGACACCCATCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCTCCATCATTGGCT‑g < 1:686767/147‑1 (MQ=255) cGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACTCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GAtttt < 1:172679/151‑1 (MQ=255) cTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGGTTAGGTACTTCCCCCCCCCACCCCACCGGCGGCGGCATCCAGCGCCCCCC‑CAAAGACCCCCCCCCGCCAAAAACTTTGCGCCCCCGTGGCCCCCCA‑CATTGGCT‑GGttttt > 1:358609/1‑150 (MQ=255) gCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GATTTTTaaa < 1:195464/150‑1 (MQ=255) gCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GATTTTTAAAg < 2:643785/151‑1 (MQ=255) cGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GATTTTTAAAgg < 2:278119/151‑1 (MQ=255) ttCATCTTCACTGATGGTTTCGTGGCGGGGTTAGGTACTTCCCCCCCCCCCCAAACCCGCGGCGGCATCCCACGCCCCCGCCACAGACCCC‑CCCCGCCAAAAGCTTTTCG‑CCCCGGTGCCCCCCCACATTGGGTGGTTTTTTAAAGGccgc > 2:411592/1‑149 (MQ=255) | CGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACG‑CAATGACCCCGCCGTGCAAAATGCTTTGCG‑CCCAGTTGCCCACCATCATTGGCT‑GATTTTTAAAGGCCAG > NC_000913/4004304‑4004574 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |