Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,285,091	C→T	7.5%	P239L (CCA→CTA)	yejM →	essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,285,091	0	C	T	7.5%	166.5 / 4.1	53	P239L (CCA→CTA)	yejM	essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (2/2);  total (27/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 = 5.47e-01

GCCAACTTCTATCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAA  >  NC_000913/2284952‑2285239                                                                                                                                            |                                                                                                                                                     gCCAACTTCTATCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGa                                                                                                                                                   <  2:463963/143‑1 (MQ=255) gCCAACTTCTATCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGa                                                                                                                                                   >  1:463963/1‑143 (MQ=255)   cAACTTCTATCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGttt                                                                                                                                           <  1:18721/149‑1 (MQ=255)            tCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGt                                                                                                                                  <  2:444259/149‑1 (MQ=255)                   aTCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGtt                                                                                                                           >  2:323710/1‑149 (MQ=255)                      aCCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAg                                                                                                                        >  2:16627/1‑149 (MQ=255)                        cATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCg                                                                                                                      >  2:239846/1‑149 (MQ=255)                        cATGCAGCACGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCg                                                                                                                      >  4:41326/1‑149 (MQ=255)                                      aCCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATcgc                                                                                                        >  3:170662/1‑149 (MQ=255)                                      aCCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATcgc                                                                                                        <  1:239846/149‑1 (MQ=255)                                       ccTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATcgcg                                                                                                       >  4:265072/1‑149 (MQ=255)                                        cTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGa                                                                                                      >  3:33067/1‑149 (MQ=255)                                        cTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGa                                                                                                      >  2:277692/1‑149 (MQ=255)                                        cTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGa                                                                                                      >  2:219064/1‑149 (MQ=255)                                          gCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGtt                                                                                                                                      <  1:1161/115‑1 (MQ=255)                                          gCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGtt                                                                                                                                      >  2:1161/1‑115 (MQ=255)                                                                   cgACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGtt                                                                           <  2:215328/149‑1 (MQ=255)                                                                    gACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCTAGACGCCGTTTCCGTTCAGTATCCGt                                                                                                                            <  3:452359/99‑1 (MQ=255)                                                                    gACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCTAGACGCCGTTTCCGTTCAGTATCCGt                                                                                                                            >  4:268523/1‑99 (MQ=255)                                                                    gACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCTAGACGCCGTTTCCGTTCAGTATCCGt                                                                                                                            >  4:452359/1‑99 (MQ=255)                                                                    gACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCTAGACGCCGTTTCCGTTCAGTATCCGt                                                                                                                            <  3:268523/99‑1 (MQ=255)                                                                     tcGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCTTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGa                                                                         <  3:41326/148‑1 (MQ=255)                                                                     aCGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGa                                                                         <  4:125075/149‑1 (MQ=255)                                                                      cGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGAt                                                                        <  3:265072/149‑1 (MQ=255)                                                                       gTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTgc                                                                                                                >  3:5373/1‑108 (MQ=255)                                                                       gTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTgc                                                                                                                <  4:5373/108‑1 (MQ=255)                                                                       gTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGAtt                                                                       <  2:233553/149‑1 (MQ=255)                                                                       gTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGAtt                                                                       <  4:96010/149‑1 (MQ=255)                                                                        ttttCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAAt                                                                                   >  2:318468/1‑136 (MQ=255)                                                                        ttttCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAAt                                                                                   <  1:318468/136‑1 (MQ=255)                                                                           tCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTcc                                                                                                                                         >  1:16647/1‑79 (MQ=255)                                                                           tCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTcc                                                                                                                                         <  2:16647/79‑1 (MQ=255)                                                                           tCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTg                                                                   >  1:322444/1‑149 (MQ=255)                                                                             ttGAGAAGCATGGTCTGCTTGATGCGCAGGGGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTc                                                                                                                                     >  4:143817/1‑81 (MQ=255)                                                                             ttGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTc                                                                                                                                     <  1:455490/81‑1 (MQ=255)                                                                             ttGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTc                                                                                                                                     <  3:143817/81‑1 (MQ=255)                                                                             ttGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTc                                                                                                                                     >  2:455490/1‑81 (MQ=255)                                                                               gagaAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGa                                                               >  2:195063/1‑149 (MQ=255)                                                                                    gCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGcc                                                          <  2:189588/149‑1 (MQ=255)                                                                                        ggTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGaa                                                      <  4:146067/149‑1 (MQ=255)                                                                                           cTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAacta                                                   <  2:322444/149‑1 (MQ=255)                                                                                            tGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAactac                                                  >  1:168638/1‑149 (MQ=255)                                                                                              cTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTc                                                <  2:168638/149‑1 (MQ=255)                                                                                              cTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTc                                                >  4:32777/1‑149 (MQ=255)                                                                                               ttGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCa                                               >  4:359833/1‑149 (MQ=255)                                                                                                          gAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGaa                                    <  1:219064/149‑1 (MQ=255)                                                                                                               tCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGa                               >  1:406794/1‑149 (MQ=255)                                                                                                                      cGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTgc                        >  4:150754/1‑149 (MQ=255)                                                                                                                        tCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTgcgc                      <  4:33067/149‑1 (MQ=255)                                                                                                                                aGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGaa                                                      <  4:246124/109‑1 (MQ=255)                                                                                                                                aGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGaa                                                      >  3:246124/1‑109 (MQ=255)                                                                                                                                  cAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGa                                                                         <  4:202839/88‑1 (MQ=255)                                                                                                                                  cAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGa                                                                         >  3:202839/1‑88 (MQ=255)                                                                                                                                            cAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCaa  >  1:243682/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     GCCAACTTCTATCGCCCGATCACCATGCAGCGCGCTAACCTGCCGCTTTCGTACCCGATGACGGCGCGACGTTTTCTTGAGAAGCATGGTCTGCTTGATGCGCAGGAGTATCAACGCCGTCTTATTGAGCAAGGTAATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAA  >  NC_000913/2284952‑2285239

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