Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 495,020 T→A 6.7% E137V (GAA→GTA)  aes ← acetyl esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273495,0200TA 6.7% 776.1 / 32.9 235E137V (GAA→GTA) aesacetyl esterase
Reads supporting (aligned to +/- strand):  ref base T (105/114);  new base A (11/5);  total (116/119)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

TACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAA  >  CP009273/494871‑495167
                                                                                                                                                     |                                                                                                                                                   
tACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAAttt                                                                                                                                                     >  2:2707861/1‑150 (MQ=255)
tACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAAttt                                                                                                                                                     <  1:1729730/150‑1 (MQ=255)
tACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAAttt                                                                                                                                                     <  2:1650580/150‑1 (MQ=255)
tACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAAttt                                                                                                                                                     <  1:389634/150‑1 (MQ=255)
   cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAGTAACAACAAGCAGCCACAATTTCct                                                                                                                                                  <  1:4592522/150‑1 (MQ=255)
     cAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCctct                                                                                                                                                >  1:4514160/1‑150 (MQ=255)
        tCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATc                                                                                                                                             >  1:885410/1‑150 (MQ=255)
         cGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGTCACCAGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGAATTTATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATct                                                                                                                                            >  1:2143352/1‑149 (MQ=255)
           aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCt                                                                                                                                          >  2:3178593/1‑150 (MQ=255)
                tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCgg                                                                                                                                     >  1:3463594/1‑150 (MQ=255)
                tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCgg                                                                                                                                     >  1:4244363/1‑150 (MQ=255)
                 ttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGa                                                                                                                                    <  2:1344967/150‑1 (MQ=255)
                  tATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                                                                                                                                   <  2:3722709/150‑1 (MQ=255)
                     cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                                                                                                                                >  1:1767085/1‑150 (MQ=255)
                     cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                                                                                                                                <  1:2754973/150‑1 (MQ=255)
                     cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                                                                                                                                <  1:2859267/150‑1 (MQ=255)
                      aCGCAACCACACCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcg                                                                                                                               <  2:2131749/150‑1 (MQ=255)
                      aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcg                                                                                                                               <  1:715274/150‑1 (MQ=255)
                              acaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAg                                                                                                                       <  2:2114018/150‑1 (MQ=255)
                                aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATACCAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTGCAGGt                                                                                                                     >  2:3990068/1‑150 (MQ=255)
                                aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGt                                                                                                                     >  2:117516/1‑150 (MQ=255)
                                   gCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaa                                                                                                                  >  2:3377090/1‑150 (MQ=255)
                                   gCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaa                                                                                                                  >  2:3517445/1‑150 (MQ=255)
                                     aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg                                                                                                                <  1:46855/150‑1 (MQ=255)
                                     aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg                                                                                                                <  1:2356134/150‑1 (MQ=255)
                                     aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg                                                                                                                >  2:1000305/1‑150 (MQ=255)
                                     aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGaacg                                                                                                                >  2:2276840/1‑148 (MQ=255)
                                     aCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg                                                                                                                >  2:596569/1‑150 (MQ=255)
                                       tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAggg                                                                                                              >  2:82806/1‑150 (MQ=255)
                                            aGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATCTTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                         <  1:4381492/150‑1 (MQ=255)
                                            aGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                         >  2:1418963/1‑150 (MQ=255)
                                                 cAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                    >  2:2413250/1‑150 (MQ=255)
                                                 cAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                    <  1:698100/150‑1 (MQ=255)
                                                    cATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATAcc                                                                                                 <  2:353363/150‑1 (MQ=255)
                                                     aTGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCa                                                                                                >  1:3634855/1‑150 (MQ=255)
                                                     aTGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCa                                                                                                >  1:4541366/1‑150 (MQ=255)
                                                       ggCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAAt                                                                                              >  1:3014505/1‑150 (MQ=255)
                                                       ggCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAAt                                                                                              <  2:3934027/150‑1 (MQ=255)
                                                         cACCTGCGGAATCACCTGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCATGTGAAAGGGTGTAATCAATACCAATCa                                                                                            <  2:1262061/150‑1 (MQ=255)
                                                          aCCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCAc                                                                                           >  2:721663/1‑150 (MQ=255)
                                                            cTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCg                                                                                         >  2:3601534/1‑150 (MQ=255)
                                                               cGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTAc                                                                                      <  2:2969105/150‑1 (MQ=255)
                                                               cGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTAc                                                                                      <  1:2358414/150‑1 (MQ=255)
                                                                ggAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACa                                                                                     >  2:4213707/1‑150 (MQ=255)
                                                                ggAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACa                                                                                     >  1:1555840/1‑150 (MQ=255)
                                                                 gAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAt                                                                                    >  2:4300048/1‑150 (MQ=255)
                                                                 gAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAt                                                                                    >  2:4517932/1‑150 (MQ=255)
                                                                         ggCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTa                                                                            <  2:2578749/150‑1 (MQ=255)
                                                                         ggCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTa                                                                            <  2:2978958/150‑1 (MQ=255)
                                                                            aaaGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                         <  1:1189534/150‑1 (MQ=255)
                                                                            aaaGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCt                                                                         <  1:411074/150‑1 (MQ=255)
                                                                               gCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGc                                                                      >  1:3769470/1‑150 (MQ=255)
                                                                                  aaTGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcag                                                                   >  1:131191/1‑150 (MQ=255)
                                                                                  aaTGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATAACCGTACATTGGCTTTAGCTTGCcag                                                                   >  2:1916019/1‑150 (MQ=255)
                                                                                    tGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcagca                                                                 <  2:1297310/150‑1 (MQ=255)
                                                                                      cGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                               >  2:3066714/1‑150 (MQ=255)
                                                                                      cGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                               <  1:1974175/150‑1 (MQ=255)
                                                                                      cGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                               >  1:2565604/1‑150 (MQ=255)
                                                                                        ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcg                                                             >  1:789676/1‑150 (MQ=255)
                                                                                        ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcg                                                             <  2:714864/150‑1 (MQ=255)
                                                                                        ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcg                                                             <  1:3937795/150‑1 (MQ=255)
                                                                                        ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcg                                                             <  1:579045/150‑1 (MQ=255)
                                                                                         gACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcgc                                                            >  2:4396306/1‑150 (MQ=255)
                                                                                         gACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcgc                                                            >  1:1048235/1‑150 (MQ=255)
                                                                                           cATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                          >  1:2724892/1‑150 (MQ=255)
                                                                                           cATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                          >  2:2565753/1‑150 (MQ=255)
                                                                                            atatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatg                                                         >  1:4590140/1‑150 (MQ=255)
                                                                                              atTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatgat                                                       >  1:1409697/1‑150 (MQ=255)
                                                                                                  atttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                   <  2:3652582/150‑1 (MQ=255)
                                                                                                        tAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTgg                                             >  1:3573909/1‑150 (MQ=255)
                                                                                                         aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTggg                                            >  1:198991/1‑150 (MQ=255)
                                                                                                            cctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTAt                                         <  1:895161/150‑1 (MQ=255)
                                                                                                             ctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATc                                        <  1:8523/150‑1 (MQ=255)
                                                                                                              tacGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCg                                       <  1:3975545/148‑1 (MQ=255)
                                                                                                               ccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCga                                      >  2:982887/1‑150 (MQ=255)
                                                                                                                 gCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgaga                                    <  1:1892695/150‑1 (MQ=255)
                                                                                                                  ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGAt                                   >  2:2639506/1‑150 (MQ=255)
                                                                                                                  ccTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGAt                                   >  1:3387614/1‑150 (MQ=255)
                                                                                                                   ctgctgGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGAtt                                  >  2:3401612/1‑150 (MQ=255)
                                                                                                                   ctgctgGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGAtt                                  <  2:3476859/150‑1 (MQ=255)
                                                                                                                    tgctgGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                                 >  1:1864931/1‑150 (MQ=255)
                                                                                                                    tgctgGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                                 >  2:1732310/1‑150 (MQ=255)
                                                                                                                      ctgGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                               >  1:3094093/1‑150 (MQ=255)
                                                                                                                         gtggAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTCTCGAGATTGCCgag                            >  2:4056671/1‑150 (MQ=255)
                                                                                                                          tggAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCgaga                           >  1:1172758/1‑150 (MQ=255)
                                                                                                                          tggAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCgaga                           <  2:4243739/150‑1 (MQ=255)
                                                                                                                            gAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                         >  2:1921880/1‑150 (MQ=255)
                                                                                                                            gAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAAt                         >  1:2573346/1‑150 (MQ=255)
                                                                                                                              aaTAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTAGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaa                       >  1:4456987/1‑150 (MQ=255)
                                                                                                                                  acaacaAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAcc                   <  2:1055876/150‑1 (MQ=255)
                                                                                                                                    aacaaGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc                 >  1:542538/1‑150 (MQ=255)
                                                                                                                                     acaaGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGcc                >  2:2903493/1‑150 (MQ=255)
                                                                                                                                       aaGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTc              <  1:916859/150‑1 (MQ=255)
                                                                                                                                         gcagcCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCa            <  2:3471956/150‑1 (MQ=255)
                                                                                                                                           agcCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATg          <  1:3206327/150‑1 (MQ=255)
                                                                                                                                               acaATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaa      >  2:3166828/1‑150 (MQ=255)
                                                                                                                                                caATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaat     >  1:1074585/1‑150 (MQ=255)
                                                                                                                                                 aaTTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaata    <  1:4145252/150‑1 (MQ=255)
                                                                                                                                                   tttCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataaa  >  1:2491422/1‑150 (MQ=255)
                                                                                                                                                   tttCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaaataaa  >  1:2828514/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAATAAA  >  CP009273/494871‑495167

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: