Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,291,564	T→A	7.3%	L119M (TTG→ATG) ‡	estC →	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,291,564	0	T	A	7.3%	160.1 / 3.3	55	L119M (TTG→ATG) ‡	estC	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

TTCGCGGCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGCGAG  >  NC_002947/1291423‑1291710                                                                                                                                              |                                                                                                                                                   ttCGCGGCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGc                                                                                                                                             >  7:68431/1‑149 (MQ=255)    gcgGCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTTGACCTCGATGCCCCGGAGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGACGATCACCTTGCGTCAGCTGCTCCGCCACCCTGCCGGCTTGCCGGCCAt                                                                                                                                          >  6:213701/1‑149 (MQ=255)    gcgGCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCa                                                                                                                                           <  1:42364/148‑1 (MQ=255)       gCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCg                                                                                                                                       >  2:163875/1‑149 (MQ=255)        ccGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGcacggtgg                                                                                                                                      <  8:68431/149‑8 (MQ=255)               cGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGc                                                                                                                                                 >  6:167192/1‑131 (MQ=255)               cGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGc                                                                                                                                                 <  5:167192/131‑1 (MQ=255)                  ccTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCtgctgc                                                                                                                            >  2:47615/1‑149 (MQ=255)                   cTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGcc                                                                                                                           >  2:256056/1‑149 (MQ=255)                             tCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGcc                                                                                                                 >  3:94589/1‑149 (MQ=255)                               ggCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt                                                                                                               <  5:217457/149‑1 (MQ=255)                               atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACTGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt                                                                                                               >  6:194516/3‑149 (MQ=255)                               atcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCtt                                                                                                               >  8:198031/3‑149 (MQ=255)                                tcGAGGGCCAGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTg                                                                                                              >  4:180255/2‑149 (MQ=255)                                 cGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGt                                                                                                             <  4:269628/149‑1 (MQ=255)                                  gAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGcacggt                                                                                                                                        <  3:276716/121‑6 (MQ=255)                                  gAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCAAGCCGGcacggt                                                                                                                                        >  4:276716/1‑116 (MQ=255)                                       ccaGCTTGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACt                                                                                                       <  3:180255/147‑1 (MQ=255)                                               cccTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGAt                                                                                                >  8:229325/1‑148 (MQ=255)                                               cccTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATc                                                                                               <  6:217108/149‑1 (MQ=255)                                                 cTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCAt                                                                                             <  7:198031/149‑1 (MQ=255)                                                    gATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGt                                                                                          <  1:163875/149‑1 (MQ=255)                                                     aTGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGt                                                                                                                                      <  8:41307/104‑1 (MQ=255)                                                     aTGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGt                                                                                                                                      >  7:41307/1‑104 (MQ=255)                                                        ccccGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTc                                                                                                                                        <  3:111367/99‑1 (MQ=255)                                                        ccccGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTc                                                                                                                                        >  4:111367/1‑99 (MQ=255)                                                        ccccGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGAt                                                                                      <  2:219705/149‑1 (MQ=255)                                                         cccGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATg                                                                                     >  2:218542/1‑149 (MQ=255)                                                         cccGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATg                                                                                     >  2:327485/1‑149 (MQ=255)                                                             gTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACTGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCt                                                                                 <  5:194516/149‑1 (MQ=255)                                                                gCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTgg                                                                               <  1:47615/148‑1 (MQ=255)                                                                   aaCTACTGGCCCGAGTTCGACCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGc                                                                           >  2:18826/1‑149 (MQ=255)                                                                    actactGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCTGATCATGGTCGATGCCCTGGCAGcc                                                                          <  3:115185/149‑1 (MQ=255)                                                                        ctGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCgaga                                                                      >  5:152902/1‑149 (MQ=255)                                                                        ctGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCga                                                                        <  7:229325/147‑1 (MQ=255)                                                                           gCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGAc                                                                     <  1:327485/147‑1 (MQ=255)                                                                             ccGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGcc                                                                                                                 >  3:287556/1‑101 (MQ=255)                                                                             ccGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGcc                                                                                                                 <  4:287556/101‑1 (MQ=255)                                                                             ccGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACATTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCg                                                                 >  1:46357/1‑149 (MQ=255)                                                                                      cccAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACg                                                        >  7:42518/1‑149 (MQ=255)                                                                                          ggCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAg                                                    <  8:42518/149‑1 (MQ=255)                                                                                                   gCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTcc                                                                  <  1:253828/126‑1 (MQ=255)                                                                                                   gCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTcc                                                                  >  2:253828/1‑126 (MQ=255)                                                                                                               cTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTAcgccg                               >  1:173200/1‑149 (MQ=255)                                                                                                                            ctCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACgccgcc                             <  6:181372/138‑1 (MQ=255)                                                                                                                            ctCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACgccgcc                             >  5:181372/1‑138 (MQ=255)                                                                                                                            ctCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTtg                   <  6:152902/148‑1 (MQ=255)                                                                                                                               aGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTtggct                <  1:256056/148‑1 (MQ=255)                                                                                                                                gCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTtggctgg              >  5:353765/1‑149 (MQ=255)                                                                                                                                      gTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATc        >  2:120770/1‑149 (MQ=255)                                                                                                                                         ccGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGc     <  1:218542/149‑1 (MQ=255)                                                                                                                                         ccGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGc     <  2:173200/149‑1 (MQ=255)                                                                                                                                         ccGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGc     <  1:251432/149‑1 (MQ=255)                                                                                                                                           ggCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGCGa   <  1:120770/149‑1 (MQ=255)                                                                                                                                            gCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGCGAg  >  6:52134/1‑149 (MQ=255)                                                                                                                                              ttGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTAcgcc                                <  6:158924/117‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   TTCGCGGCCGTTACCGCCCTGCAGTTGGTCGGCGAGGGCAAGCTGGCCCTCGATGCCCCGGTGGCCAACTACTGGCCCGAGTTCGCCCAGGCGGGTAAGCAGTCGATCACCTTGCGTCAGCTGCTCAGCCACCGTGCCGGCTTGCCGGCCATTCGTGAGCTGCTGCCGGCTGAAGCCTTGTATGACTGGCAGATCATGGTCGATGCCCTGGCAGCCGAGACTCCGTGGTGGACGCCAGGCACAGAACACGGCTACGCCGCCATCACCTTTGGCTGGCTGATCGGCGAG  >  NC_002947/1291423‑1291710

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