Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,482,041	T→C	8.0%	S152G (AGT→GGT)	estB ←	carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,482,041	0	T	C	8.0%	138.1 / 2.6	51	S152G (AGT→GGT)	estB	carboxylesterase
Reads supporting (aligned to +/- strand):  ref base T (26/20);  new base C (2/2);  total (29/22)
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.46e-01

CTTCGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCCGCGTCAGG  >  NC_002947/1481894‑1482183                                                                                                                                                     |                                                                                                                                                cTTCGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACt                                                                                                                                                 >  4:58925/1‑148 (MQ=255) cTTCGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGAACACCGAGTCATGCACGCCATGCAGGCACA‑GGGAAGGGGTGCGTTGCTGGCAGGAACt                                                                                                                                                 >  8:241628/1‑148 (MQ=255)  ttCGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTc                                                                                                                                                >  6:135283/1‑148 (MQ=255)    cGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCCAc                                                                                                                                             >  6:91427/1‑149 (MQ=255)     gtgtTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAAc                                                                                                                                             <  6:330465/148‑1 (MQ=255)            aTCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTg                                                                                                                                    >  2:177311/1‑149 (MQ=255)                 gTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCAt                                                                                                                               >  4:302034/1‑149 (MQ=255)                     tCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGAACGACATACAGGAACA‑GGGCCGGGGTGCGTTGATGGAAGGCAAACAACT‑GGGGGTGGTAATTg                                                                                                                             >  8:207273/1‑147 (MQ=255)                       gTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAg                                                                                                                          <  8:96599/148‑1 (MQ=255)                          ccAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGACGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATTCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTg                                                                                                                       <  7:241628/148‑1 (MQ=255)                                ggCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTc                                                                                                                                 <  3:136452/132‑1 (MQ=255)                                ggCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTc                                                                                                                                 >  4:136452/1‑132 (MQ=255)                                 gCAGCGACGCCCCAGGTGTTCTGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCa                                                                                                                 >  7:161199/1‑147 (MQ=255)                                  cAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAg                                                                                                              <  5:135283/149‑1 (MQ=255)                                     cGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTg                                                                                                           >  1:325405/1‑149 (MQ=255)                                        cGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGAc                                                                                                        <  2:126513/149‑1 (MQ=255)                                            ccAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAggg                                                                                                    >  6:42487/1‑149 (MQ=255)                                              aGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCg                                                                                                  <  4:99208/149‑1 (MQ=255)                                                gtgtTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACCCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGAt                                                                                                >  8:99658/1‑149 (MQ=255)                                                      aGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCAcccc                                                                                          >  6:329939/1‑149 (MQ=255)                                                      aGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑CGTGCTGGTCATTGAAGGTGGGGGCATATGTGGACAGGGCGATCAcccc                                                                                          >  5:44201/1‑149 (MQ=255)                                                        gTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGc                                                                                        >  8:297974/1‑149 (MQ=255)                                                           cTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGcc                                                                                       >  5:121303/1‑147 (MQ=255)                                                            tCGAAGGCGGTACGGCCCATGGTCGGGATTTCCACCTATTCATGGACGCCATTCAGGCACATGTGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGa                                                                                                                            <  1:74815/110‑1 (MQ=255)                                                            tCGAAGGCGGTACGGCCCATGGACGGGATGACAACCGAGTAATGAACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGa                                                                                                                            >  2:74815/1‑109 (MQ=255)                                                               aaGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGc                                                                                 >  7:13624/1‑149 (MQ=255)                                                                  gCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGACCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTc                                                                              >  8:163795/1‑149 (MQ=255)                                                                   cGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACTCCGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTc                                                                              <  1:267648/149‑1 (MQ=255)                                                                             cATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTta                                                                   <  3:302034/149‑1 (MQ=255)                                                                               tGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTtata                                                                 >  4:12855/1‑149 (MQ=255)                                                                                 gACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTtatat                                                                <  7:297974/148‑1 (MQ=255)                                                                                    gggATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACCCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAgg                                                             >  7:310727/1‑148 (MQ=255)                                                                                     ggATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCa                                                           <  1:177311/149‑1 (MQ=255)                                                                                     ggATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACCCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCa                                                           <  8:310727/149‑1 (MQ=255)                                                                                      gATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAg                                                          >  5:135588/1‑149 (MQ=255)                                                                                      gATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGGGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCA‑CCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAgt                                                         >  4:247609/1‑149 (MQ=255)                                                                                          accaccGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGtgtg                                                      >  6:271570/1‑149 (MQ=255)                                                                                           ccaccGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTAAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTgc                                                     <  7:207273/149‑1 (MQ=255)                                                                                             accGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGAAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTgcag                                                   >  3:96268/1‑149 (MQ=255)                                                                                               cGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTgcagca                                                 >  4:187754/1‑149 (MQ=255)                                                                                                 aGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGcacc                                               >  6:113209/1‑149 (MQ=255)                                                                                                   tCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGcaccac                                             <  5:329939/149‑1 (MQ=255)                                                                                                     aTGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGGCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGATTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACcg                                           >  5:162398/1‑149 (MQ=255)                                                                                                      tGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACcgc                                          <  3:12855/149‑1 (MQ=255)                                                                                                                 cAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTggg                               <  8:161199/149‑1 (MQ=255)                                                                                                                        a‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCA‑CCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAAAAGCCg                       <  3:247609/149‑1 (MQ=255)                                                                                                                          gggCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCg                       <  5:42487/149‑1 (MQ=255)                                                                                                                             ccGGGGTGCGTTGCTGGCAGGCACCCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGcc                    <  7:99658/149‑1 (MQ=255)                                                                                                                                 ggTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAgg                 >  8:254800/1‑148 (MQ=255)                                                                                                                                          cTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCcgc        <  7:163795/148‑1 (MQ=255)                                                                                                                                             gcaggcaCTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCCGCGTc     <  5:113209/148‑1 (MQ=255)                                                                                                                                               aggcaCTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCCGCGTCAgg  >  5:174510/1‑149 (MQ=255)                                                                                                                                               aggcaCTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCCGCGTCAgg  >  5:208627/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                CTTCGTGTTCCATCGGGTATTCGTACCAGCGGGCAGCGACGCCCCAGGTGTTCAGGTACTCGAAGGCGGTACGGCCCATGGACGGGATGACCACCGAGTCATGCACGCCATGCAGGCACA‑GGGCCGGGGTGCGTTGCTGGCAGGCACTCAACT‑GGTGCTGGTCATTGAAGGTGGGGGCATAGGTGGACAGGGCGATCACCCCGCCCAGTGCTTCCTGCCACTTTATATAGGCAGTGTGCAGCACCACCGCGCCGCCCTGGGAGAAGCCGGCCAGGAAGATCCGCGTCAGG  >  NC_002947/1481894‑1482183

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑