Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,279,619 T→A 7.7% L239Q (CTG→CAG)  fadH → 2,4‑dienoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,279,6190TA 7.7% 154.6 / 4.8 52L239Q (CTG→CAG) fadH2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base A (2/2);  total (27/25)
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.54e-01

GTTCCTGGCCGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGCGGCCTTCA  >  NC_002947/2279470‑2279765
                                                                                                                                                     |                                                                                                                                                  
gTTCCTGGCCGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCt                                                                                                                                                    >  2:586522/1‑150 (MQ=255)
     tGGCCGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCg                                                                                                                                               >  1:519740/1‑150 (MQ=255)
     tGGCAGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCg                                                                                                                                               >  1:514682/1‑150 (MQ=255)
      ggCCGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGa                                                                                                                                              >  1:398059/1‑150 (MQ=255)
       gccgccCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAg                                                                                                                                             >  1:1066965/1‑150 (MQ=255)
              acacCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGgca                                                                                                                                      <  2:607024/150‑1 (MQ=255)
              acacCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGgca                                                                                                                                      >  2:853965/1‑150 (MQ=255)
              acacCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGgca                                                                                                                                      <  2:519740/150‑1 (MQ=255)
               cacCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGgcag                                                                                                                                     >  2:145441/1‑150 (MQ=255)
                accaaccaCCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGgcagc                                                                                                                                    >  1:926550/1‑150 (MQ=255)
                       accgcaccgACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCAGGTCGAgcatcga                                                                                                                                       >  1:1083428/1‑134 (MQ=255)
                       accgcaccgACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCAGGTCGAgcatcga                                                                                                                                       <  2:1083428/140‑7 (MQ=255)
                           caccgACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACga                                                                                                                         >  2:456578/1‑150 (MQ=255)
                             ccgaccgCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACgaga                                                                                                                       <  1:939716/150‑1 (MQ=255)
                                   gCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAgc                                                                                                                 <  1:805050/150‑1 (MQ=255)
                                   gCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAgc                                                                                                                 >  2:744144/1‑150 (MQ=255)
                                         gcggcAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTgg                                                                                                           <  1:744144/150‑1 (MQ=255)
                                          cggcAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGc                                                                                                          >  1:764821/1‑150 (MQ=255)
                                              aGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAg                                                                                                      >  1:18476/1‑150 (MQ=255)
                                                ttATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGc                                                                                                    <  1:145441/150‑1 (MQ=255)
                                                tgaTGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGc                                                                                                    >  2:1109618/3‑150 (MQ=255)
                                                    gAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATc                                                                                                >  2:1120039/1‑150 (MQ=255)
                                                    gAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATc                                                                                                >  2:607123/1‑150 (MQ=255)
                                                      aaaCCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGa                                                                                              >  1:1043012/1‑150 (MQ=255)
                                                      aaaCCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGa                                                                                              <  2:160657/150‑1 (MQ=255)
                                                         ccGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCa                                                                                           <  2:18476/150‑1 (MQ=255)
                                                             aTGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGcc                                                                                       >  2:1009578/1‑150 (MQ=255)
                                                                              aaTCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGAc                                                                                                                           <  1:813766/97‑1 (MQ=255)
                                                                              aaTCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGAc                                                                                                                           >  2:813766/1‑97 (MQ=255)
                                                                                 cGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAacac                                                                   <  1:586522/150‑1 (MQ=255)
                                                                                   tCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCg                                                                 <  1:456578/150‑1 (MQ=255)
                                                                                           gTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGt                                                         <  2:1131912/150‑1 (MQ=255)
                                                                                           gTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGt                                                         <  1:1120039/150‑1 (MQ=255)
                                                                                           gTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGt                                                         >  2:1151983/1‑150 (MQ=255)
                                                                                           gTGCGCGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGt                                                         <  2:235655/150‑1 (MQ=255)
                                                                                                        tAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGgcgc                                            <  2:872543/150‑1 (MQ=255)
                                                                                                               gAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTcc                                     >  1:690223/1‑150 (MQ=255)
                                                                                                               gAACTTCATCATCAGCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTcc                                     >  2:540245/1‑150 (MQ=255)
                                                                                                                   ttCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGAcc                                 >  1:654342/1‑150 (MQ=255)
                                                                                                                     catcatcatCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCAt                               <  1:25779/150‑1 (MQ=255)
                                                                                                                                   gCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGaa                                                              <  2:200177/105‑1 (MQ=255)
                                                                                                                                   gCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGaa                                                              >  1:200177/1‑105 (MQ=255)
                                                                                                                                        tCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCTCGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCcgcac            <  1:1151983/150‑3 (MQ=255)
                                                                                                                                          gATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGc          >  2:668655/1‑150 (MQ=255)
                                                                                                                                           acccTCGACCAGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTa                                         <  1:48320/115‑1 (MQ=255)
                                                                                                                                           acccTCGACCAGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTa                                         >  2:48320/4‑118 (MQ=255)
                                                                                                                                             gCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGAc                                  >  1:933716/1‑123 (MQ=255)
                                                                                                                                             gCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGAc                                  <  2:933716/123‑1 (MQ=255)
                                                                                                                                             gCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGCGGc       <  2:690223/150‑1 (MQ=255)
                                                                                                                                               tCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCcgcac            <  1:311415/143‑3 (MQ=255)
                                                                                                                                               tCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCcgcac            >  2:311415/1‑141 (MQ=255)
                                                                                                                                                cGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGCGGCCtt    <  1:668655/150‑1 (MQ=255)
                                                                                                                                                  aCCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGCGGCCTTCa  <  2:764821/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTTCCTGGCCGCCCACACCAACCACCGCACCGACCGCTGGGGCGGCAGTTATGAAAACCGCATGCGCCTGGCAGTGGAAATCGTCAGCCGGGTGCGTGGCGCGGTAGGGCCGAACTTCATCATCATCTTCCGCCTGTCGATGCTCGACCTGGTCGAGGGTGGCAGCACCTGGGACGAGATCGAGCTGCTGGCCAAGGCCATCGAGCAGGCCGGCGCGACCTTGATCAACACCGGAATCGGTTGGCACGAGGCGCGTATTCCGACCATCGCCACCAAAGTGCCGCGTGCGGCCTTCA  >  NC_002947/2279470‑2279765

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