Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 673,246 G→A 7.0% A361A (GCG→GCA)  fadI → 3‑ketoacyl‑CoA thiolase FadI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122673,2460GA 7.0% 165.6 / 3.7 57A361A (GCG→GCA) fadI3‑ketoacyl‑CoA thiolase FadI
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base A (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGC  >  ACB122/673127‑673379
                                                                                                                       |                                                                                                                                     
cTGACGCTGCAGGACTTTGACAACTCCGAAATCCCCGAAGCCTTTGCGGCCCAGGTGCTGGGCACGCTCAAGGCCTGGGAATATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATc                                                                                                                        <  2:152245/135‑1 (MQ=255)
 tGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGAtag                                                                                                                       >  1:87289/1‑133 (MQ=255)
     gcggCAGGACTTTGACTACTACGAAAGCCACGAAGCCTTTGCGGCCCAGGGGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACGGCCGTGAACGGTGGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCg                                                                                                                   <  6:340294/132‑1 (MQ=255)
       tGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGCACGCCCgca                                                                                                                 >  1:128827/1‑135 (MQ=255)
       tGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCgca                                                                                                                 <  5:188655/135‑1 (MQ=255)
       tGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCgca                                                                                                                 >  5:96410/1‑135 (MQ=255)
       tGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCgca                                                                                                                 <  6:96410/135‑1 (MQ=255)
           ggACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGCGGCGCCGCTTGGCGCGCTCGACAGCAGCaa                                                                                                             >  8:358576/1‑135 (MQ=255)
            gACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCACCGCTTGGCGCGATCGACCGCAGCaaa                                                                                                            >  5:356515/1‑135 (MQ=255)
            gACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCACCGCTTGGCGCGATCGACCGCAGCaa                                                                                                             >  5:362400/1‑134 (MQ=255)
             aCTTTGACTACTACGAAATCCACGAAGCGTTTGCGGCCCAGGTTCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGTGGCGCCGCTTGGCGCGATCGACCGCAGCAAAc                                                                                                           <  6:325854/135‑1 (MQ=255)
                    ctactaCGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATg                                                                                                     >  7:341125/1‑134 (MQ=255)
                    ctactaCGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCAGCTTTTCGCGGAAGACCGCACCAAACACAAAGt                                                                                                    >  7:396735/1‑135 (MQ=255)
                    ctactaCGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCAGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATg                                                                                                     >  3:100373/1‑134 (MQ=255)
                        tacGAAATCCACGAAGCCTTTGCGGCCCAGGGGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAg                                                                                                 <  8:357417/134‑1 (MQ=255)
                           gAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTTGGCCTGGAGGCGCCGCTTTGCGCGATAGACCGCCGCAAACTCAAAGTCAggg                                                                                               >  5:455111/1‑133 (MQ=255)
                           gAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGCCCGCAGCCAACTCAATGTCAGGGGc                                                                                             >  2:171112/1‑135 (MQ=255)
                                ccccGACGCCTTTGCGGCCCAGGTGCTGTGCACGCTCACGGCCGGGGAAGATGCCTATTCCTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTc                                                                                        <  5:236886/132‑1 (MQ=255)
                                     aaGCCTGTGCGGCCCAGGTGCTGGGCACGCTCAAGGCGTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCGGGCGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTgg                                                                                   <  4:422007/135‑1 (MQ=255)
                                             gCGGCCCAGGGGCTGTGCACGCGCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTGGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCgg                                                                            <  4:215449/134‑1 (MQ=255)
                                             agacagcAGGTGCTGTGCACGCGCAAGGCCTGGGACGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCggg                                                                           <  8:435375/129‑1 (MQ=255)
                                              gggCCCTGGTGCTGTGCACGCTCAAGGCGTGGGAAGATGCCGATTACTGCCGTTACCGGTTGGGCCTGGTGGCGCCGCTGGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCggg                                                                           <  5:401917/133‑1 (MQ=255)
                                                     ggTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCCCTTGGGGCGATCCCCCGCAGCAAACTAAATCTGAGGGGCCGTTCGCTGGCGGCCCGGCAccccgt                                                                   >  2:204111/1‑132 (MQ=255)
                                                              cACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCCAACTCAATGTCAGGGGCAGTTCGCGGGCGGCCGGGCACCCCGTTGCGGCAAc                                                          >  5:251500/1‑135 (MQ=255)
                                                               aCGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGCTCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAAcc                                                         >  6:112790/1‑135 (MQ=255)
                                                               aCGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGCCCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAAcc                                                         >  7:261400/1‑135 (MQ=255)
                                                                 gCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGCACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGAAAGCCCGCAGCAAACACAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGGTTTCGGccagtat                                                       >  7:435375/1‑129 (MQ=255)
                                                                  cTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGGGGGGCCGGGCGGCGCCGGTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGt                                                      <  2:252119/135‑1 (MQ=255)
                                                                   tCACGGCCTGGGAAGATGCCGAGTACTGCCGTGAACGGTGGGGCCTGGAGGCGCCGCTGGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTg                                                     <  5:212534/135‑1 (MQ=255)
                                                                   tCAAGGCGTGGGAAGATGCGGTTTACTGCCGTGAACGGTTGGGCGGGGAGGCGCCGGTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTg                                                     <  2:128827/135‑1 (MQ=255)
                                                                   tCAAGGCCTGGGATGATGCCGATTACTGGCGTGAACGGTGCGGCCTGGAGGGGCCGTTTTGCGCTTTCGACCGCAGCAAACTCAATGTCTGGTGCAGTTCGGTGGCGGCCGGGCTCCCGTTTGCGGCATCCGGTg                                                     <  2:97527/135‑1 (MQ=255)
                                                                    cAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGCCCGCAGCAAACTCAATGTCAGGGGCACTTCGCTGGCGGCCGGGCCCCCGTTTGCGGGCACCGGTgg                                                    >  6:324006/1‑135 (MQ=255)
                                                                        gCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCgg                                                <  5:204330/135‑1 (MQ=255)
                                                                         ccTGGGAAGATGCCGATTACTGCCGTGCACGGTTGGGCCTGGAGGCGCCGCGTTGGGCGCACCTCCGCCGCCAACGCCATGTGACGGGGCGTTCGCGGGGGGCGGGGCCCCCGGGTGCGGGACCCCGTggggggg                                               >  6:297326/1‑134 (MQ=255)
                                                                         ccTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGc                                                                <  4:71475/118‑1 (MQ=255)
                                                                         ccTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGc                                                                >  3:71475/1‑118 (MQ=255)
                                                                              gAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGCCCGCACCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCCCCCGTTTGCGGCCACCGGGGGGGGGAttttt                                          >  5:292173/1‑134 (MQ=255)
                                                                                aGATGCCGATTACTGCCGTGAACGGTTTGGCCTGGAGGCGCCGCTTGGGCCGATACACCGAAGCAAACTCAATGTCAGGGGCAGTTGGCTGGCGGCCGGGCACCCCTTTGCGGCAACCGGTTGGCGTAtttttg                                         >  2:71228/1‑132 (MQ=255)
                                                                                aGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTgg                                         >  3:207744/1‑134 (MQ=255)
                                                                                aGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTTGc                                        >  6:103129/1‑135 (MQ=255)
                                                                                 gATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCACCGCTTGGCGCGATCGACCGCAGCAAACTCACTGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGcc                                       <  6:362400/135‑1 (MQ=255)
                                                                                 gATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCACCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGcc                                       <  6:356515/135‑1 (MQ=255)
                                                                                     ccGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAAcc                                   >  4:82438/1‑135 (MQ=255)
                                                                                      cGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCCACCGGTGGGCGGATATTGGCCAACCt                                  >  6:142869/1‑135 (MQ=255)
                                                                                       gATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTg                                 >  6:262713/1‑135 (MQ=255)
                                                                                          tACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCCACCGGGGGGCGGAAATTGGCCACCCTggg                               >  8:127391/1‑133 (MQ=255)
                                                                                                  tGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGtt                      >  7:280042/1‑135 (MQ=255)
                                                                                                  tGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGtt                      <  1:104443/135‑1 (MQ=255)
                                                                                                     aCGGGTGGGCCTGTAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGc                   <  1:204111/135‑1 (MQ=255)
                                                                                                       gggtGGGCGTGTAGGGGACGTTTGGGGCTATCTACAGCAGAAAACTCTATGTCAGGGGCAGTTCGCTGGCGGCCGGGCTCCCGTTTGCGGCAACCGGTGGGCGGATTTTGGCCAACCTGGCGAAGATGTTGGCaa                 <  2:176342/132‑1 (MQ=255)
                                                                                                       ggTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCaa                 <  4:207744/135‑1 (MQ=255)
                                                                                                       ggTGGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCaa                 <  2:97017/135‑1 (MQ=255)
                                                                                                            ggCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCACcggcg            >  3:82757/1‑135 (MQ=255)
                                                                                                              ccTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCACACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCCACCTGGCCAAGCTGTTGGCATCGGCggg          >  8:260282/1‑135 (MQ=255)
                                                                                                              ccTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCggg          <  3:127054/135‑1 (MQ=255)
                                                                                                                 ggAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTaa       >  1:349516/1‑135 (MQ=255)
                                                                                                                    ggCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAggg    <  6:152342/135‑1 (MQ=255)
                                                                                                                      cgccgcTTGTCGGTATCGACCGCAGAAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTTGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAgggg   <  1:71228/134‑1 (MQ=255)
                                                                                                                      cgccgcTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGgc  >  4:382111/1‑135 (MQ=255)
                                                                                                                       |                                                                                                                                     
CTGACGCTGCAGGACTTTGACTACTACGAAATCCACGAAGCCTTTGCGGCCCAGGTGCTGTGCACGCTCAAGGCCTGGGAAGATGCCGATTACTGCCGTGAACGGTTGGGCCTGGAGGCGCCGCTTGGCGCGATCGACCGCAGCAAACTCAATGTCAGGGGCAGTTCGCTGGCGGCCGGGCACCCGTTTGCGGCAACCGGTGGGCGGATATTGGCCAACCTGGCCAAGCTGTTGGCAACGGCGGGTAAGGGGC  >  ACB122/673127‑673379

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