Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	499,459	A→C	6.5%	L172R (CTC→CGC) ‡	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,459	0	A	C	6.5%	210.9 / 7.0	62	L172R (CTC→CGC) ‡	aes	acetyl esterase
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base C (2/2);  total (28/34)
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.26e-01

CTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGT  >  NC_000913/499317‑499604                                                                                                                                               |                                                                                                                                                  cTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCa                                                                                                                                             <  4:211431/149‑1 (MQ=255)  tGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAg                                                                                                                                            >  4:578019/1‑149 (MQ=255)               cATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCg                                                                                                                               >  4:481036/1‑149 (MQ=255)                aTCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCgg                                                                                                                              >  1:158780/1‑149 (MQ=255)                        cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCg                                                                                                                      <  4:72864/149‑1 (MQ=255)                               cAACAGACGACGAGTCACGGAATCCCGTAATCCGTATAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGc                                                                                                               <  3:390494/149‑1 (MQ=255)                               cAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGc                                                                                                               <  3:220494/149‑1 (MQ=255)                                aaCAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAgc                                                                                                                                                 <  3:164498/114‑1 (MQ=255)                                aaCAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAgc                                                                                                                                                 >  4:164498/1‑114 (MQ=255)                                aaCAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGcc                                                                                                              >  3:415485/1‑149 (MQ=255)                                   aGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAAt                                                                                                           >  4:287269/1‑149 (MQ=255)                                    gacgacGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATg                                                                                                          >  2:265504/1‑149 (MQ=255)                                       gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCgg                                                                                                       <  4:271304/149‑1 (MQ=255)                                             tCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACata                                                                                                 >  4:205462/1‑149 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGCTCGCCAgt                                                                                                                                           >  1:28941/1‑104 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGCTCGCCAgt                                                                                                                                           <  2:28941/105‑2 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGcc                                                                                                              >  1:451463/1‑134 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGcc                                                                                                              <  2:451463/134‑1 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACat                                                                                                  >  2:144615/1‑146 (MQ=255)                                               aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACat                                                                                                  <  1:144615/146‑1 (MQ=255)                                                  gAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgat                                                                                            <  4:557036/149‑1 (MQ=255)                                                  gAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgat                                                                                            <  3:481036/149‑1 (MQ=255)                                                   aaTCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatt                                                                                           <  3:178771/149‑1 (MQ=255)                                                      cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttga                                                                                        <  1:13280/149‑1 (MQ=255)                                                       ccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttgat                                                                                       >  4:175203/1‑149 (MQ=255)                                                           aaTCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtc                                                                                   <  1:265504/149‑1 (MQ=255)                                                            aTCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttga                                                                                        <  4:544914/143‑1 (MQ=255)                                                            aTCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttga                                                                                        >  3:544914/1‑143 (MQ=255)                                                              ccGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcctc                                                                                <  4:367525/149‑1 (MQ=255)                                                                gTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCaaa                                                                                                                 >  3:132138/1‑114 (MQ=255)                                                                gTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCaaa                                                                                                                 <  4:132138/114‑1 (MQ=255)                                                                gTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCg                                                                              >  3:529183/1‑149 (MQ=255)                                                                   aaGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCTCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCct                                                                           <  3:505724/149‑1 (MQ=255)                                                                      cccATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCctgct                                                                        >  4:131124/1‑149 (MQ=255)                                                                            ccACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGa                                                                  >  1:355630/1‑149 (MQ=255)                                                                               caGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAAt                                                               <  1:270537/149‑1 (MQ=255)                                                                                aGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATa                                                              >  3:406806/1‑149 (MQ=255)                                                                                      aCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAacaaca                                                        <  1:365222/149‑1 (MQ=255)                                                                                        gCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAag                                                      >  4:489724/1‑149 (MQ=255)                                                                                         cccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagc                                                     <  2:28533/149‑1 (MQ=255)                                                                                              aaCTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGTAGCca                                                >  4:377303/1‑149 (MQ=255)                                                                                                  ttACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAAt                                            <  4:415485/149‑1 (MQ=255)                                                                                                               tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCt                               <  3:284038/149‑1 (MQ=255)                                                                                                               tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCt                               >  1:243588/1‑149 (MQ=255)                                                                                                                 tGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTg                             >  3:376432/1‑149 (MQ=255)                                                                                                                  gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                            <  3:157156/149‑1 (MQ=255)                                                                                                                  gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                            >  2:407291/1‑149 (MQ=255)                                                                                                                   tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCg                           >  2:570877/1‑149 (MQ=255)                                                                                                                      aTCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                        >  3:186781/1‑149 (MQ=255)                                                                                                                              aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCt                <  3:511756/149‑1 (MQ=255)                                                                                                                                    aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGt          <  2:158780/149‑1 (MQ=255)                                                                                                                                     catgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                        <  2:476474/114‑1 (MQ=255)                                                                                                                                     catgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                        >  1:476474/5‑118 (MQ=255)                                                                                                                                      gCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaaa                                                                >  3:164491/1‑93 (MQ=255)                                                                                                                                      gCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaaa                                                                <  4:164491/93‑1 (MQ=255)                                                                                                                                            gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgt  <  1:396370/149‑1 (MQ=255)                                                                                                                                            gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgt  <  3:3452/149‑1 (MQ=255)                                                                                                                                            gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgt  <  1:564435/149‑1 (MQ=255)                                                                                                                                            gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgt  <  3:205462/149‑1 (MQ=255)                                                                                                                                            gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgt  <  3:287269/149‑1 (MQ=255)                                                                                                                                              gAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                         >  2:45283/1‑108 (MQ=255)                                                                                                                                              gAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                         <  1:45283/108‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  CTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGT  >  NC_000913/499317‑499604

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