Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,464,945	A→T	7.7%	L58Q (CTG→CAG)	lipB ←	octanoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,464,945	0	A	T	7.7%	150.0 / 3.4	52	L58Q (CTG→CAG)	lipB	octanoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base T (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTC  >  NC_002947/5464798‑5465090                                                                                                                                                      |                                                                                                                                                     aGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACcagcag                                                                                                                                                     <  2:143378/149‑1 (MQ=255)  gCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAgg                                                                                                                                                    <  5:272517/149‑1 (MQ=255)             cTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTCACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCtt                                                                                                                                         <  7:243626/149‑1 (MQ=255)                aCCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGc                                                                                                                                       <  2:32540/148‑1 (MQ=255)                aCCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCCGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTg                                                                                                                                        <  4:34078/147‑1 (MQ=255)                  cAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccggg                                                                                                                                             >  7:86594/1‑136 (MQ=255)                  cAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccggg                                                                                                                                             <  8:86594/140‑5 (MQ=255)                     cTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCt                                                                                                                                 <  4:248771/149‑1 (MQ=255)                     cTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCt                                                                                                                                 <  2:105423/149‑1 (MQ=255)                      tCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAgg                                                                                                                                                    <  4:124616/129‑1 (MQ=255)                      tCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAgg                                                                                                                                                    >  3:124616/1‑129 (MQ=255)                      tCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTg                                                                                                                                >  8:255636/1‑149 (MQ=255)                      tCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTg                                                                                                                                >  2:22982/1‑149 (MQ=255)                            acgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑ccc                                                                                                                          >  4:342328/1‑146 (MQ=255)                            acgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑ccc                                                                                                                          <  3:342328/146‑1 (MQ=255)                             cgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTgg                                                                                                                      <  7:255636/149‑1 (MQ=255)                             cgcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTgg                                                                                                                      <  6:188373/149‑1 (MQ=255)                              gcCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTggg                                                                                                                     <  8:96202/149‑1 (MQ=255)                                      ccAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccg                                                                                                                                               >  6:375151/1‑116 (MQ=255)                                      ccAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCTGCTGGTGccg                                                                                                                                               <  5:375151/118‑3 (MQ=255)                                         gCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCa                                                                                                          <  2:276784/149‑1 (MQ=255)                                         gCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCa                                                                                                          <  8:136886/149‑1 (MQ=255)                                           cGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAgg                                                                                                        <  5:165013/149‑1 (MQ=255)                                             gcgcACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGt                                                                                                      <  3:263807/149‑1 (MQ=255)                                             gcgcACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGt                                                                                                      >  2:115454/1‑149 (MQ=255)                                                 acaTCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGCATGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCt                                                                                                   >  1:105827/1‑149 (MQ=255)                                                     ccAGCAGCAGGTAGGCCACCAGCTGCCCGGGGCTCATGGTAGGTGTCCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGAc                                                                                               >  1:45023/1‑149 (MQ=255)                                                               gTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGAttt                                                                                    >  6:230082/1‑149 (MQ=255)                                                                  ggCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCa                                                                                  <  1:184642/148‑1 (MQ=255)                                                                                       cATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCAt                                                                                 >  2:323578/1‑129 (MQ=255)                                                                                       cATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCAt                                                                                 <  1:323578/129‑1 (MQ=255)                                                                                       cATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGgctgcgc                                                             >  8:60231/1‑149 (MQ=255)                                                                                            tAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGGCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTc                                                        >  1:260584/1‑149 (MQ=255)                                                                                              ggTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCg                                                       <  2:105827/148‑1 (MQ=255)                                                                                               gTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGt                                                     <  2:260584/149‑1 (MQ=255)                                                                                                ggACCTGGCCGG‑CGGGGGCGGGCTGCACCACCGGAATGGCCCCCGGCACCACCATGTGCTTGGCCTTGCCGGCCTGCCCCCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCg                                                       <  1:115454/148‑1 (MQ=255)                                                                                                  aCCTGGCAGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCCCCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCAAGATTACATCCTGGCTGTCCGGGCTGCGCTGCTCGGGaa                                                   >  6:205578/1‑148 (MQ=255)                                                                                                    cTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAac                                                >  4:36511/1‑149 (MQ=255)                                                                                                        ccgc‑cgCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAacgacg                                            <  5:230082/149‑1 (MQ=255)                                                                                                         cgc‑cgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCt                                                                              <  6:188720/114‑1 (MQ=255)                                                                                                         cgc‑cgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCt                                                                              >  5:188720/1‑114 (MQ=255)                                                                                                          gc‑cgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCACAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCa                                          >  2:229719/1‑149 (MQ=255)                                                                                                           c‑cgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCAt                                         >  4:251022/1‑149 (MQ=255)                                                                                                           c‑cgcGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCACAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCAt                                         <  1:229719/149‑1 (MQ=255)                                                                                                                ggtcggtcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGcc                                     >  7:256473/1‑149 (MQ=255)                                                                                                                 gtcggtcTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCt                                    >  5:298218/1‑149 (MQ=255)                                                                                                                          caccacCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCg                           >  1:172906/1‑149 (MQ=255)                                                                                                                             caccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCa                                                                                  <  1:25815/91‑1 (MQ=255)                                                                                                                             caccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGtt                        >  4:230298/1‑149 (MQ=255)                                                                                                                             aaccGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCATGGTGCTCGACCAGAAAGATTTCa                                                                                  >  2:25815/2‑91 (MQ=255)                                                                                                                              accGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGtt                        >  1:47271/1‑148 (MQ=255)                                                                                                                                 ggAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATa                    <  7:60231/149‑1 (MQ=255)                                                                                                                                          cccGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAg           >  6:360225/1‑149 (MQ=255)                                                                                                                                                   agcagGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTc  >  1:36273/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGCTCGATACGGCTGACCAGCTCACGCACGCCAAAGCCCAGCCGGCGCACATCCAGCAGCAGGTAGGCCACCAGCTGCCCGGGGC‑CATGGTAGGTGACCTGGCCGC‑CGCGGTCGGTCTGCACCACCGGAATGTCGCCCGGCACCAGCAGGTGCTCGGCCTTGCCGGCCTG‑‑‑CCCCTGGGTGAAGACCGCAGGGTGCTCGACCAGCCAGATTTCATCCTGGCTGTCCGGGCTGCGCTGCTCGGTAAAACGACGCATGGCCTCCAGCACCGGTTCATAGGGCTGCAGGCCAAGCTC  >  NC_002947/5464798‑5465090

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