Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 479,244 C→A 62.2% A206S (GCC→TCC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273479,2440CA62.2% 216.2 / inf 323A206S (GCC→TCC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (58/64);  new base A (120/81);  total (178/145)
Fisher's exact test for biased strand distribution p-value = 3.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 324 total aligned reads displayed.

CTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAG  >  CP009273/479095‑479392
                                                                                                                                                     |                                                                                                                                                    
  gATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTTTTTGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAgg                                                                                                                                                    <  1:5399377/150‑1 (MQ=255)
  gATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAgg                                                                                                                                                    >  1:1362260/1‑150 (MQ=255)
   aTTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGt                                                                                                                                                   >  2:1223242/1‑150 (MQ=255)
      cACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAAt                                                                                                                                                >  1:1010169/1‑150 (MQ=255)
       aCTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATg                                                                                                                                               >  1:520979/1‑150 (MQ=255)
            cagcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATc                                                                                                                                          <  2:1587088/150‑1 (MQ=255)
            cagcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATc                                                                                                                                          <  2:3640667/150‑1 (MQ=255)
             agcagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCa                                                                                                                                         >  1:3158839/1‑150 (MQ=255)
                agGATTTTGCCGTACTCTTCAGTAGAGTTCTGAAGCGTCTGAGCAATAATAGAGGCGTTAAGATGTTGGCCTTTCACCGGCGGAGTACCACCGAGCTGACCCGCCGCAACCTGGTAGTTCTGCGCTTTGATGGCGGTAATGACATCAAcc                                                                                                                                      <  1:740551/150‑1 (MQ=255)
                 ggATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCg                                                                                                                                     <  1:1740856/150‑1 (MQ=255)
                   aTTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGc                                                                                                                                   >  2:2515977/1‑150 (MQ=255)
                     tttGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGt                                                                                                                                 >  1:5604259/1‑150 (MQ=255)
                      ttGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGtt                                                                                                                                <  2:282639/150‑1 (MQ=255)
                       tGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTa                                                                                                                               >  2:1512785/1‑150 (MQ=255)
                         ccGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGc                                                                                                                             >  1:5542117/1‑150 (MQ=255)
                          cGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCt                                                                                                                            <  1:3069598/150‑1 (MQ=255)
                            aaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTATCTgg                                                                                                                          >  2:1824280/1‑150 (MQ=255)
                            aaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTgg                                                                                                                          >  1:1700328/1‑150 (MQ=255)
                            aaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTgg                                                                                                                          >  1:5410577/1‑150 (MQ=255)
                              ctctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGAACCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGATAGCTGGaa                                                                                                                        >  2:5086166/1‑150 (MQ=255)
                               tctTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAAt                                                                                                                       >  1:1242269/1‑150 (MQ=255)
                                 ttCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAttt                                                                                                                     >  1:5287105/1‑150 (MQ=255)
                                 ttCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAttt                                                                                                                     <  1:2214098/150‑1 (MQ=255)
                                   cAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt                                                                                                                   >  1:569297/1‑150 (MQ=255)
                                         agGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                                                                                             <  2:742711/150‑1 (MQ=255)
                                         agGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                                                                                             <  2:4775256/150‑1 (MQ=255)
                                         agGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCt                                                                                                             <  2:2726108/150‑1 (MQ=255)
                                             cAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCAtt                                                                                                         <  2:833036/150‑1 (MQ=255)
                                              aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTc                                                                                                        <  1:5154701/150‑1 (MQ=255)
                                                aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                                                                                                      <  1:1821763/150‑1 (MQ=255)
                                                 cgcgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCggg                                                                                                     >  2:571341/1‑150 (MQ=255)
                                                    gTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTc                                                                                                  >  2:3567492/1‑150 (MQ=255)
                                                    gTCTGAGCAAGAAGAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTc                                                                                                  <  2:3862949/150‑1 (MQ=255)
                                                     tCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCa                                                                                                 <  2:5131196/150‑1 (MQ=255)
                                                      cTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCAt                                                                                                >  2:1407378/1‑150 (MQ=255)
                                                      cTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAAAGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCAt                                                                                                >  2:2306314/1‑150 (MQ=255)
                                                       tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                               <  2:2572435/150‑1 (MQ=255)
                                                       tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                               >  1:323512/1‑150 (MQ=255)
                                                       tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                               >  1:2131524/1‑150 (MQ=255)
                                                       tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                               >  1:1158987/1‑150 (MQ=255)
                                                        gAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATcc                                                                                              <  1:4230665/150‑1 (MQ=255)
                                                         aGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCa                                                                                             >  2:4951432/1‑150 (MQ=255)
                                                         aGCAATAAGAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCa                                                                                             <  1:662929/150‑1 (MQ=255)
                                                          gCAATCATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAg                                                                                            <  2:3973076/150‑1 (MQ=255)
                                                           cAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGa                                                                                           >  1:3287847/1‑150 (MQ=255)
                                                            aataatAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGAt                                                                                          >  2:4549279/1‑150 (MQ=255)
                                                             ataataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATa                                                                                         >  1:2424448/1‑150 (MQ=255)
                                                               aataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACg                                                                                       >  2:551472/1‑150 (MQ=255)
                                                               aataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACg                                                                                       <  1:359389/150‑1 (MQ=255)
                                                                  agagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCAt                                                                                    >  1:3633620/1‑150 (MQ=255)
                                                                    agGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATcg                                                                                  <  2:1675157/150‑1 (MQ=255)
                                                                    agGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATcg                                                                                  <  2:4376632/150‑1 (MQ=255)
                                                                         ttAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTAc                                                                             <  2:692570/150‑1 (MQ=255)
                                                                          tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                            <  1:4301181/150‑1 (MQ=255)
                                                                          tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                            <  1:4138281/150‑1 (MQ=255)
                                                                          tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                            >  2:4692501/1‑150 (MQ=255)
                                                                            aGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgt                                                                          <  1:4548145/150‑1 (MQ=255)
                                                                              cTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                        >  2:666251/1‑150 (MQ=255)
                                                                              cTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                        <  1:1019318/150‑1 (MQ=255)
                                                                              cTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                        >  1:1959889/1‑150 (MQ=255)
                                                                               tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGaa                                                                       >  1:2187302/1‑150 (MQ=255)
                                                                               tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGaa                                                                       >  1:4644322/1‑150 (MQ=255)
                                                                                 ttGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAAcc                                                                     >  2:450246/1‑150 (MQ=255)
                                                                                  tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                    <  1:2247152/150‑1 (MQ=255)
                                                                                  tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                    >  1:2407913/1‑150 (MQ=255)
                                                                                  tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                    >  2:3378775/1‑150 (MQ=255)
                                                                                   ggCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGa                                                                   >  1:4097026/1‑150 (MQ=255)
                                                                                    gCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaa                                                                  >  1:3480555/1‑150 (MQ=255)
                                                                                    gCCTGTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaa                                                                  <  2:5482981/150‑1 (MQ=255)
                                                                                     ccTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaac                                                                 >  2:5601185/1‑150 (MQ=255)
                                                                                        ttCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaacaac                                                              >  2:2211039/1‑150 (MQ=255)
                                                                                        ttCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaacaac                                                              >  1:4351042/1‑150 (MQ=255)
                                                                                            ccGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGaa                                                          >  1:1182310/1‑150 (MQ=255)
                                                                                            ccGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGATAGCTGGAATTTGTTCAGCTAATACGGGTTCAACCAGATAAGCATCGCGTACTTTGAACCGAACAACTGaa                                                          >  1:2979839/1‑150 (MQ=255)
                                                                                              ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACa                                                        <  2:237475/150‑1 (MQ=255)
                                                                                                 ggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCa                                                     <  2:1010169/150‑1 (MQ=255)
                                                                                                 ggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCa                                                     >  2:3238296/1‑150 (MQ=255)
                                                                                                      accaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCAc                                                >  2:5233892/1‑150 (MQ=255)
                                                                                                        caccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGc                                              >  2:1271170/1‑150 (MQ=255)
                                                                                                             aGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGAc                                         <  1:1436533/150‑1 (MQ=255)
                                                                                                              gCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACg                                        <  2:227492/150‑1 (MQ=255)
                                                                                                                              cTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATc                        >  2:1158608/1‑150 (MQ=255)
                                                                                                                                 ggCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCttt                     <  1:1194603/150‑1 (MQ=255)
                                                                                                                                 ggCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCttt                     >  2:4290837/1‑150 (MQ=255)
                                                                                                                                 ggCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCttt                     >  2:4089044/1‑150 (MQ=255)
                                                                                                                                  gCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTc                    <  2:4130101/150‑1 (MQ=255)
                                                                                                                                   cGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCa                   >  2:2556338/1‑150 (MQ=255)
                                                                                                                                   cGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCa                   <  2:1885532/150‑1 (MQ=255)
                                                                                                                                     ttCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCata                 >  1:303646/1‑150 (MQ=255)
                                                                                                                                       cTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATAtt               >  2:1107510/1‑150 (MQ=255)
                                                                                                                                        tGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTc              >  1:2966954/1‑150 (MQ=255)
                                                                                                                                          cgcTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgc            <  1:3831915/150‑1 (MQ=255)
                                                                                                                                           gcTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgcc           >  1:4905070/1‑150 (MQ=255)
                                                                                                                                            cTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgccg          >  2:5283513/1‑150 (MQ=255)
                                                                                                                                              ttGATGGCGGTAATGACATCAACCGGCGTGAGATGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGACTGATGGCATCTTTCATATTCgccgcc        <  2:2980724/150‑1 (MQ=255)
                                                                                                                                               tGATGGATGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCa       <  1:2665893/150‑1 (MQ=255)
                                                                                                                                                gATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCAc      <  1:2567263/150‑1 (MQ=255)
                                                                                                                                                gATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCAc      >  1:4610597/1‑150 (MQ=255)
                                                                                                                                                  tGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGt    >  1:398011/1‑150 (MQ=255)
                                                                                                                                                  tGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGt    >  1:1538332/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAG  >  CP009273/479095‑479392

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: