Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 802,799 C→G 9.8% intergenic (‑62/+90) ybhC ← / ← ybhB acyl‑CoA thioesterase, lipoprotein/kinase inhibitor homolog, UPF0098 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273802,7990CG 9.8% 366.0 / 24.2 109intergenic (‑62/+90)ybhC/ybhBacyl‑CoA thioesterase, lipoprotein/kinase inhibitor homolog, UPF0098 family
Reads supporting (aligned to +/- strand):  ref base C (72/26);  new base G (4/6);  total (76/33)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 139 total aligned reads displayed.

TCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGC  >  CP009273/802734‑802866
                                                                 |                                                                   
tCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGtt                                                                   <  1:29882584/68‑1 (MQ=255)
  acaGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTgg                                                                 <  1:17675007/68‑2 (MQ=255)
  acaGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTgg                                                                 <  1:20366901/68‑2 (MQ=255)
    aGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggca                                                               >  1:32242767/1‑65 (MQ=255)
       ttGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggc                                                                <  1:10682462/64‑3 (MQ=255)
       ttGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggc                                                                <  1:28759982/64‑3 (MQ=255)
       ttGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAAGTTggc                                                                <  1:27073399/64‑3 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCtt                    >  1:19166730/4‑68 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCtt                    >  1:6584953/4‑68 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCtt                    >  1:24936711/4‑68 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCt                     <  1:3550727/64‑1 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCt                     <  1:22628683/64‑1 (MQ=255)
                                               tcagCTTTTTTATACTAACTTGAGCGAAACGGGAAGATAAAAAGACAAAAAGTTGTTTTTAATACCt                     <  1:21947118/64‑1 (MQ=255)
                                                  gCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTaa                 >  1:30776610/1‑68 (MQ=255)
                                                   cTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                >  1:34487446/1‑68 (MQ=255)
                                                   cTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                >  1:10551937/1‑68 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTaa                 <  1:12653225/66‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTaa                 <  1:12544647/66‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                <  1:9570362/67‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                <  1:21548439/67‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                <  1:9234802/67‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAg                <  1:11792256/67‑1 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGt               >  1:19490718/1‑68 (MQ=255)
                                                    ttttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGt               >  1:7523124/1‑68 (MQ=255)
                                                     tttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTg              <  1:19311092/68‑1 (MQ=255)
                                                     tttttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTg              <  1:24292823/68‑1 (MQ=255)
                                                       tttATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGAt            >  1:28135197/1‑68 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:6965934/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:25692975/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:31989631/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:28754957/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:31980127/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:26534627/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:23270241/1‑67 (MQ=255)
                                                         tataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATa           >  1:22419038/1‑67 (MQ=255)
                                                          attcTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAcc         <  1:10221450/65‑1 (MQ=255)
                                                          ataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAcc         <  1:22898535/68‑1 (MQ=255)
                                                          ataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAcc         <  1:25489294/68‑1 (MQ=255)
                                                          ataCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATAcc         <  1:21427647/68‑1 (MQ=255)
                                                           taCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCa        >  1:21444333/1‑68 (MQ=255)
                                                            tatAATTTGAGCCTAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       <  1:1668243/66‑1 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:31523175/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:13155529/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:31523190/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:3063050/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:30521948/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:29421646/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:13606803/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:13056752/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:32991791/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:33101211/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:33643111/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:33677645/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       <  1:33807569/68‑1 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:34330524/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:12369161/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:4971110/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:5870141/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:5972360/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:11272446/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:9018005/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:25828932/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:20641294/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:23687339/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:14780708/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       <  1:25154898/68‑1 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:20521706/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:25588898/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:25633437/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       <  1:18459993/68‑1 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:14737580/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:14933642/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:20809763/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:21001757/1‑68 (MQ=255)
                                                            aCTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAg       >  1:15720524/1‑68 (MQ=255)
                                                             cTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGa      >  1:27338425/1‑68 (MQ=255)
                                                                aCTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATg    >  1:21027905/1‑67 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  <  1:7836298/68‑1 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:6694240/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:22777719/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:808002/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:8584884/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:17317223/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:11219997/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:9235794/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:17688352/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:4273149/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:22524165/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:23291995/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:23406015/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:32391264/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:32380571/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:31633245/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:28219061/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:13991150/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:29401420/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:30454195/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGc  >  1:30356940/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTAATACCAGATGGc  >  1:7293070/1‑68 (MQ=255)
                                                                 cTTGAGCGAAACGGCAACGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTCATACAAGATCGc  >  1:5250616/1‑68 (MQ=255)
                                                                 |                                                                   
TCACAGGTTGCTCCGGGCTATGAAATAGAAAAATGAATCCGTTGAAGCCTGCTTTTTTATACTAACTTGAGCGAAACGGGAAGGTAAAAAGACAAAAAGTTGTTTTTAATACCTTTAAGTGATACCAGATGGC  >  CP009273/802734‑802866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: