Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 775,120 G→A A123A (GCG→GCA ybgC → acyl‑CoA thioester hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913775,1200GA98.0% 153.4 / ‑4.1 50A123A (GCG→GCAybgCacyl‑CoA thioester hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/28);  total (21/29)
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.76e-01

ACAGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCAT  >  NC_000913/774985‑775253
                                                                                                                                       |                                                                                                                                     
aCAGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACtt                                                                                                                                    <  1:561588/139‑1 (MQ=255)
  aGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTcc                                                                                                                                  >  1:761942/1‑139 (MQ=255)
                tCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcg                                                                                                                    >  1:544211/1‑139 (MQ=255)
                 cAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  2:597773/139‑1 (MQ=255)
                  aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgcg                                                                                                                  <  1:391309/139‑1 (MQ=255)
                  aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgcg                                                                                                                  <  1:391310/139‑1 (MQ=255)
                  aaTGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGAGCACTTCCCAAGTCTATTGTcgcg                                                                                                                  >  2:446351/1‑139 (MQ=255)
                         ggCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   <  1:128958/131‑1 (MQ=255)
                         ggCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTcgc                                                                                                                   >  2:128958/1‑131 (MQ=255)
                             ccTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCTCGGAGTTTAAGCa                                                                                                       <  1:536852/139‑1 (MQ=255)
                                 tttGGTTTTCACGTAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   <  1:528928/139‑1 (MQ=255)
                                 tttGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   <  1:220981/139‑1 (MQ=255)
                                  ttGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtgac                                                                                                  <  2:761942/139‑1 (MQ=255)
                                        ttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACAt                                                                                            >  2:76284/1‑139 (MQ=255)
                                        ttCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACAt                                                                                            >  1:665591/1‑139 (MQ=255)
                                             gCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTATGCGTTGACCCACTCAATATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAata                                                                                       <  1:488856/139‑1 (MQ=255)
                                                             gCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCt                                                                       >  2:670511/1‑139 (MQ=255)
                                                                          tGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          <  2:105776/139‑1 (MQ=255)
                                                                          tGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          <  2:544211/139‑1 (MQ=255)
                                                                          tGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          <  2:752421/139‑1 (MQ=255)
                                                                           gCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTc                                                         <  1:209156/139‑1 (MQ=255)
                                                                           gCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGAAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCATTc                                                         >  2:468667/1‑139 (MQ=255)
                                                                              gaatgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTgg                                                      >  2:727695/1‑139 (MQ=255)
                                                                               aatgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGt                                                     <  1:670511/139‑1 (MQ=255)
                                                                                 tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          >  2:319946/1‑132 (MQ=255)
                                                                                 tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCtt                                                          <  1:319946/132‑1 (MQ=255)
                                                                                 tgaaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTa                                                   >  1:149442/1‑139 (MQ=255)
                                                                                   aaGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTaaa                                                 <  2:268234/139‑1 (MQ=255)
                                                                                      cAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   <  1:782820/86‑1 (MQ=255)
                                                                                      cAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGtga                                                                                                   >  2:782820/1‑86 (MQ=255)
                                                                                        gagGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTAt                                            <  2:261496/139‑1 (MQ=255)
                                                                                         agGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATc                                           >  2:221302/1‑139 (MQ=255)
                                                                                           gTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTa                                                   <  2:468784/129‑1 (MQ=255)
                                                                                           gTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTa                                                   >  1:468784/1‑129 (MQ=255)
                                                                                           gTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCAt                                         >  2:352910/1‑139 (MQ=255)
                                                                                                     tttGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                         <  1:558604/81‑1 (MQ=255)
                                                                                                     tttGCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGaa                                                                                         >  2:558604/1‑81 (MQ=255)
                                                                                                        gCGTTGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGttgattttgatt                            <  1:637169/139‑1 (MQ=255)
                                                                                                           ttGACCCACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGt                         >  2:710697/1‑139 (MQ=255)
                                                                                                               cccACTCAAAATGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGttttt                     >  1:359928/1‑139 (MQ=255)
                                                                                                                          tGAAGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCt          <  2:221024/139‑1 (MQ=255)
                                                                                                                             aGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTgg       >  1:562357/1‑139 (MQ=255)
                                                                                                                             aGCCTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTgg       >  2:61521/1‑139 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTTACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:61521/139‑1 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATTAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:483781/139‑1 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTGAATCGCATCTTGGGc     <  1:318546/139‑1 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:727695/139‑1 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:221302/139‑1 (MQ=255)
                                                                                                                               ccTCGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGc     <  1:299604/139‑1 (MQ=255)
                                                                                                                                  cGTGCACTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCat  >  2:494329/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                     
ACAGACTGAAATAACATCAATGCGTGGCACCTCTTTGGTTTTCACGCAACGTATTGTCAACGCCGAGAATACTTTGCTGAATGAAGCAGAGGTTCTGGTTGTTTGCGTTGACCCACTCAAAATGAAGCCTCGTGCGCTTCCCAAGTCTATTGTCGCGGAGTTTAAGCAGTGACTGACATGAATATCCTTGATTTGTTCCTGAAGGCTAGCCTTCTGGTTAAACTTATCATGTTGATTTTGATTGGTTTTTCAATCGCATCTTGGGCCAT  >  NC_000913/774985‑775253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: