Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,605,176 | T→C | 12.5% | Q163R (CAG→CGG) | fhuF ← | hydroxamate siderophore iron reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,605,176 | 0 | T | C | 12.5% | 68.3 / 2.3 | 33 | Q163R (CAG→CGG) | fhuF | hydroxamate siderophore iron reductase |
Reads supporting (aligned to +/- strand): ref base T (13/14); new base C (2/2); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01 |
GATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGC > NC_000913/4605038‑4605315 | gATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCCCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCt < 1:432823/143‑1 (MQ=255) tCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGGATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCGTCTAGTGCTTGCACAACCGGACCCAGCGCCTGGCGGAt > 2:139888/1‑143 (MQ=255) gCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGCTACCAGTTGATCAACTAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTc < 1:387647/143‑1 (MQ=255) ccTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTcc < 1:524428/143‑1 (MQ=255) ccTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTcc < 1:523297/143‑1 (MQ=255) cggcGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCAtt < 1:203702/140‑1 (MQ=255) ccGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCGTCTAGTGCTTGCACAACCGGAAAGTGCGACAGGCTGATTAACGTGTCCATTCg > 2:548338/1‑143 (MQ=255) ttGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTgcg < 1:280720/143‑1 (MQ=255) cTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTACCGTTTCCATTCGATGCTgcggcg > 1:227036/1‑143 (MQ=255) ttCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGa > 1:402611/1‑143 (MQ=255) aGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGCATGCCGTTGAACTCTCCGGGCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATCAGATGCTGCGGCGAATgtggtg > 2:243102/1‑143 (MQ=255) aCCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGCACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGttt > 2:482394/1‑143 (MQ=255) ccAGTTGATCAAATAACCGGTATTACTCCAGAGACCGTTCCATTTTATTTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGTCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGtttt < 2:118853/143‑1 (MQ=255) aGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAGTGTGGTGTTGCGTTTGTa > 1:195315/1‑143 (MQ=255) aTCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACAGTTCCATGCGACGCTACGGCGAGTGTGGTGTTGCGTTTTTATCGTc > 2:273490/1‑143 (MQ=255) tAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTGTAGTGATTGCACAACCGGAACCAGCGCCTGGCTGATGAACGTTTCCCTCCGATGCTGAGCCGCAGGGGGTGTTGCGGTTTGATCTTcccacac > 2:185680/1‑143 (MQ=255) aCCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTCTATCTTCACTCTCAt > 1:286845/1‑143 (MQ=255) ccGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCCCACACATc > 1:122198/1‑143 (MQ=255) ccGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTGTTTATCTTCACACACCTc > 1:378148/1‑143 (MQ=255) cgTATTACTCCAGCTAAGTGGGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCTGCGGCTGGCTGATTAACGTTTCCTTTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGa < 2:463345/142‑1 (MQ=255) aTTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGAccc < 1:420959/143‑1 (MQ=255) aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGa < 2:180955/143‑1 (MQ=255) ccGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTcc > 1:165473/1‑143 (MQ=255) cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa < 1:561122/128‑1 (MQ=255) cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa < 1:410555/128‑1 (MQ=255) cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa > 2:410555/1‑128 (MQ=255) cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa > 2:561122/1‑128 (MQ=255) cctcCGGTCGCGTCTCGTGCGTGCGCAACCGGAACCAGCGCTTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACAGCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTg < 2:214497/142‑1 (MQ=255) ctcCGGTCGCTTCTAGTGCTGGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGa < 1:326386/143‑1 (MQ=255) cgcgTCTAGTGCTTGCACAAGCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGATGGTGGTGTGGCGTTTTGATCTTCAGACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAg < 2:121624/143‑1 (MQ=255) aGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGCTTTTATCTTCACACCCATCGACCCAGCAACAGGCGACGCCTCCCGTCTCGTGCCACCCAGCACGGaa > 2:207394/1‑143 (MQ=255) aacGGAACGAGAGCCTGGCTGATTAACGCTTCCATGCGATGCTGCGGCGAATGTGGTGTTGCGGTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGa < 1:19441/141‑1 (MQ=255) gcCGGGCTGATTAACGTTCCCATTCGCTGCTGCGGCGAAGCTGGTGTTGCGCTTTTATCTTCACACACATCGATCCACAAACAGGCGTCGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGc < 1:361827/143‑1 (MQ=255) | GATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGC > NC_000913/4605038‑4605315 |
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: ‑ |