Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,149,603 A→T 91.0% E287V (GAA→GTA)  mdtB → multidrug efflux system, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,149,6030AT91.0% 1084.2 / 74.7 321E287V (GAA→GTA) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base A (10/19);  new base T (152/140);  total (162/159)
Fisher's exact test for biased strand distribution p-value = 8.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

TTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGCACCACCAATA  >  CP009273/2149457‑2149749
                                                                                                                                                  |                                                                                                                                                  
ttCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACa                                                                                                                                                 <  2:1723873/150‑1 (MQ=255)
  ccGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAgg                                                                                                                                               >  2:2465866/1‑150 (MQ=255)
      gAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCAt                                                                                                                                           <  1:875829/150‑1 (MQ=255)
      gAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCAt                                                                                                                                           >  1:1719774/1‑150 (MQ=255)
           aCCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGtga                                                                                                                                      <  2:1359377/150‑1 (MQ=255)
            ccAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGtgat                                                                                                                                     >  2:3288351/1‑150 (MQ=255)
            ccAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGtgat                                                                                                                                     <  2:740215/150‑1 (MQ=255)
              aGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGtgatga                                                                                                                                   >  2:1309332/1‑150 (MQ=255)
              aGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGtgatga                                                                                                                                   >  2:3334838/1‑150 (MQ=255)
                    aaTCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTc                                                                                                                             >  1:3380561/1‑150 (MQ=255)
                      tCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAg                                                                                                                           <  1:4102429/150‑1 (MQ=255)
                        cgccgAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAgcg                                                                                                                         <  2:1472001/150‑1 (MQ=255)
                           cgAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCa                                                                                                                      >  1:4304824/1‑150 (MQ=255)
                            gAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAg                                                                                                                     >  2:3551553/1‑150 (MQ=255)
                            gAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAg                                                                                                                     >  2:3593390/1‑150 (MQ=255)
                             aaGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGc                                                                                                                    >  1:2102134/1‑150 (MQ=255)
                              agagTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGcc                                                                                                                   >  2:4034687/1‑150 (MQ=255)
                               gagTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGccc                                                                                                                  <  2:1977714/150‑1 (MQ=255)
                               gagTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGccc                                                                                                                  <  2:3425100/150‑1 (MQ=255)
                                agTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACCGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCg                                                                                                                 <  1:740459/150‑1 (MQ=255)
                                agTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCg                                                                                                                 <  1:1272146/150‑1 (MQ=255)
                                      gCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTa                                                                                                           >  2:3400750/1‑150 (MQ=255)
                                        cAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc                                                                                                         >  2:2897639/1‑150 (MQ=255)
                                        cAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc                                                                                                         <  1:1644342/150‑1 (MQ=255)
                                         aGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACa                                                                                                        >  1:1660211/1‑150 (MQ=255)
                                         aGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACa                                                                                                        <  1:70800/150‑1 (MQ=255)
                                           cTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACAtt                                                                                                      <  2:2266005/150‑1 (MQ=255)
                                             aaTCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTAt                                                                                                    <  2:1346782/150‑1 (MQ=255)
                                              atcatcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTAtc                                                                                                   >  1:1869273/1‑150 (MQ=255)
                                                 atcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTcc                                                                                                <  1:3215129/150‑1 (MQ=255)
                                                 atcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTcc                                                                                                >  1:1833710/1‑150 (MQ=255)
                                                   cGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAc                                                                                              <  1:1337788/150‑1 (MQ=255)
                                                    gCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAcc                                                                                             >  2:4333000/1‑150 (MQ=255)
                                                    gCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGAGCTAACATTATCTCCAca                                                                                             >  2:3007758/1‑149 (MQ=255)
                                                    gCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAcc                                                                                             <  1:4013638/150‑1 (MQ=255)
                                                    gCCTACCAGAACGGCGCGCCAAGTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAcc                                                                                             <  1:3030592/150‑1 (MQ=255)
                                                     ccTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAccg                                                                                            >  1:3781929/1‑150 (MQ=255)
                                                      cTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAccgc                                                                                           >  2:1644843/1‑150 (MQ=255)
                                                       tACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAccgcc                                                                                          <  2:1869273/150‑1 (MQ=255)
                                                       tACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAccgcc                                                                                          >  2:4000974/1‑150 (MQ=255)
                                                         acAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGa                                                                                        >  1:1783284/2‑150 (MQ=255)
                                                           aGAACGGCGCGCCACTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACa                                                                                      <  2:3935573/150‑1 (MQ=255)
                                                           aGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACa                                                                                      <  1:2554235/150‑1 (MQ=255)
                                                                ggCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCAtt                                                                                 <  2:243307/150‑1 (MQ=255)
                                                                ggCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCAtt                                                                                 >  1:1115678/1‑150 (MQ=255)
                                                                  cgcgcCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCg                                                                               <  2:2594382/150‑1 (MQ=255)
                                                                  cgcgcCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCg                                                                               >  1:2906556/1‑150 (MQ=255)
                                                                  cgcgcCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCg                                                                               >  1:272878/1‑150 (MQ=255)
                                                                      ccAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAg                                                                           <  1:3649742/150‑1 (MQ=255)
                                                                        aaTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAt                                                                         <  1:3557253/150‑1 (MQ=255)
                                                                        aaTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAt                                                                         >  2:1841020/1‑150 (MQ=255)
                                                                        aaTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAt                                                                         <  1:2283076/150‑1 (MQ=255)
                                                                         aTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtg                                                                        >  2:2690032/1‑150 (MQ=255)
                                                                         aTTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtg                                                                        >  2:1203787/1‑150 (MQ=255)
                                                                           tCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtgct                                                                      <  1:346532/150‑1 (MQ=255)
                                                                           tCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtgct                                                                      <  1:3898515/150‑1 (MQ=255)
                                                                           tCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtgct                                                                      <  1:1709838/150‑1 (MQ=255)
                                                                             gTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGAtgctgc                                                                    >  1:3732813/1‑150 (MQ=255)
                                                                               cTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGACAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCca                                                                  <  2:3846175/150‑1 (MQ=255)
                                                                               cTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCca                                                                  <  2:1719297/150‑1 (MQ=255)
                                                                                 gggCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCcaca                                                                <  1:3487089/150‑1 (MQ=255)
                                                                                  ggCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAg                                                               <  2:1004297/150‑1 (MQ=255)
                                                                                     gATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTc                                                            >  2:1103138/1‑150 (MQ=255)
                                                                                     gATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTc                                                            <  2:3081784/150‑1 (MQ=255)
                                                                                     gATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTc                                                            <  2:234503/150‑1 (MQ=255)
                                                                                      aTGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCa                                                           <  1:1655280/150‑1 (MQ=255)
                                                                                       tGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCAc                                                          >  1:140265/1‑150 (MQ=255)
                                                                                        gTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACt                                                         <  1:1189642/150‑1 (MQ=255)
                                                                                          cGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTga                                                       >  2:2045522/1‑150 (MQ=255)
                                                                                            cAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTgaga                                                     >  1:448530/1‑150 (MQ=255)
                                                                                            cAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTgaga                                                     >  2:2743410/1‑150 (MQ=255)
                                                                                             aaCTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGagag                                                    >  2:79120/1‑150 (MQ=255)
                                                                                             aaCTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGagag                                                    <  1:3593390/150‑1 (MQ=255)
                                                                                             aaCTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGagag                                                    <  1:2446388/150‑1 (MQ=255)
                                                                                              aCTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGt                                                   >  1:3787359/1‑150 (MQ=255)
                                                                                              aCTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGt                                                   >  1:1514846/1‑150 (MQ=255)
                                                                                                tGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCt                                                 >  2:3354515/1‑150 (MQ=255)
                                                                                                        aaGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAAt                                         >  2:2525798/1‑150 (MQ=255)
                                                                                                        aaGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAAt                                         >  1:51371/1‑150 (MQ=255)
                                                                                                            tGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGt                                     <  1:3251549/150‑1 (MQ=255)
                                                                                                             gCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTg                                    <  2:3883631/150‑1 (MQ=255)
                                                                                                             gCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTg                                    <  1:25081/150‑1 (MQ=255)
                                                                                                               aGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGaa                                  >  1:2016016/1‑150 (MQ=255)
                                                                                                                gAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCAGCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAg                                 <  2:859492/150‑1 (MQ=255)
                                                                                                                      aGCTGGCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACa                           >  1:2808313/1‑150 (MQ=255)
                                                                                                                          ggCTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGc                       >  1:2984352/1‑150 (MQ=255)
                                                                                                                            cTCGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCtt                     <  2:1692731/150‑1 (MQ=255)
                                                                                                                              cGGCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTc                   >  2:2365236/1‑150 (MQ=255)
                                                                                                                               ggCGCGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTcc                  >  1:2352494/1‑150 (MQ=255)
                                                                                                                                gcgcGTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCg                 <  1:1565553/150‑1 (MQ=255)
                                                                                                                                   cgTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATc              >  1:3083916/1‑150 (MQ=255)
                                                                                                                                   cgTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATc              <  2:4274408/150‑1 (MQ=255)
                                                                                                                                    gTGGGCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCg             >  1:1932607/1‑150 (MQ=255)
                                                                                                                                      gggCGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGca           >  1:2578303/1‑150 (MQ=255)
                                                                                                                                       ggCGAACAAAGTACAGGCCAGTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGcac          <  2:734255/150‑1 (MQ=255)
                                                                                                                                         cGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCTGATCGcacca        >  1:3423061/1‑150 (MQ=255)
                                                                                                                                         cGAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGcacca        <  2:3539239/150‑1 (MQ=255)
                                                                                                                                          gAACAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGcaccac       >  2:597022/1‑150 (MQ=255)
                                                                                                                                           aaCAAAGTACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGCaccacc      <  2:476928/150‑1 (MQ=255)
                                                                                                                                               aaGTACAGGCCATTGTGATGACTGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGCACCACCAata  >  2:439000/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTCCACCGCCGACAGCATTCGGCAGATGCTGCCACAGCTCACTGAGAGTCTGCCGAAATCGGTGAAGGTGACAGTGCTTTCCGATCGCACCACCAATA  >  CP009273/2149457‑2149749

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: