Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 2,394,308 | G→A | 9.5% | A85V (GCA→GTA) | moaB ← | molybdopterin cofactor biosynthesis protein MoaB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 2,394,308 | 0 | G | A | 9.5% | 111.2 / 4.0 | 42 | A85V (GCA→GTA) | moaB | molybdopterin cofactor biosynthesis protein MoaB |
Reads supporting (aligned to +/- strand): ref base G (18/20); new base A (2/2); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01 |
TGTACCTGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTCATCTGCTTGTAAT > CP000730/2394169‑2394448 | tGTCCCTGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGc < 2:289078/141‑1 (MQ=255) ccTGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTACAATa > 2:273447/1‑141 (MQ=255) ccTGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTACAATa > 2:286683/1‑141 (MQ=255) tGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATAcc < 2:324444/141‑1 (MQ=255) cTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTg < 1:92972/141‑1 (MQ=255) aCGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAACCAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTccacc > 2:115313/1‑141 (MQ=255) cGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTccacca > 1:165319/1‑141 (MQ=255) gaCAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAACGTGGCCTTACTGACTCACTCGCAACATCAAGTTGTGCAATACCTGTTCCACCAGTc > 1:233660/1‑141 (MQ=255) aCAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCg < 2:226769/141‑1 (MQ=255) cgTGCCAACAGCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATc < 1:138865/141‑1 (MQ=255) aaCATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCACTTTAGTTAAATGTGGTTTTACTGCTGTAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATAtc > 1:64598/1‑141 (MQ=255) tCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATAtcttct < 2:12264/141‑1 (MQ=255) ttCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATAtcttcttc > 1:254446/1‑141 (MQ=255) tCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATAtcttcttct < 1:45609/141‑1 (MQ=255) aaCATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAAGACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTaa > 1:95943/1‑141 (MQ=255) aaaTATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACttcttc > 1:67747/1‑141 (MQ=255) cAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGcgtcgt > 1:38463/1‑141 (MQ=255) cccAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTGACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGc > 1:255955/1‑141 (MQ=255) ccAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCt < 2:170495/141‑1 (MQ=255) ccAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCt < 2:67747/141‑1 (MQ=255) aaGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACt < 1:130081/141‑1 (MQ=255) aaGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACt < 2:243356/141‑1 (MQ=255) gCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACttt > 2:202059/1‑141 (MQ=255) tttAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTcac < 2:183082/141‑1 (MQ=255) tAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACaa > 2:63476/1‑141 (MQ=255) tAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACaa > 2:63463/1‑141 (MQ=255) gTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTACAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAAtt < 1:286683/141‑1 (MQ=255) gTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTACAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAAtt < 1:273447/141‑1 (MQ=255) aaaGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGtat < 2:195653/141‑1 (MQ=255) aaaGTGGTTTTACTGCTTCAATCGTCACATCACGTTGAGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGtat < 2:250117/141‑1 (MQ=255) tGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAAtg < 1:241238/141‑1 (MQ=255) gTTTTACTGTGTCAATCGTCACATCACGTTGTGCAATACCTGTTCAACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCCACTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTTTAAtgtg > 2:158623/1‑141 (MQ=255) gTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAAtgtg > 2:224750/1‑141 (MQ=255) gTTTTACTGCTGCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAAtgtg < 1:153646/141‑1 (MQ=255) ttACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGt > 1:327931/1‑141 (MQ=255) tactacTTCAATCGTCACATAACGTTGTGCAATACCTGTTCCACCAGTTGTAATGATGAGATGCATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTCTTCCATCTTTCACACTTGTATAATGTGCGtc > 1:193114/1‑141 (MQ=255) ttCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGtcactcac > 2:200831/1‑141 (MQ=255) aTCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCa < 2:77570/141‑1 (MQ=255) aTCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCa > 2:330975/1‑141 (MQ=255) cacaTCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTc < 1:334459/141‑1 (MQ=255) acaTCACGGTGGGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTCa < 2:233660/141‑1 (MQ=255) gTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTCATCTGCtt < 2:103015/141‑1 (MQ=255) gCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTCATCTGCTTGTAAt > 1:63795/1‑141 (MQ=255) | TGTACCTGCTACAGCACGAGACAATAATGCACGCGTGCCAACATCTTCAACATAACTCAAATATCTAAACAATTCCCCAAAGCCATCTATCTCTTTAGTTAAAAGTGGTTTTACTGCTTCAATCGTCACATCACGTTGTGCAATACCTGTTCCACCAGTCGTAATGATGACATCAATATCTTCTTCTAACCACTTCTTCACCTGCGTCGTGATGGCTACTTTTTCATCTTTCACAATTGTATAATGTGCGTCACTCACTTCAACGTCATCTGCTTGTAAT > CP000730/2394169‑2394448 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |