Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,136,104	A→T	10.7%	L216Q (CTG→CAG)	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,136,104	0	A	T	10.7%	123.5 / 4.3	45	L216Q (CTG→CAG)	fadE	medium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (2/3);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTTGCAGATGATGCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGGATGATCAGGGCGAAAAA  >  NC_002947/2135959‑2136252                                                                                                                                                  |                                                                                                                                                     gCCTTTGCAGATGATGCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTgcagcagc                                                                                                                                                    <  2:328481/148‑1 (MQ=255)      tGCAGATGATGCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTc                                                                                                                                               <  6:108831/148‑1 (MQ=255)                gCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGcc                                                                                                                                      >  1:41996/1‑147 (MQ=255)                gCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGcc                                                                                                                                      >  2:332462/1‑147 (MQ=255)                   ggTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGTGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAggg                                                                                                                                 <  3:318985/149‑1 (MQ=255)                    gtgtCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAgg                                                                                                                                  >  8:257283/1‑147 (MQ=255)                    gtgtCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAgg                                                                                                                                  >  2:316005/1‑147 (MQ=255)                    gtgtCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGa                                                                                                                                <  7:257283/149‑1 (MQ=255)                       tCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGtt                                                                                                                             >  6:147986/1‑149 (MQ=255)                       tCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGtt                                                                                                                             <  8:101303/149‑1 (MQ=255)                       tCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGtt                                                                                                                             >  2:99878/1‑149 (MQ=255)                         gggCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAgg                                                                                                                                  >  4:85641/1‑142 (MQ=255)                         gggCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAgg                                                                                                                                  <  3:85641/142‑1 (MQ=255)                          ggCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGAACTGCTGCTGCACTACGGCac                                                                                                                                          <  6:111117/133‑2 (MQ=255)                          ggCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGAACTGCTGCTGCACTACGGCac                                                                                                                                          >  5:111117/1‑132 (MQ=255)                                  gcCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAg                                                                                                                                   >  1:298277/1‑132 (MQ=255)                                  gcCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAg                                                                                                                                   <  2:298277/132‑1 (MQ=255)                                    cAGCGACCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGCGATTTCCTCGCCACGGGCCAGCCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCAcc                                                                                                                >  5:234869/1‑149 (MQ=255)                                            gAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGc                                                                                                        <  8:221167/149‑1 (MQ=255)                                              aCCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCa                                                                                                      <  7:193215/149‑1 (MQ=255)                                              aCCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCa                                                                                                      <  2:291381/149‑1 (MQ=255)                                                cGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCTGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAgg                                                                                                    >  5:272223/1‑149 (MQ=255)                                                  gCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTc                                                                                                  <  1:332462/149‑1 (MQ=255)                                                    cAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTcgc                                                                                                >  5:83172/1‑149 (MQ=255)                                                            cGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCgg                                                                                        <  2:41996/149‑1 (MQ=255)                                                                   ggCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGc                                                                                   >  8:120716/1‑147 (MQ=255)                                                                     cGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGc                                                                               >  4:86873/1‑149 (MQ=255)                                                                           aGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATc                                                                         <  6:170779/149‑1 (MQ=255)                                                                           aGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATc                                                                         <  8:83356/149‑1 (MQ=255)                                                                           aGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATc                                                                         >  1:69500/1‑149 (MQ=255)                                                                                     cTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCTGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTggg                                                               <  6:272223/149‑1 (MQ=255)                                                                                      tCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTgg                                                                >  3:311295/1‑147 (MQ=255)                                                                                      tCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCAGGGTGGCCAGCTTCATCGCCACCTgg                                                                >  6:316540/1‑147 (MQ=255)                                                                                          cACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCCGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTgg                                                          >  8:258237/1‑149 (MQ=255)                                                                                               gCCATGGGCGTCAGGTTGTGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTa                                                     <  6:234869/149‑1 (MQ=255)                                                                                                    gcgcGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCg                                                 >  1:215028/1‑148 (MQ=255)                                                                                                                   gTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTgg                                                                                                                           >  6:345194/1‑59 (MQ=255)                                                                                                                   gTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTgg                                                                                                                           <  5:345194/59‑1 (MQ=255)                                                                                                                      gCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGccgccg                               <  8:145612/148‑1 (MQ=255)                                                                                                                        gttgttCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTAt                            >  4:316241/1‑149 (MQ=255)                                                                                                                               gtcgGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGTGAGTGGGCGTAGGCGGAGCAGCCCTTGCCGCCGTATTCCTTgg                     <  5:147986/149‑1 (MQ=255)                                                                                                                               gtcgGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTgg                     <  2:215028/149‑1 (MQ=255)                                                                                                                               gtcgGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTg                      <  7:120716/148‑1 (MQ=255)                                                                                                                                 cgGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTggg                    >  2:60913/1‑148 (MQ=255)                                                                                                                                      ccGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGgatgat              <  2:69500/149‑1 (MQ=255)                                                                                                                                           gTGCAGCTGCAGTTCGGCAGGGCCCAGGGAGTTGGTCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACGTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGGATGATCAgg          <  5:226106/148‑1 (MQ=255)                                                                                                                                                  agcagTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGGATGATCAGGGCGaaaaa  <  3:110902/149‑1 (MQ=255)                                                                                                                                                  agcagTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGGATGATCAGGGCGaaaaa  <  3:115673/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCCTTTGCAGATGATGCCGGTGTCGGGCATGGCGCCAGCGTCCGAACCGGCCAGCGGCCCGGTGAGGGCGAAACAGGGGATTTCCTCGCCACGGGCCAGGCGCGGCAGGTAGCGGTTGCGTTGTTCGTCGGTGCCGTAGTGCAGCAGCAGTTCGGCAGGGCCCAGGGAGTTGGGCACCATCACCGTGGAGGCCAGGTCGCCGCTCCGGGTGGCCAGCTTCATCGCCACCTGGGAGTGGGCGTAGGCGGAGAAGCCCTTGCCGCCGTATTCCTTGGGGATGATCAGGGCGAAAAA  >  NC_002947/2135959‑2136252

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: ‑