Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 479,244 C→A 85.4% 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,2440CA85.4% 694.6 / 124.4 289A206S (GCC→TCC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (16/26);  new base A (130/115);  total (148/141)
Fisher's exact test for biased strand distribution p-value = 9.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 289 total aligned reads displayed.

GATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGT  >  CP009273/479097‑479390
                                                                                                                                                   |                                                                                                                                                  
gATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAgg                                                                                                                                                  >  1:1177369/1‑150 (MQ=255)
gATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAgg                                                                                                                                                  >  2:1881573/1‑150 (MQ=255)
 aTTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGt                                                                                                                                                 >  2:4530078/1‑150 (MQ=255)
 aTTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGt                                                                                                                                                 >  2:243749/1‑150 (MQ=255)
   tCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTaa                                                                                                                                               <  2:1014938/150‑1 (MQ=255)
    cACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAAt                                                                                                                                              <  1:2505643/150‑1 (MQ=255)
     aCTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATg                                                                                                                                             >  1:1767960/1‑150 (MQ=255)
      cTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCCGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTTGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCCCCGCAACCTGGGCGTTCTGCGCTGTGATGGCGGGAAAGa                                                                                                                                            >  1:28250/1‑150 (MQ=255)
      cTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCCGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCAGGGCGGTCTGCGCGTTGATGGAGGTAATGa                                                                                                                                            >  1:765071/1‑150 (MQ=255)
        ttCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACa                                                                                                                                          <  2:816394/150‑1 (MQ=255)
             cagGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAAc                                                                                                                                     >  2:2321502/1‑150 (MQ=255)
              agGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAAcc                                                                                                                                    >  2:4374869/1‑150 (MQ=255)
               ggATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCg                                                                                                                                   >  2:3505020/1‑150 (MQ=255)
                gATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCgg                                                                                                                                  <  1:1869430/150‑1 (MQ=255)
                gATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCgg                                                                                                                                  >  2:2117403/1‑150 (MQ=255)
                 aTTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGc                                                                                                                                 >  2:2816227/1‑150 (MQ=255)
                 aTTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGc                                                                                                                                 <  2:1313582/150‑1 (MQ=255)
                 aTTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAAGAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGc                                                                                                                                 <  2:1313365/150‑1 (MQ=255)
                     tGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTa                                                                                                                             >  1:190911/1‑150 (MQ=255)
                     tGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTa                                                                                                                             >  2:671675/1‑150 (MQ=255)
                       ccGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGc                                                                                                                           >  2:3449458/1‑150 (MQ=255)
                         gAACTCTTCAGTAGAGGTCAGACACGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTg                                                                                                                         <  1:4028382/150‑1 (MQ=255)
                          aaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATTACATAAACCCGCGTTAGCTgg                                                                                                                        >  1:1976382/1‑150 (MQ=255)
                          aaCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTgg                                                                                                                        >  1:2985316/1‑150 (MQ=255)
                           aCTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGa                                                                                                                       >  2:2961864/1‑150 (MQ=255)
                               ttCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAAttt                                                                                                                   <  1:305150/150‑1 (MQ=255)
                                 cAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCCACCGGCGTTAGCTGGAATTTGt                                                                                                                 >  2:1108906/1‑150 (MQ=255)
                                 cAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGt                                                                                                                 >  2:4429032/1‑150 (MQ=255)
                                   gTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCGGGGCGTTCTGCGCTTTTATGGAGGTAATGACATCAACCGGCGTTTGCTGGAAtttttta                                                                                                               >  1:479161/1‑149 (MQ=255)
                                    tAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCa                                                                                                              <  1:3715630/150‑1 (MQ=255)
                                     agagGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAg                                                                                                             <  1:3244344/150‑1 (MQ=255)
                                     agagGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAg                                                                                                             >  2:3321102/1‑150 (MQ=255)
                                        ggTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTc                                                                                                          >  1:3349433/1‑150 (MQ=255)
                                         gTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCa                                                                                                         >  2:2187164/1‑150 (MQ=255)
                                           cAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCAtt                                                                                                       <  1:2546938/150‑1 (MQ=255)
                                           cAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCAtt                                                                                                       <  1:2264568/150‑1 (MQ=255)
                                            aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTTTTCAGCTCATTc                                                                                                      >  1:3650145/1‑150 (MQ=255)
                                            aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTc                                                                                                      >  1:3698785/1‑150 (MQ=255)
                                            aGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAAACGGCGTTAGCTGGAATTTTTTCAGCTAATTc                                                                                                      >  1:103711/1‑150 (MQ=255)
                                              aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCgg                                                                                                    <  1:1443331/150‑1 (MQ=255)
                                                 cgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCCGCGGCGTACCACCGAGCTGACCCGCCCCAACCAGGGCGTTGTGCGCGTTTATGGAGGTAAATACAACAACCCGCGGTAGCTGGAAAATTTTAAGCACAATAGGGtt                                                                                                 >  1:3420249/1‑150 (MQ=255)
                                                     tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                             <  1:3240544/150‑1 (MQ=255)
                                                     tGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATc                                                                                             >  2:4063854/1‑150 (MQ=255)
                                                      gAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATcc                                                                                            >  2:1197187/1‑150 (MQ=255)
                                                       aGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCa                                                                                           <  1:4135023/150‑1 (MQ=255)
                                                        gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAg                                                                                          <  2:1076144/150‑1 (MQ=255)
                                                        gCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAg                                                                                          <  2:1925002/150‑1 (MQ=255)
                                                         cAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGa                                                                                         >  2:1778988/1‑150 (MQ=255)
                                                          aataatAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGAt                                                                                        <  2:3339908/150‑1 (MQ=255)
                                                          aataatAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGAt                                                                                        <  2:2370549/150‑1 (MQ=255)
                                                             aataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACg                                                                                     <  2:1390732/150‑1 (MQ=255)
                                                                agagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCAt                                                                                  <  1:3854135/150‑1 (MQ=255)
                                                                 gagGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATc                                                                                 <  2:492534/150‑1 (MQ=255)
                                                                     cGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt                                                                             >  2:2295005/1‑150 (MQ=255)
                                                                      gTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTa                                                                            >  2:2924456/1‑150 (MQ=255)
                                                                        tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                          >  2:2604854/1‑150 (MQ=255)
                                                                        tAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCACAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                          <  2:3222053/150‑1 (MQ=255)
                                                                        tAAGCTGTTGGCCTTTCACCGGCGGCATACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACt                                                                          <  1:189588/150‑1 (MQ=255)
                                                                          aGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgt                                                                        >  2:361407/1‑150 (MQ=255)
                                                                          aGCTGTTGGACTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACtgt                                                                        <  1:3146844/150‑1 (MQ=255)
                                                                            cTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                      >  1:1006283/1‑150 (MQ=255)
                                                                                tGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                  >  2:884023/1‑150 (MQ=255)
                                                                                   ccTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaac                                                               >  2:1159979/1‑150 (MQ=255)
                                                                                    cTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaaca                                                              >  2:1941958/1‑150 (MQ=255)
                                                                                       tCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACt                                                           >  2:2632586/1‑150 (MQ=255)
                                                                                        cACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTg                                                          >  1:2185936/1‑150 (MQ=255)
                                                                                         aCCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGa                                                         >  2:3966522/1‑150 (MQ=255)
                                                                                         aCCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGa                                                         >  1:4506886/1‑150 (MQ=255)
                                                                                         aCCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGa                                                         <  2:2738910/150‑1 (MQ=255)
                                                                                         aCCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGGTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGa                                                         <  2:2688762/150‑1 (MQ=255)
                                                                                            ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACa                                                      >  1:2672856/1‑150 (MQ=255)
                                                                                              cggcgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACaac                                                    >  1:4010029/1‑148 (MQ=255)
                                                                                                 cgTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCAcc                                                 >  2:3371461/1‑150 (MQ=255)
                                                                                                  gTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGGACTGTGAACCGAACAACTGAACATCAccc                                                >  1:849202/1‑150 (MQ=255)
                                                                                                    accaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCAc                                              <  2:2418683/150‑1 (MQ=255)
                                                                                                    accaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCAc                                              <  2:600371/150‑1 (MQ=255)
                                                                                                        ccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGccc                                          >  1:1366309/1‑150 (MQ=255)
                                                                                                        ccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGccc                                          <  1:671675/150‑1 (MQ=255)
                                                                                                           aGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGAc                                       >  2:3723235/1‑150 (MQ=255)
                                                                                                               gACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTAc                                   <  1:4312753/150‑1 (MQ=255)
                                                                                                                   cgccgcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCt                               <  1:1881573/150‑1 (MQ=255)
                                                                                                                     ccgcAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGa                             >  1:1613856/1‑150 (MQ=255)
                                                                                                                        cAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                          <  1:1097420/150‑1 (MQ=255)
                                                                                                                        cAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                          <  1:2428829/150‑1 (MQ=255)
                                                                                                                        cAACCTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                          <  2:1612557/150‑1 (MQ=255)
                                                                                                                           ccTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCAt                       >  1:678835/1‑150 (MQ=255)
                                                                                                                            cTGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATc                      >  1:618186/1‑150 (MQ=255)
                                                                                                                             tGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCt                     <  1:2083292/150‑1 (MQ=255)
                                                                                                                             tGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCt                     <  1:3440680/150‑1 (MQ=255)
                                                                                                                             tGGGCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCGACGTACGGCTGATGGCATCtt                    >  2:113940/1‑150 (MQ=255)
                                                                                                                              gggCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCtt                    <  1:3478705/150‑1 (MQ=255)
                                                                                                                               ggCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCttt                   <  2:4226820/150‑1 (MQ=255)
                                                                                                                               ggCGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCttt                   <  1:2471310/150‑1 (MQ=255)
                                                                                                                                 cGTTCTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTc                  <  1:4449676/149‑1 (MQ=255)
                                                                                                                                     cTGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATAtt             <  2:3033/150‑1 (MQ=255)
                                                                                                                                      tGCGCTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTc            <  2:3283838/150‑1 (MQ=255)
                                                                                                                                        cgcTTTGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgc          >  2:4196302/1‑150 (MQ=255)
                                                                                                                                             tGATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCa     >  2:1695751/1‑150 (MQ=255)
                                                                                                                                              gATGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCAc    >  2:39009/1‑150 (MQ=255)
                                                                                                                                               aTGGAGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACg   <  1:2105895/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATTCACTTTCAGCAGGATTTTGCCGAACTCTTCAGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGT  >  CP009273/479097‑479390

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: