Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,272,345	T→G	14.7%	K157N (AAA→AAC)	U712_16830 ←	Carboxylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,272,345	0	T	G	14.7%	117.2 / 2.3	56	K157N (AAA→AAC)	U712_16830	Carboxylesterase
Reads supporting (aligned to +/- strand):  ref base T (13/34);  new base G (7/2);  total (20/36)
Fisher's exact test for biased strand distribution p-value = 7.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAATATTGGCGCTTTCGGTATTAATCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTGCGCACATTGGGA  >  CP006881/3272207‑3272481                                                                                                                                           |                                                                                                                                          aataatATTGGCGCTTTAGGTATTAATCATGTTGTCATGAAGGGCCTGCACCAAAAATGTCTGTGAATAAATCATATCGACATTATTCCGCACATCAGAAATTAAGTCTTTCATCGCTTTTAGGGTATTCATCTGCGtt                                                                                                                                           >  8:344286/1‑139 (MQ=255)          ggCGCTTTCGGTATTAATCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAAt                                                                                                                                  >  6:155996/1‑139 (MQ=255)           gcgcTTTCGGTATTAATCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAAt                                                                                                                                  >  8:56999/1‑138 (MQ=255)                     tagtAATCCTGGGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTc                                                                                                                       <  1:334465/136‑1 (MQ=255)                     tATTAATCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTc                                                                                                                       <  4:252075/139‑1 (MQ=255)                         aaTCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTGGAATTCTTTCATTTccacg                                                                                                                   >  5:314479/1‑136 (MQ=255)                                   tCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCt                                                                                                         <  7:56999/139‑1 (MQ=255)                                     aTGACGGTCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTcc                                                                                                       <  4:120968/139‑1 (MQ=255)                                            gCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCttt                                                                                                <  2:271550/139‑1 (MQ=255)                                             ccTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTGGAATTGTTTTATTTAGTGTTGAAATTGGTGCGGGCtttt                                                                                               >  3:65735/1‑139 (MQ=255)                                             ccTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGAAGCGCCTTGAGGGGATTCATCGGCGTGTTTTCGAATTTTTTCATTTCCGCGTCAATTTTCTCCGGGCtttt                                                                                               >  8:285138/1‑139 (MQ=255)                                                   ccacAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTTTTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTc                                                                                         <  7:216811/139‑1 (MQ=255)                                                   ccacAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTc                                                                                         <  3:297880/139‑1 (MQ=255)                                                      caAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTGATTTACTCTTAAATTTGCTCCGGGCTTTTGCCCTCaaa                                                                                      >  7:299933/1‑139 (MQ=255)                                                      caAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTCGAAATAGGTGATTTCATATTTAAATTGCTACGTGCTTTTACTCTCaa                                                                                       >  4:355408/1‑138 (MQ=255)                                                      caAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTAGAATTCTTTGATTTCCTGATAAATTTTGTCCGGGGGTTTCCCCTCaaa                                                                                      >  8:152913/1‑139 (MQ=255)                                                           gTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTGGAATTCTTTCATTTGCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACtttt                                                                                 >  6:360032/1‑139 (MQ=255)                                                              ggTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTCGAAATGTTTCATTTACTCTTGAAATTCGAACGCGCTTTTCCCCTCAAACTTTTTg                                                                               >  8:367154/1‑138 (MQ=255)                                                               gTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCTGCGCCTTGAGGGTATTCATCGGCGTTTTTTAGAATTATTTCATTTGATCTTGAAATTTCTCCGGGCGTTTGCCCTCAAACttttttt                                                                              >  4:62488/1‑136 (MQ=255)                                                                      aaTCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCTGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGa                                                                     <  3:62488/139‑1 (MQ=255)                                                                      aaTCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGa                                                                     <  5:326598/139‑1 (MQ=255)                                                                          atatCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCg                                                                  <  3:5580/138‑1 (MQ=255)                                                                           tatCGACATTATTCCGCACATCAGCAATTAATTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTa                                                                <  7:356519/139‑1 (MQ=255)                                                                            atCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCtt                                                                                                                              <  4:213370/77‑1 (MQ=255)                                                                            atCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCtt                                                                                                                              >  3:213370/1‑77 (MQ=255)                                                                                 cATTATTACGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGGTTTTTCGAATTCTTTAATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAg                                                          <  5:236936/139‑1 (MQ=255)                                                                                    tattCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTAGAATTATTTGATTTACTCTTGAATTTTATACCGGCTTTTGACCTCAAACATTTTTTTA‑TTTAGATCGTTTGAAAGAAc                                                       >  2:320349/1‑139 (MQ=255)                                                                                     attCCGCACATCAGCAATTAAGTCTTGCAGCGCCGTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAAc                                                       <  7:367154/138‑1 (MQ=255)                                                                                         cGCACATGAGCAATTAAGTCTTTCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCAATTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCtt                                                  <  2:209819/139‑1 (MQ=255)                                                                                         cGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCtt                                                  <  4:256103/139‑1 (MQ=255)                                                                                           cacaTCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTGGAATTGTTTGATTTCCTCTTCAATTTGGTCCGGGGGTTTGCCCTCAAACTTTTTGTTA‑TTGGGCGCGTATGAAAGAACGCCTTGa                                                >  2:38196/1‑139 (MQ=255)                                                                                           cacaTCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGa                                                <  8:187658/139‑1 (MQ=255)                                                                                            acaTCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCGTTTCCCCTCAAACTTTTTGTTA‑TTGCGAGCGTATGAAAGAACGCCGTGa                                                >  2:5742/1‑138 (MQ=255)                                                                                                 aGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATg                                          <  4:19514/139‑1 (MQ=255)                                                                                                   cAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGAc                                        <  2:206847/139‑1 (MQ=255)                                                                                                        aaGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTGTTTTGGAATTATTTGATTTAGTATTCAATTTTCTCCGGGCGGTTTCCCGCCAACATTTTTTAATTTGTGGGCGT‑TGAAAAAAAGCCTTGATAGAAGACcactg                                   >  2:321184/1‑138 (MQ=255)                                                                                                        aaGTCTTGCAGCGCCGTGAGGGGATTCATCGGCGTGTTATCGAATTCTTTCAATTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCtt                                   <  4:4827/139‑1 (MQ=255)                                                                                                         aGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCtt                                   <  7:152913/138‑1 (MQ=255)                                                                                                         aGTCTAGCAGCGCCTTGAGGGTAATCATCGGCGTTTTTTCGAATTCTTTCCTTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCATACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCtt                                   <  3:35978/138‑1 (MQ=255)                                                                                                          gTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCa                                 <  4:65735/139‑1 (MQ=255)                                                                                                           tCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCAc                                <  7:255924/139‑1 (MQ=255)                                                                                                           tCTTGCAGCGCATTGAGGGTATTCATCGGCGTTTTTTAGAATTCTTTCATTTCCTCTTCAATTTGATCCGTGCTTTTAACCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGAACTCTTCAc                                <  7:344286/139‑1 (MQ=255)                                                                                                                cAGCGCCTTGAGGGTATTCATCGGCGTTTTTTGGAATTATTTCATTTACTCTTCAATTTGCTCCGGGGTTTTCCCCGCAAACATTTTTTTA‑TTGCGAGCGTATGAAACAACGCCGTGATACATGACCTCTTCACTCtt                            >  7:104346/1‑138 (MQ=255)                                                                                                                cAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTa                           <  4:164524/139‑1 (MQ=255)                                                                                                                     ccTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATg                      <  2:132372/139‑1 (MQ=255)                                                                                                                     ccTTGAGGGTATTCATCGGCGTTTTTTCGAATTATTTAATTTACTCTTCAATTTTCTACGGGCTTTTACCCTCAAACTTTTTTTAA‑TTTCGAGCGTATGAAAGAACGCCTAGATACATGACCTCTTCAATCTTAatat                       >  2:234040/1‑138 (MQ=255)                                                                                                                     acTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATg                      <  8:299933/138‑1 (MQ=255)                                                                                                                         gAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCAt                   <  6:314479/138‑1 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTgcg            <  1:240473/138‑1 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTgcg            <  1:256724/138‑1 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTgcg            <  2:367299/138‑1 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTgcg            <  1:234040/138‑1 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTACCCTCAAACTTTTTGTAA‑TTTCGAGCGTATGAAAGAAAGCCATGATACATGACCGCGTCACTCTTAAAATTCATAGGTgcgc           >  7:107169/1‑139 (MQ=255)                                                                                                                                ttCATCGGCGTTTTTTCGAATTATTTGATTTCCACTTGAATTTGCTCCGGGCTTTTCCCCGCCAACTTTTTTTAA‑TTGCGAGCGTATGAAACTACGCCATGATACATGACCTCTTCAATCTAAATATGCAACGGTgcg            >  3:186456/1‑138 (MQ=255)                                                                                                                                ttCATCGGCGTGTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTgcg            <  6:311040/138‑1 (MQ=255)                                                                                                                                 tCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTACACTCTTAATATGCATCGGTgcgc           <  7:93749/138‑1 (MQ=255)                                                                                                                                         gTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTGCGCACATTGGGa  <  5:155996/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          AATAATATTGGCGCTTTCGGTATTAATCATGTGGTCATGACGGGCCTGCACCACAAATGTCGGTGAATAAATCATATCGACATTATTCCGCACATCAGCAATTAAGTCTTGCAGCGCCTTGAGGGTATTCATCGGCGTTTTTTCGAATTCTTTCATTTCCTCTTCAATTTGCTCCGGGCTTTTCCCCTCAAACTTTTTGTAA‑TTGCGAGCGTATGAAAGAACGCCTTGATACATGACCTCTTCACTCTTAATATGCATCGGTGCGCACATTGGGA  >  CP006881/3272207‑3272481

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