Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,397,366 A→T 18.8% S49S (TCT→TCA abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,397,3660AT18.8% 480.5 / 104.3 197S49S (TCT→TCAabgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (77/83);  new base T (28/9);  total (105/92)
Fisher's exact test for biased strand distribution p-value = 3.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 197 total aligned reads displayed.

GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCT  >  CP009273/1397220‑1397508
                                                                                                                                                  |                                                                                                                                              
gTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAg                                                                                                                                               <  1:342578/148‑1 (MQ=255)
    gTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTc                                                                                                                                           <  1:3238137/148‑1 (MQ=255)
    gTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCAGTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTc                                                                                                                                           >  2:929944/1‑148 (MQ=255)
    gTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTGGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTc                                                                                                                                           <  1:3428376/148‑1 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAg                                                                                                                                >  2:664317/1‑148 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATACCAGCGCCGAAg                                                                                                                                >  2:2072266/1‑148 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCTATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  1:2160091/148‑1 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  2:2526414/148‑1 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  1:557970/148‑1 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  1:405146/148‑1 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  2:606780/148‑1 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                            <  1:231290/148‑1 (MQ=255)
                      gctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATg                                                                                                                         <  1:792248/148‑1 (MQ=255)
                        tGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCt                                                                                                                       <  2:308857/148‑1 (MQ=255)
                        tGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCt                                                                                                                       <  2:380939/148‑1 (MQ=255)
                         gACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTc                                                                                                                      >  2:1634720/1‑148 (MQ=255)
                          aCTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTcc                                                                                                                     >  1:714008/1‑148 (MQ=255)
                                aCCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGAc                                                                                                               >  2:2344000/1‑148 (MQ=255)
                                    gCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGa                                                                                                           <  2:3662222/148‑1 (MQ=255)
                                       aGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAct                                                                                                        >  2:373536/1‑148 (MQ=255)
                                       aGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAct                                                                                                        <  1:2183453/148‑1 (MQ=255)
                                       aGGGCGTCATATTCGCCCAGCAGGGCGATAACCGGTTGGCCTTGACCACACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAct                                                                                                        <  2:2224336/148‑1 (MQ=255)
                                        gggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActc                                                                                                       <  2:1792953/148‑1 (MQ=255)
                                         ggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActct                                                                                                      >  1:2911413/1‑148 (MQ=255)
                                         ggCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAActct                                                                                                      >  2:3649751/1‑148 (MQ=255)
                                          gCGTCATAGTCTCCCAGCAGGGCGATAACCGGTGTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCtt                                                                                                     <  2:2610928/148‑1 (MQ=255)
                                           cGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTc                                                                                                    <  1:3032529/148‑1 (MQ=255)
                                              cATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCaaa                                                                                                 <  1:3813544/148‑1 (MQ=255)
                                                tatTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAAc                                                                                               >  2:3800847/1‑148 (MQ=255)
                                                 atTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACg                                                                                              <  1:2285073/148‑1 (MQ=255)
                                                  ttCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgt                                                                                             >  2:869212/1‑148 (MQ=255)
                                                   tctcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgtg                                                                                            <  2:2832785/148‑1 (MQ=255)
                                                    ctcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgtgt                                                                                           <  2:1165547/148‑1 (MQ=255)
                                                    ctcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgtgt                                                                                           <  2:1208748/148‑1 (MQ=255)
                                                     tcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGtt                                                                                          <  2:2520338/148‑1 (MQ=255)
                                                     tcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGtt                                                                                          <  2:3656552/148‑1 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                         >  2:3523378/1‑148 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                         <  2:1308178/148‑1 (MQ=255)
                                                      cccAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGttt                                                                                         >  1:235681/1‑148 (MQ=255)
                                                         agcagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTg                                                                                      <  2:3229808/148‑1 (MQ=255)
                                                         agcagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTg                                                                                      <  2:1324963/148‑1 (MQ=255)
                                                          gcagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTgg                                                                                     >  1:256242/1‑148 (MQ=255)
                                                          gcagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTgg                                                                                     >  2:3045234/1‑148 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                    >  2:1446408/1‑148 (MQ=255)
                                                            agGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgat                                                                                   <  1:2727089/148‑1 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                 >  1:129977/1‑148 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                 >  1:1888621/1‑148 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                 >  1:2197136/1‑148 (MQ=255)
                                                               gCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGACAGATGCTCCGCTGACCAGAAATCTTCAAAACGTGTTTCTGgatgat                                                                                <  1:896824/148‑1 (MQ=255)
                                                                  aTAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATccc                                                                             >  2:3626982/1‑148 (MQ=255)
                                                                    aaCCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAg                                                                           >  2:1087425/1‑148 (MQ=255)
                                                                     aCCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGa                                                                          <  1:3765812/148‑1 (MQ=255)
                                                                      ccGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAt                                                                         <  1:280604/148‑1 (MQ=255)
                                                                       cGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAtt                                                                        <  1:1303227/148‑1 (MQ=255)
                                                                           tttccTTCACAAAACGAAGCAATAAATGAATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGAGACGAAGCCAGATGCTCAGCTGACCAGAACTCTTCAAAACGAGTTACTGGATGATCCCAGATTTGAt                                                                    >  2:2018240/4‑148 (MQ=255)
                                                                           ttGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGAt                                                                    >  1:2957599/1‑148 (MQ=255)
                                                                                  gACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAat                                                             <  1:1791448/148‑1 (MQ=255)
                                                                                  gACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAat                                                             >  2:3292946/1‑148 (MQ=255)
                                                                                    ccAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAatat                                                           >  2:3432388/1‑148 (MQ=255)
                                                                                      aaaCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCg                                                         >  2:3266271/1‑148 (MQ=255)
                                                                                      aaaCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCg                                                         <  2:2794223/148‑1 (MQ=255)
                                                                                        aCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGt                                                       >  2:3606565/1‑148 (MQ=255)
                                                                                           aaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGtata                                                    >  1:2540221/1‑148 (MQ=255)
                                                                                           aaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGtata                                                    >  1:1592853/1‑148 (MQ=255)
                                                                                           aaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGtata                                                    <  2:71845/148‑1 (MQ=255)
                                                                                           aaGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGtata                                                    >  1:104734/1‑148 (MQ=255)
                                                                                              cAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACg                                                 <  1:1883106/148‑1 (MQ=255)
                                                                                                aTAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGtt                                               >  1:2819986/1‑148 (MQ=255)
                                                                                                 tAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTg                                              >  1:2266989/1‑148 (MQ=255)
                                                                                                 tAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTg                                              >  1:2266084/1‑148 (MQ=255)
                                                                                                 tAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTg                                              >  1:2141090/1‑148 (MQ=255)
                                                                                                 tAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTg                                              >  1:600627/1‑148 (MQ=255)
                                                                                                  aaaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTgg                                             >  1:1327035/1‑148 (MQ=255)
                                                                                                   aaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGc                                            >  1:1369987/1‑148 (MQ=255)
                                                                                                   aaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGc                                            >  2:2019508/1‑148 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                           >  1:3600264/1‑148 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                           >  1:1703473/1‑148 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                           >  2:478106/1‑148 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                           >  2:98787/1‑148 (MQ=255)
                                                                                                      gCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGCGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGAt                                         <  1:1969055/148‑1 (MQ=255)
                                                                                                         tttGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCgg                                      >  1:2841163/1‑148 (MQ=255)
                                                                                                            gggATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCtt                                   >  2:487647/1‑148 (MQ=255)
                                                                                                               atatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCaa                                >  1:2421465/1‑148 (MQ=255)
                                                                                                                tatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAAt                               <  1:1394109/148‑1 (MQ=255)
                                                                                                                tatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAAt                               <  1:694765/148‑1 (MQ=255)
                                                                                                                 atTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATc                              >  2:2407894/1‑148 (MQ=255)
                                                                                                                  ttGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCg                             >  1:1587717/1‑148 (MQ=255)
                                                                                                                  ttGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCg                             >  1:301414/1‑148 (MQ=255)
                                                                                                                       tACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcg                        <  2:2933727/148‑1 (MQ=255)
                                                                                                                       tACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcg                        >  1:3010871/1‑148 (MQ=255)
                                                                                                                        aCGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                       <  2:2632716/148‑1 (MQ=255)
                                                                                                                           ttGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGa                    >  1:1260581/1‑148 (MQ=255)
                                                                                                                            tGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGAt                   <  1:3571310/148‑1 (MQ=255)
                                                                                                                               gggTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaa                <  2:2148249/148‑1 (MQ=255)
                                                                                                                                ggTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaaa               <  1:3477049/148‑1 (MQ=255)
                                                                                                                                    aCGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGAt           >  1:3086452/1‑148 (MQ=255)
                                                                                                                                    aCGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTAGATAAAACGAt           <  1:1249974/148‑1 (MQ=255)
                                                                                                                                         gAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGAtt      <  2:2667662/148‑1 (MQ=255)
                                                                                                                                           aGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTc    <  2:1922981/148‑1 (MQ=255)
                                                                                                                                             ccTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCt  >  2:1606782/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCT  >  CP009273/1397220‑1397508

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: