Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,868,018	A→G	10.9%	L182L (CTT→CTC)	lgt ←	prolipoprotein diacylglyceryl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,868,018	0	A	G	10.9%	140.5 / 7.8	55	L182L (CTT→CTC)	lgt	prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base G (4/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAGTTGCCGATACGCCCGGCCCCCA  >  NC_002947/5867875‑5868160                                                                                                                                                 |                                                                                                                                               ccGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTTGAGGGc                                                                                                                                           >  2:213530/1‑149 (MQ=255)  cGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGc                                                                                                                                           >  1:336702/1‑148 (MQ=255)  cGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCg                                                                                                                                          <  7:274275/149‑1 (MQ=255)      ttGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTTGAGGGCgtggc                                                                                                                                      <  1:213530/149‑5 (MQ=255)      ttGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGaa                                                                                                                                        >  8:236235/1‑147 (MQ=255)      ttGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGaa                                                                                                                                        <  7:236235/147‑1 (MQ=255)             tCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTTGAGGGCgtgg                                                                                                                                       >  2:214637/1‑138 (MQ=255)                   aCCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGa                                                                                                                         >  3:20271/1‑149 (MQ=255)                       ggACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGAt                                                                                                                        >  8:334271/1‑146 (MQ=255)                       ggACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGAt                                                                                                                        <  7:334271/146‑1 (MQ=255)                             aCTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGgcgc                                                                                                               <  2:336702/149‑1 (MQ=255)                              cTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGaa                                                                                                                                        <  1:106771/123‑1 (MQ=255)                              cTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGaa                                                                                                                                        >  2:106771/1‑123 (MQ=255)                              cTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGcgcg                                                                                                              >  5:357177/1‑149 (MQ=255)                                                    ccATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCCTACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTg                                                                                          >  1:104112/1‑148 (MQ=255)                                                           acaACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAAcggcgg                                                                                  >  3:269218/1‑148 (MQ=255)                                                             aaCGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCggg                                                                                 >  7:132922/1‑147 (MQ=255)                                                               cGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTa                                                                                                                                  <  2:232469/96‑1 (MQ=255)                                                               cGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTa                                                                                                                                  >  1:232469/1‑96 (MQ=255)                                                                gagaACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCt                                                                                           >  7:215295/1‑134 (MQ=255)                                                                gagaACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCt                                                                                           <  8:215295/134‑1 (MQ=255)                                                                gagaACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGaac                                                                              >  1:297473/1‑146 (MQ=255)                                                                  gaACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGatca                                                                          >  7:179744/1‑149 (MQ=255)                                                                          ccGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCAt                                                                  >  8:317638/1‑149 (MQ=255)                                                                             acacCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATggg                                                               <  4:269218/149‑1 (MQ=255)                                                                                   ccATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAAcggc                                                                                    >  3:245739/1‑122 (MQ=255)                                                                                   ccATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAAcggc                                                                                    <  4:245739/122‑1 (MQ=255)                                                                                   ccATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCg                                                         >  1:144462/1‑149 (MQ=255)                                                                                      tGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTg                                                      >  7:137489/1‑149 (MQ=255)                                                                                        gTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTgg                                                     <  6:7576/148‑1 (MQ=255)                                                                                         tCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGt                                                   <  2:297473/149‑1 (MQ=255)                                                                                             ccGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGc                                               <  7:317638/149‑1 (MQ=255)                                                                                               gTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGccc                                             <  4:20271/149‑1 (MQ=255)                                                                                                  gCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCcac                                          <  8:137489/149‑1 (MQ=255)                                                                                                  gCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCca                                           <  8:132922/148‑1 (MQ=255)                                                                                                    ttGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATggg                                                               >  5:335956/1‑126 (MQ=255)                                                                                                    ttGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATggg                                                               <  6:335956/126‑1 (MQ=255)                                                                                                    ttGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAg                                        >  4:97788/1‑149 (MQ=255)                                                                                                        gcgAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGtt                                      <  4:1018/147‑1 (MQ=255)                                                                                                        gcgAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGtt                                      >  3:1018/1‑147 (MQ=255)                                                                                                             aaCAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGCGCGAACTGGAAGACCTCCGATTGGTGGCCCGGCAACTCGGCCCGGTCGCAGAACGCCCGGAAGATCATCGCCCATGCGACGACGCAGGGATTGcac                                             >  5:116472/1‑133 (MQ=255)                                                                                                              aCAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTtg                              <  2:104112/149‑1 (MQ=255)                                                                                                                     acaGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAg                        <  6:357177/149‑1 (MQ=255)                                                                                                                         gATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAGTTGc                    >  4:55489/1‑149 (MQ=255)                                                                                                                             aCGAACATCTCCACGCCCTAAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCTGCGTTGATGAAGTTGCCGa                 <  8:179744/148‑1 (MQ=255)                                                                                                                                   accagttcGCCCTTGAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAGTTGCCGATACGccc          >  5:48557/8‑149 (MQ=255)                                                                                                                                       gttcGCCCTTGAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGc                                               >  7:144662/4‑108 (MQ=255)                                                                                                                                       gttcGCCCTTGAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGc                                               <  8:144662/105‑1 (MQ=255)                                                                                                                                           gCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAGTTGCCGATACGCCCGGCCCCCa  <  3:55489/149‑1 (MQ=255)                                                                                                                                             ccTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCcac                                          >  8:194144/1‑107 (MQ=255)                                                                                                                                             ccTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCcac                                          <  7:194144/107‑1 (MQ=255)                                                                                                                                             ccTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCca                                           >  1:222416/1‑106 (MQ=255)                                                                                                                                             ccTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCca                                           <  2:222416/106‑1 (MQ=255)                                                                                                                                              cTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGtt                               <  6:172982/117‑1 (MQ=255)                                                                                                                                              cTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGtt                               >  5:172982/1‑117 (MQ=255)                                                                                                                                                 |                                                                                                                                               CCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGCCATAGCACAACGAGAACATTCCCGACACCGCCATGGTCGGCCGTGGCTTGCGCGAGAACAGCC‑ACAGGATGACGAACAGCGCCACGCCCTCAAGGGCGAACTGGTAGAGCTGCGATGGGTGGCGCGGCAACTGGGCCGGGTCGCTGAACGGCGGGAAGATCATCGCCCATGGGACGTCGGTGGGTTTGCCCCACAGTTCGGCGTTGATGAAGTTGCCGATACGCCCGGCCCCCA  >  NC_002947/5867875‑5868160

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