Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,805,420 C→A 12.7% A212E (GCG→GAG)  emrA → multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,805,4200CA12.7% 942.6 / 75.3 283A212E (GCG→GAG) emrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (102/145);  new base A (25/11);  total (127/156)
Fisher's exact test for biased strand distribution p-value = 2.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 308 total aligned reads displayed.

CCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGAC  >  CP009273/2805274‑2805567
                                                                                                                                                  |                                                                                                                                                   
  gcgACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTg                                                                                                                                                 >  2:3817754/1‑149 (MQ=255)
  gcgACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTg                                                                                                                                                 <  1:2408633/149‑1 (MQ=255)
      cGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGc                                                                                                                                             <  2:4186330/149‑1 (MQ=255)
        ccGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGt                                                                                                                                           <  2:2278495/149‑1 (MQ=255)
           tCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACACTGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACt                                                                                                                                        <  2:1218118/149‑1 (MQ=255)
           tCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACt                                                                                                                                        <  2:853932/149‑1 (MQ=255)
                   gCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGtattat                                                                                                                                >  1:4095917/1‑149 (MQ=255)
                        ggCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTc                                                                                                                           <  1:121340/149‑1 (MQ=255)
                        ggCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTc                                                                                                                           <  2:2526902/149‑1 (MQ=255)
                         gcgcAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTcc                                                                                                                          <  1:3033784/149‑1 (MQ=255)
                         gcgcAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTcc                                                                                                                          <  1:954658/149‑1 (MQ=255)
                            cAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGAt                                                                                                                       <  1:3579448/149‑1 (MQ=255)
                              aCTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGa                                                                                                                     <  2:4263264/149‑1 (MQ=255)
                                tGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGAcc                                                                                                                   <  2:3008392/149‑1 (MQ=255)
                                 ggACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCg                                                                                                                  >  1:1427978/1‑149 (MQ=255)
                                 ggACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCg                                                                                                                  <  2:1827344/149‑1 (MQ=255)
                                 ggACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCg                                                                                                                  >  1:2709160/1‑149 (MQ=255)
                                   aCGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGt                                                                                                                <  1:1935891/149‑1 (MQ=255)
                                     gTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTa                                                                                                              >  2:3887513/1‑149 (MQ=255)
                                     gTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTa                                                                                                              <  2:2058518/149‑1 (MQ=255)
                                      tCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTAt                                                                                                             >  2:1451561/1‑149 (MQ=255)
                                         cgATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGtc                                                                                                          <  2:2843911/149‑1 (MQ=255)
                                           aTTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGtctc                                                                                                        >  1:4741874/1‑149 (MQ=255)
                                            ttCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTcc                                                                                                       >  1:2527382/1‑149 (MQ=255)
                                            ttCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTcc                                                                                                       <  1:1423268/149‑1 (MQ=255)
                                            ttCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTcc                                                                                                       <  1:767355/149‑1 (MQ=255)
                                                acaatacaatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCccgcc                                                                                                   >  2:1535233/1‑149 (MQ=255)
                                                acaatacaatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCccgcc                                                                                                   >  2:4484431/1‑149 (MQ=255)
                                                  aatacaatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCcgccgc                                                                                                 <  1:3525208/149‑1 (MQ=255)
                                                    tacaatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCcgcgc                                                                                               >  2:695437/1‑149 (MQ=255)
                                                     acaatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCcgcgcg                                                                                              <  1:3161728/149‑1 (MQ=255)
                                                      caatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCgg                                                                                             <  2:4290817/149‑1 (MQ=255)
                                                       aatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGt                                                                                            <  1:1339970/149‑1 (MQ=255)
                                                       aatGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGt                                                                                            >  1:1824968/1‑149 (MQ=255)
                                                         tGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTAc                                                                                          <  1:3767159/149‑1 (MQ=255)
                                                         tGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTAc                                                                                          >  2:2742778/1‑149 (MQ=255)
                                                          gCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACa                                                                                         <  1:689319/149‑1 (MQ=255)
                                                          gCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACa                                                                                         >  1:3484247/1‑149 (MQ=255)
                                                          gCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCAGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACa                                                                                         >  1:851829/1‑149 (MQ=255)
                                                             aaTCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGcc                                                                                      <  2:317257/149‑1 (MQ=255)
                                                              aTCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCt                                                                                     >  1:1933465/1‑149 (MQ=255)
                                                               tCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTg                                                                                    <  2:1628581/149‑1 (MQ=255)
                                                                 aGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTggg                                                                                  >  2:3433877/1‑149 (MQ=255)
                                                                   gcgcTGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTTACGCCTGGCTGGCGCTGGTGCGTTATAGTATTATCCGTCCGATGACCGTTTATGTCTCCCGCCGCGCGGTACATCCGGGGgc                                                                                <  2:4684251/149‑1 (MQ=255)
                                                                     gatgatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCACGCCGCGCGGTACAGCCTGGGgcgc                                                                              >  1:4241939/1‑149 (MQ=255)
                                                                        gatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCaaa                                                                           <  1:2117337/149‑1 (MQ=255)
                                                                        gatTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCaaa                                                                           <  1:1966935/149‑1 (MQ=255)
                                                                         atTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAAt                                                                          >  1:1067312/1‑149 (MQ=255)
                                                                         atTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAAt                                                                          >  2:2704344/1‑149 (MQ=255)
                                                                              ggggACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                                     >  1:1030720/1‑149 (MQ=255)
                                                                              ggggACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                                     <  2:1735534/149‑1 (MQ=255)
                                                                              ggggACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGcc                                                                     >  2:2227819/1‑149 (MQ=255)
                                                                                 gACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCaa                                                                  <  1:2805360/149‑1 (MQ=255)
                                                                                 gACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCaa                                                                  >  1:148838/1‑149 (MQ=255)
                                                                                  aCTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAac                                                                 <  1:1662934/149‑1 (MQ=255)
                                                                                    tAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacga                                                               <  2:1380729/149‑1 (MQ=255)
                                                                                     aaaCTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAacgac                                                              <  1:3185137/149‑1 (MQ=255)
                                                                                        cTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGAcgcc                                                           >  1:580512/1‑149 (MQ=255)
                                                                                           gAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCt                                                        >  2:567766/1‑149 (MQ=255)
                                                                                              gATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGAt                                                     >  1:4349924/1‑149 (MQ=255)
                                                                                               aTCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                                    <  1:2368794/149‑1 (MQ=255)
                                                                                               aTCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGTCCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATg                                                    >  2:1744667/1‑149 (MQ=255)
                                                                                                 cagccagcCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                                  <  2:3340998/149‑1 (MQ=255)
                                                                                                 cagccagcCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCGGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGc                                                  <  2:3629285/149‑1 (MQ=255)
                                                                                                  agccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCCAATTAGCCCAACGACGCCGCTGATGGCg                                                 >  2:2441176/1‑149 (MQ=255)
                                                                                                    ccagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGt                                               <  2:2622818/149‑1 (MQ=255)
                                                                                                     cagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTc                                              >  2:3381103/1‑149 (MQ=255)
                                                                                                      agccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCg                                             >  2:2335953/1‑149 (MQ=255)
                                                                                                       gccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGt                                            >  2:2493238/1‑149 (MQ=255)
                                                                                                       gccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGt                                            >  2:438403/1‑149 (MQ=255)
                                                                                                       gccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGt                                            <  2:1771471/149‑1 (MQ=255)
                                                                                                        ccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGtt                                           <  2:613295/149‑1 (MQ=255)
                                                                                                         cgagCACCAGGCTGCCACCGACGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTCTCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTc                                          <  2:2322303/146‑1 (MQ=255)
                                                                                                            gCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAg                                       >  1:3408377/1‑149 (MQ=255)
                                                                                                                   gCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCaa                                <  1:219205/149‑1 (MQ=255)
                                                                                                                    cTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAat                               <  2:1332759/149‑1 (MQ=255)
                                                                                                                     tGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAata                              <  1:1868588/149‑1 (MQ=255)
                                                                                                                     tGCCACCGAAGTACGTAACGCCGGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAata                              <  1:4152375/149‑1 (MQ=255)
                                                                                                                     tGCCACCGAAGTACGTAACGCATGGCTGGCGCTGTAGCGTACTCTTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAata                              <  1:3900191/149‑1 (MQ=255)
                                                                                                                      gCCACCGAAGTACGTAACGCCTGGCTGGAGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAatat                             <  2:2148332/149‑1 (MQ=255)
                                                                                                                       ccACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtg                            <  1:2738548/149‑1 (MQ=255)
                                                                                                                       ccACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATAtg                            <  1:3943525/149‑1 (MQ=255)
                                                                                                                          ccGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgg                         <  1:329669/149‑1 (MQ=255)
                                                                                                                           cGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCACTATgtggg                        >  2:817165/1‑149 (MQ=255)
                                                                                                                           cGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtggg                        <  2:110900/149‑1 (MQ=255)
                                                                                                                             aaGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATgtgggtg                      <  1:3603261/149‑1 (MQ=255)
                                                                                                                                tACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGAt                   <  2:2946330/149‑1 (MQ=255)
                                                                                                                                  cGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGTGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGc                 <  2:1953893/149‑1 (MQ=255)
                                                                                                                                  cGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGc                 <  2:2656311/149‑1 (MQ=255)
                                                                                                                                  cGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGc                 <  1:4521470/149‑1 (MQ=255)
                                                                                                                                   gTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGcc                <  2:3919721/149‑1 (MQ=255)
                                                                                                                                     aaCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCaa              >  1:439716/1‑149 (MQ=255)
                                                                                                                                      aCGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAAc             <  1:4539275/149‑1 (MQ=255)
                                                                                                                                       cGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACt            <  2:1067312/149‑1 (MQ=255)
                                                                                                                                         ccTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACttt          <  2:4112633/149‑1 (MQ=255)
                                                                                                                                           tggctggaGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaa        >  2:4037255/1‑149 (MQ=255)
                                                                                                                                           gggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaa        <  2:40031/148‑1 (MQ=255)
                                                                                                                                            ggctggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTaaa       <  2:523844/149‑1 (MQ=255)
                                                                                                                                                ggcGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAgaga   >  1:1942935/1‑149 (MQ=255)
                                                                                                                                                 gcgcTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGAc  >  2:2025474/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCGCGACGCCGTCACCAGTGCCCAGGCGCAACTGGACGTCGCGATTCAACAATACAATGCCAATCAGGCGATGATTCTGGGGACTAAACTGGAAGATCAGCCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCGTATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGACGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGAC  >  CP009273/2805274‑2805567

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: