Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,134,916	Δ1 bp	7.7%	coding (1835/2448 nt)	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,916	0	T	.	7.7%	157.5 / 7.4	52	coding (1835/2448 nt)	fadE	medium‑long chain acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base . (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 = 1.00e+00

TCCGTCTGGCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCG  >  NC_002947/2134772‑2135062                                                                                                                                                 |                                                                                                                                                    tCCGTCTGGCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTcgcg                                                                                                                                                 >  7:101930/1‑149 (MQ=255)         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                                                                                                                                  <  2:226424/120‑1 (MQ=38)                                            gCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTCAAGCGC‑CGCGAGCGCCTCAGCGccc                                                                                                                                  >  1:226424/1‑120 (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                                           >  5:93820/1‑129 (MQ=255)                                                                                                                          cGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCgggg                                           <  6:93820/129‑1 (MQ=255)                                                                                                                            ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAg                                                                     >  8:66741/1‑101 (MQ=255)                                                                                                                            ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAg                                                                     <  7:66741/101‑1 (MQ=255)                                                                                                                            ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAgg                                            >  2:247124/1‑126 (MQ=255)                                                                                                                            ccAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAgg                                            <  1:247124/126‑1 (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)                                                                                                                                                 |                                                                                                                                                    TCCGTCTGGCCCAGGCTTTCTTCCATGGCCCAGCGCAGCAGCGGCTGCAGGTATTCGGGCGAGCCCTGGTCGTGATAGCGCTTCAGGGCAGCGCTGGCGAGGTACAGATAGCTCAGCACATCGCCAAGCCGGGCGCTGAGGCGCTCGCG‑GCGCTTGAGCGCCCCACCCAGCAGCATCATCGAAAGGTCGGCCATCATGGCAAACGCCGCCGCCTGGCGGTTCAGGGCGCGGAAGTAGCCTTGGCTCAGGGCATCGCCCGGCACCCGCTCAAAGCGACCCAGCCCCAGGCCG  >  NC_002947/2134772‑2135062

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