Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,460,619 T→A 7.9% F199Y (TTT→TAT)  intS → CPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,460,6190TA 7.9% 735.5 / 42.9 227F199Y (TTT→TAT) intSCPS‑53 (KpLE1) prophage; putative prophage CPS‑53 integrase
Reads supporting (aligned to +/- strand):  ref base T (94/115);  new base A (10/8);  total (104/123)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 236 total aligned reads displayed.

GATCGCGGTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCGATTTTG  >  CP009273/2460471‑2460768
                                                                                                                                                    |                                                                                                                                                     
gATCGCGGTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGttt                                                                                                                                                      <  2:2511843/150‑1 (MQ=255)
      ggTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCt                                                                                                                                                >  1:3800086/1‑150 (MQ=255)
       gTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTg                                                                                                                                               <  1:2363810/150‑1 (MQ=255)
        tGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGc                                                                                                                                              <  1:1691880/150‑1 (MQ=255)
         gCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCa                                                                                                                                             >  2:278317/1‑150 (MQ=255)
          cAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAg                                                                                                                                            <  2:889408/150‑1 (MQ=255)
          cAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAg                                                                                                                                            <  1:2042167/150‑1 (MQ=255)
           aaTGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGa                                                                                                                                           >  2:2273098/1‑150 (MQ=255)
             tGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGAcc                                                                                                                                         <  1:2503507/150‑1 (MQ=255)
              ggAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACACCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCa                                                                                                                                        <  1:2888873/150‑1 (MQ=255)
                aaCGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGa                                                                                                                                      <  1:2435499/150‑1 (MQ=255)
                aaCGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGa                                                                                                                                      >  2:4445522/1‑150 (MQ=255)
                 aCGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGAt                                                                                                                                     >  1:967053/1‑150 (MQ=255)
                   gAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATcc                                                                                                                                   >  1:1525655/1‑150 (MQ=255)
                    aGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATccc                                                                                                                                  >  2:2960276/1‑150 (MQ=255)
                    aGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATccc                                                                                                                                  <  1:254741/150‑1 (MQ=255)
                    aGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATccc                                                                                                                                  >  2:1933793/1‑150 (MQ=255)
                      ccAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCgg                                                                                                                                <  1:2030184/150‑1 (MQ=255)
                       caacaaAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc                                                                                                                               >  1:644212/1‑150 (MQ=255)
                       caacaaAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGc                                                                                                                               >  1:629411/1‑150 (MQ=255)
                           aaaGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTc                                                                                                                           <  2:773536/150‑1 (MQ=255)
                            aaGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTca                                                                                                                          <  2:3602032/150‑1 (MQ=255)
                              gCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTcaac                                                                                                                        <  2:4567046/150‑1 (MQ=255)
                                aCGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTcaacaa                                                                                                                      <  2:325984/150‑1 (MQ=255)
                                  gCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAg                                                                                                                    <  2:1099771/150‑1 (MQ=255)
                                   cAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGc                                                                                                                   <  2:2010075/150‑1 (MQ=255)
                                    agaagaTGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCa                                                                                                                  >  1:3266591/1‑150 (MQ=255)
                                    agaagaTGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCa                                                                                                                  >  1:2637936/1‑150 (MQ=255)
                                      aagaTGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACt                                                                                                                >  1:2359770/1‑150 (MQ=255)
                                       agaTGCGGCGAGGTTTTCCGTTACGCTATTNTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTg                                                                                                               >  1:4287225/1‑150 (MQ=255)
                                           gcggcgAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCaa                                                                                                           <  2:473316/150‑1 (MQ=255)
                                           gcggcgAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCaa                                                                                                           <  2:3575631/150‑1 (MQ=255)
                                              gcgAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACAt                                                                                                        <  2:4305769/150‑1 (MQ=255)
                                               cgAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACAtt                                                                                                       <  1:3585729/150‑1 (MQ=255)
                                                gAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACAttt                                                                                                      <  2:2424185/150‑1 (MQ=255)
                                                   gTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCa                                                                                                   <  1:1621739/150‑1 (MQ=255)
                                                       tCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGaa                                                                                               <  2:37728/150‑1 (MQ=255)
                                                         cGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgt                                                                                             <  1:1792219/150‑1 (MQ=255)
                                                         cGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgt                                                                                             >  1:3000321/1‑150 (MQ=255)
                                                           ttACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgtat                                                                                           <  1:545007/150‑1 (MQ=255)
                                                              cGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGTTGCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgtatcgt                                                                                        <  2:4064449/150‑1 (MQ=255)
                                                              cGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgtatcgt                                                                                        >  2:2955123/1‑150 (MQ=255)
                                                                cTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAgtatcgtat                                                                                      <  2:1738609/150‑1 (MQ=255)
                                                                  aTTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGtatcgtatcg                                                                                    <  2:1426442/150‑1 (MQ=255)
                                                                    tGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCt                                                                                  <  2:631447/150‑1 (MQ=255)
                                                                       cACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCAt                                                                               <  2:1217652/150‑1 (MQ=255)
                                                                       cACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCAt                                                                               <  1:2822939/150‑1 (MQ=255)
                                                                       cACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTCCAGGAAGTATCGTATCGCTCAt                                                                               <  1:1387096/150‑1 (MQ=255)
                                                                        aCCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCAtt                                                                              <  2:1389071/150‑1 (MQ=255)
                                                                           ggAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCg                                                                           >  1:495321/1‑150 (MQ=255)
                                                                           ggAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCg                                                                           >  1:976194/1‑150 (MQ=255)
                                                                           ggAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCg                                                                           <  2:2304445/150‑1 (MQ=255)
                                                                            gAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGa                                                                          >  2:592045/1‑150 (MQ=255)
                                                                             aaGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGAc                                                                         >  1:3332781/1‑150 (MQ=255)
                                                                              aGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGAcc                                                                        >  2:2431589/1‑150 (MQ=255)
                                                                                   tAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGt                                                                   <  1:1924082/150‑1 (MQ=255)
                                                                                         taACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACg                                                             >  2:2721675/1‑150 (MQ=255)
                                                                                         taACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACg                                                             >  1:4407693/1‑150 (MQ=255)
                                                                                           aCCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCt                                                           >  2:4373568/1‑150 (MQ=255)
                                                                                           aCCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCt                                                           <  1:3879348/150‑1 (MQ=255)
                                                                                                  aCCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGc                                                    >  1:394664/1‑150 (MQ=255)
                                                                                                   ccTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGcc                                                   >  2:3023809/1‑150 (MQ=255)
                                                                                                   ccTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGcc                                                   >  2:279933/1‑150 (MQ=255)
                                                                                                    cTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGccc                                                  >  2:3206414/1‑150 (MQ=255)
                                                                                                      gACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTa                                                <  2:1859112/150‑1 (MQ=255)
                                                                                                      gACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTa                                                <  1:2769462/150‑1 (MQ=255)
                                                                                                        ccTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAg                                              >  1:3583147/1‑150 (MQ=255)
                                                                                                         cTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGa                                             >  2:851788/1‑150 (MQ=255)
                                                                                                         cTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGa                                             <  2:3363630/150‑1 (MQ=255)
                                                                                                         cTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTGTACGCTACACAGCCCTAAGa                                             <  2:1163617/150‑1 (MQ=255)
                                                                                                            gCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACg                                          <  1:3880663/150‑1 (MQ=255)
                                                                                                               gACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGaaa                                       <  1:4532835/150‑1 (MQ=255)
                                                                                                                aCGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAag                                      <  1:3082174/150‑1 (MQ=255)
                                                                                                                 cGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAaga                                     <  1:3663836/150‑1 (MQ=255)
                                                                                                                  gCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAagag                                    <  1:2960276/150‑1 (MQ=255)
                                                                                                                   ccATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGc                                   <  1:1112149/150‑1 (MQ=255)
                                                                                                                   acATGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGc                                   <  1:2158817/149‑1 (MQ=255)
                                                                                                                     aTGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCtt                                 >  2:3819574/1‑150 (MQ=255)
                                                                                                                     aTGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCtt                                 <  2:2485054/150‑1 (MQ=255)
                                                                                                                     aTGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCtt                                 >  1:2016091/1‑150 (MQ=255)
                                                                                                                     aTGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCtt                                 >  2:3971232/1‑150 (MQ=255)
                                                                                                                     aTGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCtt                                 >  2:3971319/1‑150 (MQ=255)
                                                                                                                      tGAAGGGATACCGCAAGAAGAACTTCCCGTATCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTc                                >  2:2178319/1‑150 (MQ=255)
                                                                                                                         aGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGtt                             <  2:2637936/150‑1 (MQ=255)
                                                                                                                           ggATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTcc                           <  1:510700/150‑1 (MQ=255)
                                                                                                                              tACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATg                        >  1:2258949/1‑150 (MQ=255)
                                                                                                                                 cGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTa                     <  1:899342/150‑1 (MQ=255)
                                                                                                                                    aagaagAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATgg                  <  2:2375216/150‑1 (MQ=255)
                                                                                                                                    aagaagAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATgg                  >  2:1500328/1‑150 (MQ=255)
                                                                                                                                    aagaagAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATgg                  >  1:4073390/1‑150 (MQ=255)
                                                                                                                                    aagaagAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATgg                  >  1:3497059/1‑150 (MQ=255)
                                                                                                                                     agaagaACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGga                 >  1:1885236/1‑150 (MQ=255)
                                                                                                                                        agaaCTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGgaaga              >  2:3475376/1‑150 (MQ=255)
                                                                                                                                          aaCTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAAc            <  2:1582567/150‑1 (MQ=255)
                                                                                                                                           aCTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACg           <  1:3494038/150‑1 (MQ=255)
                                                                                                                                             ttCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTc         <  1:3462386/150‑1 (MQ=255)
                                                                                                                                              tCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCg        <  1:592045/150‑1 (MQ=255)
                                                                                                                                              tCCCGTTTCTTCCTGAAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCTGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAATCACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCg        <  1:439833/150‑1 (MQ=255)
                                                                                                                                                 cGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCGAtt     <  1:724943/150‑1 (MQ=255)
                                                                                                                                                   tttCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCGAtttt   <  1:1286255/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATCGCGGTGCAATGGAACGAGCCAACAAAGCACGCAGAAGATGCGGCGAGGTTTTCCGTTACGCTATTGTCACCGGAAGGGCTAAATATAACCCGGCACCTGACCTTGCTGACGCCATGAAGGGATACCGCAAGAAGAACTTCCCGTTTCTTCCTGCAGACCAGATCCCGGCATTCAACAAAGCACTGGCAACATTTTCAGGAAGTATCGTATCGCTCATTGCGACCAAAGTTTTACGCTACACAGCCCTAAGAACGAAAGAGCTTCGTTCCATGCTATGGAAGAACGTCGATTTTG  >  CP009273/2460471‑2460768

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: