Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,173,244 T→C 8.9% intergenic (‑285/+47) fadD‑I ← / ← fadD‑II long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,173,2440TC 8.9% 150.0 / 5.4 56intergenic (‑285/+47)fadD‑I/fadD‑IIlong‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base C (2/3);  total (25/31)
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.26e-02

GGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGG  >  NC_002947/5173099‑5173388
                                                                                                                                                   |                                                                                                                                                  
ggCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTgg                                                                                                                                                 >  6:280239/1‑149 (MQ=255)
 gCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTggg                                                                                                                                                <  3:41589/149‑1 (MQ=255)
       aCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAgtgt                                                                                                                                          <  5:243837/149‑1 (MQ=255)
       aCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAgtgt                                                                                                                                          >  4:345504/1‑149 (MQ=255)
       aCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAgtgt                                                                                                                                          <  2:123087/149‑1 (MQ=255)
                aTGGTCACAAAAATGACTCAATGAACAATACAGTTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACACTGCCCAATGTAGGAGCg                                                                                                                                   <  5:98256/147‑1 (MQ=39)
                aTGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACACTGCCCAATGTAGGAGCgg                                                                                                                                  <  6:47030/148‑1 (MQ=255)
                  ggTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGtt                                                                                                                              >  1:147862/1‑149 (MQ=255)
                    tCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTa                                                                                                                            >  1:320965/1‑149 (MQ=255)
                     cacaAAAATGACTCAATGAACAATACAGGTCATTTTTAGATTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGGAGTGTAGGAGCGGGTTTAc                                                                                                                           >  8:133744/1‑149 (MQ=255)
                      acaAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTAcc                                                                                                                          >  7:228300/1‑149 (MQ=255)
                       caAAAATGAGTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCGTAGATTGGGGCGTGGAGGAGCGGGTTTAccc                                                                                                                         >  2:198104/1‑149 (MQ=255)
                       caAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTAccc                                                                                                                         <  7:222882/149‑1 (MQ=255)
                                cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                    >  5:259537/1‑130 (MQ=255)
                                cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGc                                                                                                                                    <  6:259537/130‑1 (MQ=255)
                                 tCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTgg                                                                                                                                                 >  6:290035/1‑116 (MQ=255)
                                  cAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGc                                                                                                               >  3:64680/1‑148 (MQ=255)
                                   aaTGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCa                                                                                                             >  5:179614/1‑149 (MQ=255)
                                       aaCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATccc                                                                                                         >  5:297529/1‑149 (MQ=255)
                                        aCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAgtgt                                                                                                                                          <  2:174849/116‑1 (MQ=255)
                                        aCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAgtgt                                                                                                                                          >  1:174849/1‑116 (MQ=255)
                                        aCAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCg                                                                                                        <  6:297529/149‑1 (MQ=255)
                                         cacaaccggTCATTTTTAGAGTAATTATCCAAAATGCCGCACATACGGCCGTATACACCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTCGCAGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                       <  4:134286/142‑1 (MQ=38)
                                         caaaccAGGTCATTTTTACACTAATTATCCAAAATGACGCAGATAAGGCTGTATACAGCCTGTTAGAGCGGTTTAAGCCATCCGCGGGT‑‑AACCCCGCTCCTACATTGGCCAGTGTAGGCGCGGGTTTACCCGCGAAGAAGCCATCCCgg                                                                                                       <  4:262751/144‑1 (MQ=255)
                                            tACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTAt                                                                                                   >  5:106571/1‑149 (MQ=255)
                                                           cattAATTATCCAAAAGCTCGCAGATCAGGCTCTATAGAGCCTGTTAGACCGGTTAAAGGCATTCGCGGGT‑‑AAACCCGCGCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgcg                                                                                     <  8:228300/145‑1 (MQ=38)
                                                                aTTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTc                                                                               >  5:17271/1‑149 (MQ=255)
                                                                     ccAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCccgcc                                                                          <  3:345504/149‑1 (MQ=255)
                                                                       aaaaTGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAAACCCCAGCTACATTGGGGAGTGTAGGAGCGGGGTTTCCCGCGAAGAAGCCATCCCGGTTTAACGCGATATAGCGCAACTCCcgccgc                                                                        >  6:238981/1‑149 (MQ=255)
                                                                             tCGCAGATAAGGCTGTATAGAGCGTGTTAGATCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATc                                                                  <  3:273689/149‑1 (MQ=255)
                                                                               gCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCtt                                                                <  1:198104/149‑1 (MQ=255)
                                                                                cAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATTCACGCTATATCGCGCAACTCCCGCCGCAAAATCtt                                                                <  1:53106/149‑1 (MQ=255)
                                                                                cAGATAAGGCGGTATAGAGCCTGTTAGAGCGGTTAAAGGCATTCGCGGGT‑‑AAACCGGCTCATACATTGCGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTg                                                               <  7:133744/149‑1 (MQ=39)
                                                                                         cTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc                                                      <  6:17271/149‑1 (MQ=255)
                                                                                               agagCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                >  8:57090/1‑149 (MQ=255)
                                                                                                gagCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATg                                                <  6:106571/148‑1 (MQ=255)
                                                                                                 agCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATggg                                              >  5:362974/1‑149 (MQ=255)
                                                                                                   agtGTTAGAGAGGTTTTAGTGTTTCGCGGGTTAAAACTCCTTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCTAAATCTTGCCCACCGGCGTCATgg                                               <  1:300147/146‑1 (MQ=255)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAAt                                                                   >  6:284027/1‑123 (MQ=255)
                                                                                                      gTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAAt                                                                   <  5:284027/123‑1 (MQ=255)
                                                                                                           agCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcg                                    >  5:281569/1‑149 (MQ=255)
                                                                                                            gCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgc                                   <  7:175279/149‑1 (MQ=255)
                                                                                                             cGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgcg                                  <  2:147862/149‑1 (MQ=255)
                                                                                                               gTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgcg                                  >  7:257158/1‑147 (MQ=255)
                                                                                                                   aaGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgc                                   >  7:109380/1‑142 (MQ=255)
                                                                                                                   aaGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTcgc                                   <  8:109380/142‑1 (MQ=255)
                                                                                                                     gccATCCGCGGTA‑‑AAACCCGCTCCTACATTGGCAAGTGACGTAGGGGTTTAACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGa                          <  5:238981/147‑1 (MQ=255)
                                                                                                                             cgGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTg                  <  5:280239/149‑1 (MQ=255)
                                                                                                                              gggT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTgg                 <  2:320965/149‑1 (MQ=255)
                                                                                                                                 t‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTa              >  1:43782/1‑149 (MQ=255)
                                                                                                                                      aCCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACtt           <  8:257158/149‑1 (MQ=255)
                                                                                                                                       cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAATAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATgg                                               >  5:54831/1‑112 (MQ=37)
                                                                                                                                       cccGCTCCTACACTGCCCAATGTAGGAGCGGGTTTACCCGCGAATAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATgg                                               <  6:54831/112‑1 (MQ=37)
                                                                                                                                        ccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg         <  7:57090/149‑1 (MQ=255)
                                                                                                                                          gCTCCTACACTGCCCAATGTAGGAGCGGAGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGt        >  1:22715/1‑149 (MQ=39)
                                                                                                                                               atcATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCgg  <  5:327608/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCCAAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGT‑‑AAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCATGGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGG  >  NC_002947/5173099‑5173388

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