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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273495,0200TA99.1% 702.7 / ‑2.7 222E137V (GAA→GTA) aesacetyl esterase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (99/120);  total (99/123)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 222 total aligned reads displayed.

CAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494876‑495164
                                                                                                                                                 |                                                                                                                                                
cAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACct                                                                                                                                               >  2:1725592/1‑148 (MQ=255)
cAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACct                                                                                                                                               <  2:2066176/148‑1 (MQ=255)
  aTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACctct                                                                                                                                             >  1:429523/1‑148 (MQ=255)
  aTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACctct                                                                                                                                             <  2:2373490/148‑1 (MQ=255)
    cGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTAt                                                                                                                                           <  1:2674474/148‑1 (MQ=255)
      aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCg                                                                                                                                         >  2:3046018/1‑148 (MQ=255)
      aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCg                                                                                                                                         >  2:67985/1‑148 (MQ=255)
       tCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGc                                                                                                                                        <  2:2256447/148‑1 (MQ=255)
        cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  1:1389417/1‑148 (MQ=255)
        cTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCt                                                                                                                                       >  1:890956/1‑148 (MQ=255)
           tttATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGc                                                                                                                                    >  2:1647715/1‑148 (MQ=255)
                cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaa                                                                                                                               <  2:1119191/148‑1 (MQ=255)
                 aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGaaaa                                                                                                                              >  1:2198285/1‑148 (MQ=255)
                    cAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAAcgc                                                                                                                           >  2:3316424/1‑148 (MQ=255)
                      aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACgcgc                                                                                                                         <  2:2744966/148‑1 (MQ=255)
                      aCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACgcgc                                                                                                                         <  2:2018608/148‑1 (MQ=255)
                       ccACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCt                                                                                                                        <  1:2815075/148‑1 (MQ=255)
                         acaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTc                                                                                                                      <  2:2940436/148‑1 (MQ=255)
                          caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     <  1:257856/148‑1 (MQ=255)
                          caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     >  2:893487/1‑148 (MQ=255)
                          caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCa                                                                                                                     >  1:1787274/1‑148 (MQ=255)
                           aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAg                                                                                                                    >  2:2913921/1‑148 (MQ=255)
                            aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAgg                                                                                                                   <  2:2619384/148‑1 (MQ=255)
                            aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAgg                                                                                                                   <  2:976812/148‑1 (MQ=255)
                             cGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGt                                                                                                                  <  2:1000505/148‑1 (MQ=255)
                                    gCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAggg                                                                                                           >  1:2790998/1‑148 (MQ=255)
                                      gAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGgtg                                                                                                         <  1:3467580/148‑1 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      <  2:2444927/148‑1 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      >  2:1970785/1‑148 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      <  2:1234437/148‑1 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      <  2:3212117/148‑1 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      <  2:3362417/148‑1 (MQ=255)
                                         cgcCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTaa                                                                                                      <  2:3508882/148‑1 (MQ=255)
                                            cAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCa                                                                                                   >  1:2650246/1‑148 (MQ=255)
                                              gCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAAt                                                                                                 <  2:2828414/148‑1 (MQ=255)
                                                    cACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAAt                                                                                           <  1:2429695/148‑1 (MQ=255)
                                                     aCCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCACACCGTACCTCTATCTATTTCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATc                                                                                          <  1:3447689/148‑1 (MQ=255)
                                                           ggAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTa                                                                                    <  2:455672/148‑1 (MQ=255)
                                                             aaTCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACa                                                                                  <  1:2968357/148‑1 (MQ=255)
                                                              aTCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAt                                                                                 >  1:655507/1‑148 (MQ=255)
                                                               tCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAtt                                                                                >  2:406780/1‑148 (MQ=255)
                                                               tCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAtt                                                                                <  1:448396/148‑1 (MQ=255)
                                                               tCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCAACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACAtt                                                                                >  2:2600377/1‑148 (MQ=255)
                                                                 aCCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTgg                                                                              >  1:1069920/1‑148 (MQ=255)
                                                                     gCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGt                                                                          <  2:710190/148‑1 (MQ=255)
                                                                      cAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTa                                                                         >  2:1780629/1‑148 (MQ=255)
                                                                        aaGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGc                                                                       <  2:2786897/148‑1 (MQ=255)
                                                                        aaGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGc                                                                       <  1:1006089/148‑1 (MQ=255)
                                                                            cAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGc                                                                   <  2:2068145/148‑1 (MQ=255)
                                                                            cAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGc                                                                   <  2:1230419/148‑1 (MQ=255)
                                                                               tGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcag                                                                >  1:1189999/1‑148 (MQ=255)
                                                                                gCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcagc                                                               >  1:1538896/1‑148 (MQ=255)
                                                                                 cGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcagca                                                              <  1:3218212/148‑1 (MQ=255)
                                                                                 cGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCcagca                                                              <  1:3221427/148‑1 (MQ=255)
                                                                                   ggACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAgg                                                            >  1:1233380/1‑148 (MQ=255)
                                                                                      cATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGgcgc                                                         <  1:339330/148‑1 (MQ=255)
                                                                                       atatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACCAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCa                                                        <  1:1588537/148‑1 (MQ=255)
                                                                                       atatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCa                                                        >  2:2760244/1‑148 (MQ=255)
                                                                                        tatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                       >  1:1025546/1‑148 (MQ=255)
                                                                                        tatTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCat                                                       >  1:1461802/1‑148 (MQ=255)
                                                                                         atggatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatg                                                      <  2:1952080/145‑1 (MQ=255)
                                                                                            gatttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCatgatg                                                   >  2:2855412/1‑148 (MQ=255)
                                                                                               ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:1743763/1‑148 (MQ=255)
                                                                                               ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:485284/1‑148 (MQ=255)
                                                                                               ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                <  1:2610747/148‑1 (MQ=255)
                                                                                               ttgatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGa                                                >  1:933554/1‑148 (MQ=255)
                                                                                                 gatAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATc                                              <  2:1389417/148‑1 (MQ=255)
                                                                                                    aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTg                                           >  2:3337543/1‑148 (MQ=255)
                                                                                                    aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTg                                           <  1:2465046/148‑1 (MQ=255)
                                                                                                    aaTCCTCCGCCTGCTGGTGGAAATAACAACAAGAAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTg                                           >  2:2551631/1‑148 (MQ=255)
                                                                                                      tcctccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTggg                                         <  2:2666733/148‑1 (MQ=255)
                                                                                                          ccGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCGTGGGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGtctc                                     >  2:1035580/1‑148 (MQ=255)
                                                                                                           cGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCg                                    <  1:1311636/148‑1 (MQ=255)
                                                                                                              ctgctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgaga                                 <  1:1345977/148‑1 (MQ=255)
                                                                                                              ctgctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCgaga                                 <  2:887833/148‑1 (MQ=255)
                                                                                                               tgctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGAt                                >  1:2726042/1‑148 (MQ=255)
                                                                                                                 ctgGTGGAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTg                              >  1:1997364/1‑148 (MQ=255)
                                                                                                                   ggtggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGcc                            >  1:3064485/1‑148 (MQ=255)
                                                                                                                    gtggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCg                           <  1:69135/148‑1 (MQ=255)
                                                                                                                     tggAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCga                          <  1:1286479/148‑1 (MQ=255)
                                                                                                                     tggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCga                          >  2:925070/1‑148 (MQ=255)
                                                                                                                      ggAAATAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCgag                         <  2:2777317/148‑1 (MQ=255)
                                                                                                                          aTAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATa                     <  1:1970785/148‑1 (MQ=255)
                                                                                                                           tAACAACAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaa                    >  1:608278/1‑148 (MQ=255)
                                                                                                                            aacaacAAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaa                   <  2:2160534/148‑1 (MQ=255)
                                                                                                                             acaacaAGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATaaaa                  <  2:2341056/148‑1 (MQ=255)
                                                                                                                              caacaaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAc                 <  1:392787/148‑1 (MQ=255)
                                                                                                                              caacaaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAAc                 <  1:392951/148‑1 (MQ=255)
                                                                                                                                 caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              <  1:3316424/148‑1 (MQ=255)
                                                                                                                                 caaGCAGCCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGc              >  1:791086/1‑148 (MQ=255)
                                                                                                                                   agcagcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAATGGGTGTAATAAATACCAATCCCCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGTTATTGCCGAGAATAAAACCGCCt            <  2:1217867/148‑1 (MQ=255)
                                                                                                                                   agcagcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCt            >  1:779579/1‑148 (MQ=255)
                                                                                                                                     cagcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTcc          >  2:133480/1‑148 (MQ=255)
                                                                                                                                       gcCACAATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCAt        >  2:1578839/1‑148 (MQ=255)
                                                                                                                                         cacaATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCTGGCGCTTGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGc      <  1:441659/148‑1 (MQ=255)
                                                                                                                                         cacaATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGc      >  1:3146885/1‑148 (MQ=255)
                                                                                                                                           caATTCCCTCTCTGGCTTGCGGAAAACGCGCTTCTGTTGAAAGGTTGTAATCAATACCATTCACGGTACTTGGGCTGTAGGTTGCCAGCCGGCGCATTATGCGTTCGTGGGTATCTAGTTTGCCGAGAATAAAACCGCCTCCATGCaa    <  1:2138176/148‑1 (MQ=255)
                                                                                                                                           caATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaa    <  1:3046018/148‑1 (MQ=255)
                                                                                                                                           caATTACCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCaa    <  1:2905232/148‑1 (MQ=255)
                                                                                                                                              aTTACCTCTATCTGTTGCTGAAAACGCGCTTCTGGTGAATGGGTGTAATCAATACAAATCACCGTAAATTGGCTTTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGATAATAAAACCGCCTCCATGCaaat  <  1:1932238/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTGTAATCAATACCAATCACCGTACATTGGCTGTAGCTTGCCAGCAGGCGCATGATGCGATCGTGGGTATCGAGATTGCCGAGAATAAAACCGCCTCCATGCAAAT  >  CP009273/494876‑495164

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: