| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,036,927 | 0 | T | A | 42.5% | 49.2 / 26.5 | 47 | L60Q (CTA→CAA) | hyaF | hydrogenase‑1 protein nickel incorporation factor |
| Reads supporting (aligned to +/- strand): ref base T (17/10); new base A (0/20); total (17/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.11e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAC > NC_000913/1036779‑1037076 | ccAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTaa < 2:334575/151‑1 (MQ=255) cAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTa > 2:905382/1‑149 (MQ=255) aGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGc < 2:958268/151‑1 (MQ=255) aGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAAACGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGc > 2:603885/1‑151 (MQ=255) cGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAAcc > 1:725094/1‑151 (MQ=255) aCGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAAccgcc < 2:709322/151‑1 (MQ=255) accccACCGCCGATCCCCTGTCAGGGAAATGATGAACCGAGTGTGGCGCCCCGGGAGCAATGTGCTCACACCCGGGAGGTAATTGCGCTTTTAAACGAGTTAAAAC‑TTCAACTAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCt < 2:377207/143‑1 (MQ=255) aaTCCACTGCCGATCCCCTGTCAGGTAAATGATGACCCGAGTAGGGCGCCCCGGGAGAAATGTGCTCACAGCCCGCGGGTGATTGCGTTGTAAAACGATTTACAAC‑ATCAACTAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAg < 2:264962/151‑1 (MQ=255) cccctcccGATCCCCTGTCGGGTGAATGAGGACCCGAGTAGGGCGCCCCGGAAGAAATGTGCCCACACCCCGGAGGTTTTTGCGCTTTAAAACGAGTTACACC‑ATCAACAAAGCAACCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTc < 1:430738/144‑1 (MQ=255) ccccggCCGAGCCCCTGTCGGGTAAATGAGGAGCGGAGTAGGGCGGCCAGGGAGCAATGTGCTCCCACCCCGCGGGGGGTGGGGCTGTAAAGGGGGTTACAAC‑AACAACAAAGCGACCGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTc < 1:565319/145‑1 (MQ=255) ccctcccaaTCCCCTGTCAGGTAAATGATGACCCGAGTATGGCGCCCCTGGAGCAGTGTGCTCCCACCCCGCGGGTGATGGCGTTGTAAAACGATTTACAAC‑ATCAACAAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGa < 2:632862/143‑1 (MQ=255) cACTGCCGATCCCCTGTCAGGTAAATGATGACCGGAGTAGGGGGGCCCGGGAGAAATGTGCTCACAGCCCGCGGGTGATTGCGCTGTAAAACGATTTAAAAC‑ATCAACAAAGCGAACGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCg < 2:787130/150‑1 (MQ=255) tGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCt > 2:353998/1‑151 (MQ=255) tGCCGACCTCCTGTCAGGTAAATGATGACCCGAGTATGGCGGCCCTGGAGCAAGGTGCTCACAGCCCGCGGGTGATTGCGTTTTAAAAGGATTTACAAC‑ATCAACAAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGAt < 2:168237/149‑1 (MQ=255) cAGGTAAATGAAGACCCGAGTTGGGCGCCCATGGAGAATTGTGCTCACACCCCGAGGGTGTTTGCGCTTTTAAGCTAGTTAAAAC‑CTCAACAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAAc < 1:584574/150‑1 (MQ=255) aaTGATGACCCGGGTATGGCGGCCGTGGAGAATTGGCCTCACGCCCCGCGGGTGTGTGCGCTGTAAAACGATTAACAAC‑ATCACCAAAGCAAACCCCACCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCcgac < 1:911341/150‑1 (MQ=255) aaTGATGAACCGAGTATGGCGCCCCTGGAGAAATGTCCTCACACCCCGCAGGTGATGGCGTTGTTAAACGGTTTACAAC‑ATAAACAAAGCGACCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCaa < 2:497596/143‑1 (MQ=255) accggagtggggcggcCCGGGAGAAAGGTCCCCACCGCCCGAGGGTGTTTGCGTTGTAAAAGGATTAACAAC‑ATCAACAAAGCGACCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGAt < 1:943065/136‑1 (MQ=255) ccGATTATGGCGCCCCGGGAAAATTGTGCTCACCCCCCGCGGGTGATTGCGCTGTAAAACGAGTTACAAC‑ATCAACTAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCa < 2:247277/149‑1 (MQ=255) ccGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACt > 2:775620/1‑151 (MQ=255) gAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACttt > 1:800908/1‑151 (MQ=255) ggcggcCCGGAATCAATGTGCTCCCCGCCGGCGGGTGATTGCCCTGTAAAACAAGTACCACC‑ACCAACAAAGCAAACGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCa < 2:867492/149‑1 (MQ=255) gcggcCTTGGAGAAATGTCCTCACCCCCCGGGGGTGTTTGCCCTGTAAAAGGAGTTACAAC‑ATCAACAAAGCGAGCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATa < 2:774715/151‑1 (MQ=255) gcCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGc > 1:850363/1‑151 (MQ=255) ggAGCAATGTGCCCACAGCCCGGGGGTGATTGCGCTTTAAAAGGATATACAAC‑ATCAACAAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCtt < 1:234163/148‑1 (MQ=255) aGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCgg > 1:907808/1‑151 (MQ=255) gtgCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaagg > 2:5280/1‑151 (MQ=255) cacaGCCCGCGGGTGAGTGCGCTGTAAAACGGGTAACAAC‑AAAAAAAAAGCGAACGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGgggaaggggaag < 1:950063/151‑1 (MQ=255) cacaCCCCGCAGGTGATTGCGCTTTAAAACGAGTAAAAAC‑ATCAACTAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaag < 2:134960/145‑1 (MQ=255) cggCCCGGGGGTAATTGCGTTGTAAAGTGAGTTACAAC‑ACCAGTTAAGGGACCGAAAGCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaag < 2:16215/141‑1 (MQ=255) cccGGAGGTGTTTGCGCTTTAAAAGGAGTAACAAC‑ATCAAAAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCa < 1:520154/151‑1 (MQ=255) cccGGAGGTGATTGCGCTTTAAAACGATTAACAAC‑ATCAACAAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaaggggaa < 2:379858/145‑1 (MQ=255) cccGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTc > 2:377982/1‑150 (MQ=255) ccGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTc > 2:133220/1‑149 (MQ=255) gCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTg > 2:615449/1‑151 (MQ=255) gTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTgcgc > 2:299277/1‑150 (MQ=255) ttGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTc > 2:354642/1‑150 (MQ=255) taaaaCGATTAACAAC‑AACAATAAAGCAAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCt < 1:599614/148‑1 (MQ=255) aaacggttaacaac‑aacaCCAAACCGACCCCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTcagcag < 1:503835/136‑1 (MQ=255) aaaCGAGTTACAAC‑ATCAAATAAGCGAACGCCACCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTgac < 1:426832/151‑1 (MQ=255) aaCGATTTACAAC‑ATCAACAAAGCGACCGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTgacg < 2:321172/151‑1 (MQ=255) cGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTgacgac > 2:199997/1‑151 (MQ=255) gattACAACAACCAACAAAGCGAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTgacgac < 1:919424/148‑1 (MQ=255) tACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACagtg > 2:909441/1‑150 (MQ=255) acaac‑aCCAACAAAGCGACCCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACagtga < 2:249478/150‑1 (MQ=255) aac‑aTCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACagtgaa > 1:844795/1‑149 (MQ=255) accaccAACTAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACagtgaaa < 1:911218/145‑1 (MQ=255) tAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAc > 2:846238/1‑150 (MQ=255) | CCAGGAACGCAACCGAACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACCGAGTATGGCGGCCCTGGAGCAATGTGCTCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAAC‑ATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAC > NC_000913/1036779‑1037076 |
| 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: ‑ |