Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,134,925	T→C	7.4%	K609R (AAG→AGG)	fadE ←	medium‑long chain acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,134,925	0	T	C	7.4%	151.1 / 2.7	54	K609R (AAG→AGG)	fadE	medium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (2/2);  total (28/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 = 7.44e-01

GCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAG  >  NC_002947/2134780‑2135071                                                                                                                                                   |                                                                                                                                                   gCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGa                                                                                                                                                 <  4:279024/149‑1 (MQ=255) gCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGa                                                                                                                                                 <  1:38799/149‑1 (MQ=255)    cAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAgcg                                                                                                                                              >  2:305125/1‑149 (MQ=255)       gatttCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGccc                                                                                                                                           >  8:204297/3‑149 (MQ=255)       gCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGccc                                                                                                                                           <  6:84553/149‑1 (MQ=255)          ttcttcCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCcccac                                                                                                                                        >  8:262466/1‑149 (MQ=255)             ttcCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTCAAGCGC‑CGCGAGCGCCTCAGCGcccggcttg                                                                                                                                     >  4:69927/1‑143 (MQ=38)              tcCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCcag                                                                                                                                    >  3:104368/1‑149 (MQ=255)               ccATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTCAAGCGC‑CGCGAGCGCCTCAGCGcccggcttg                                                                                                                                     <  3:69927/147‑7 (MQ=38)                      ccAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGGACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGcatc                                                                                                                             >  6:16355/1‑148 (MQ=255)                       cAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAgcgc                                                                                                                                             >  4:146391/1‑131 (MQ=255)                       cAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAgcgc                                                                                                                                             <  3:146391/131‑1 (MQ=255)                                  cGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCgg                                                                                                                <  5:157216/149‑1 (MQ=255)                                   ggCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGc                                                                                                               >  1:250669/1‑149 (MQ=255)                                    gCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTCAAGCGC‑CGCGAGCGCCTCAGCGccc                                                                                                                                           >  1:226424/1‑120 (MQ=38)                                    gCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTCAAGCGC‑CGCGAGCGCCTCAGCGccc                                                                                                                                           <  2:226424/120‑1 (MQ=38)                                          ggTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATg                                                                                                        <  1:305125/149‑1 (MQ=255)                                              ttCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGATCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCa                                                                                                     >  5:258367/1‑148 (MQ=255)                                                     gAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGc                                                                                                               <  2:236877/131‑1 (MQ=255)                                                     gAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGc                                                                                                               >  1:236877/1‑131 (MQ=255)                                                          cTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCt                                                                                        >  4:29905/1‑149 (MQ=255)                                                             gTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTgg                                                                                      <  2:156459/148‑1 (MQ=255)                                                             gTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGc                                                                                     >  3:154290/1‑149 (MQ=255)                                                              tCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCg                                                                                    <  7:197759/149‑1 (MQ=255)                                                                   aTAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCa                                                                               <  4:279354/149‑1 (MQ=255)                                                                      gcgcTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAgg                                                                             <  7:204297/148‑1 (MQ=255)                                                                       cgcTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATAGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCGTGCGGTTCAGGgc                                                                           >  1:138692/1‑149 (MQ=255)                                                                              gggCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAg                                                                    <  4:154290/149‑1 (MQ=255)                                                                                     gcTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCtt                                                             >  8:177097/1‑149 (MQ=255)                                                                                      cTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTg                                                            <  3:29905/149‑1 (MQ=255)                                                                                       tGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTgg                                                           >  1:44356/1‑149 (MQ=255)                                                                                         gCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTgg                                                           <  2:44356/147‑1 (MQ=255)                                                                                         gCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTgg                                                           <  7:262466/147‑1 (MQ=255)                                                                                           gAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCt                                                         >  2:123706/1‑147 (MQ=255)                                                                                            aGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAg                                                      <  8:101930/149‑1 (MQ=255)                                                                                                aCAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTgg                                                                                      <  8:24034/113‑1 (MQ=255)                                                                                                aCAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTgg                                                                                      >  7:24034/1‑113 (MQ=255)                                                                                                     tAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGcc                                             <  2:250669/149‑1 (MQ=255)                                                                                                       gCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGGCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCg                                           <  5:16355/149‑1 (MQ=255)                                                                                                             cacaTCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCAccc                                     >  4:68379/1‑149 (MQ=255)                                                                                                                 tCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCt                                  >  2:101524/1‑148 (MQ=255)                                                                                                                  cGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCgggg                                                    <  6:93820/129‑1 (MQ=255)                                                                                                                  cGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCgggg                                                    >  5:93820/1‑129 (MQ=255)                                                                                                                    ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAg                                                                              >  8:66741/1‑101 (MQ=255)                                                                                                                    ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAg                                                                              <  7:66741/101‑1 (MQ=255)                                                                                                                    ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAgg                                                     <  1:247124/126‑1 (MQ=255)                                                                                                                    ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAgg                                                     >  2:247124/1‑126 (MQ=255)                                                                                                                    ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCa                                >  6:248891/1‑147 (MQ=255)                                                                                                                         ccGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGAc                         >  8:179804/1‑149 (MQ=255)                                                                                                                                   aGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGcccc                 <  6:258367/147‑1 (MQ=255)                                                                                                                                    ggCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAgg              >  1:179875/1‑149 (MQ=255)                                                                                                                                         tCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCg           >  5:154074/1‑147 (MQ=255)                                                                                                                                           gcg‑gcgCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGa       >  1:356398/1‑149 (MQ=255)                                                                                                                                                 gcTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAg  <  1:123706/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCGTAGACCAAG  >  NC_002947/2134780‑2135071

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