Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,604,776 A→T 8.3% S238C (AGT→TGT)  hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,604,7760AT 8.3% 990.3 / 58.6 322S238C (AGT→TGT) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base A (124/171);  new base T (11/16);  total (135/187)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 324 total aligned reads displayed.

GTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAG  >  CP009273/2604627‑2604924
                                                                                                                                                     |                                                                                                                                                    
gTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATa                                                                                                                                                      <  2:3532947/150‑1 (MQ=255)
   tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                                                                   >  2:2385835/1‑150 (MQ=255)
   tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                                                                   <  1:1840722/150‑1 (MQ=255)
    cGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTAt                                                                                                                                                  >  1:4042584/1‑150 (MQ=255)
      ccATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCg                                                                                                                                                <  1:1875006/150‑1 (MQ=255)
       cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                                                               >  1:10131/1‑150 (MQ=255)
            aGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGc                                                                                                                                          >  2:1843438/1‑150 (MQ=255)
            aGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGc                                                                                                                                          >  1:4724869/1‑150 (MQ=255)
             gCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCa                                                                                                                                         <  2:5031347/150‑1 (MQ=255)
              cGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAACTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCaa                                                                                                                                        >  1:1332912/1‑150 (MQ=255)
              cGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCaa                                                                                                                                        <  2:1308512/150‑1 (MQ=255)
              cGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCaa                                                                                                                                        <  1:3985276/150‑1 (MQ=255)
               gagaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAc                                                                                                                                       >  1:4538042/1‑150 (MQ=255)
                agaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAct                                                                                                                                      <  2:4336234/150‑1 (MQ=255)
                agaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAct                                                                                                                                      <  1:810486/150‑1 (MQ=255)
                   aGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACCGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                                                                   >  2:4433735/1‑150 (MQ=255)
                   aGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                                                                   >  2:4555922/1‑150 (MQ=255)
                   aGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                                                                   <  1:2516000/150‑1 (MQ=255)
                     cgcgGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTc                                                                                                                                 >  2:2835047/1‑150 (MQ=255)
                      gcggcgGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                                                                                >  1:3332698/1‑150 (MQ=255)
                        ggcggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcgcg                                                                                                                              <  2:4013313/150‑1 (MQ=255)
                         gcggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGt                                                                                                                             <  2:3282862/150‑1 (MQ=255)
                           ggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTAc                                                                                                                           <  1:2295024/150‑1 (MQ=255)
                            gcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACt                                                                                                                          <  2:831977/150‑1 (MQ=255)
                              tGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTcc                                                                                                                        <  2:4796665/150‑1 (MQ=255)
                                gCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCgg                                                                                                                      <  2:2149894/150‑1 (MQ=255)
                                    gCTATCGTCAGGGGCGAGAAATTGGCGCTCGGCTCCTGCGCCATTTCAGCGACGATCCTCCAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGGTCAGCACCTCTTTCGCGTACTCCGGGGGt                                                                                                                  <  2:2034781/150‑1 (MQ=255)
                                      tATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTa                                                                                                                <  1:1604160/150‑1 (MQ=255)
                                       aTCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTAc                                                                                                               <  1:2818370/150‑1 (MQ=255)
                                       aTCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTAc                                                                                                               >  1:3884601/1‑150 (MQ=255)
                                         cGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACAt                                                                                                             >  1:4088749/1‑150 (MQ=255)
                                               gggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGa                                                                                                       >  2:1231113/1‑150 (MQ=255)
                                                ggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGaa                                                                                                      <  2:4571724/150‑1 (MQ=255)
                                                ggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCGCAATAGTATCGTTCAGCAACTCTTACGCGTACTCCGGGGGTTACATGACTGaa                                                                                                      >  2:3491744/1‑150 (MQ=255)
                                                  cGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAaga                                                                                                    <  2:364281/150‑1 (MQ=255)
                                                  cGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAaga                                                                                                    <  2:521698/150‑1 (MQ=255)
                                                    agaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGt                                                                                                  <  2:4247739/150‑1 (MQ=255)
                                                     gaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTg                                                                                                 <  2:842898/150‑1 (MQ=255)
                                                       aaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCg                                                                                               <  1:4126623/150‑1 (MQ=255)
                                                        aTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgg                                                                                              <  2:1935862/150‑1 (MQ=255)
                                                           tGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGa                                                                                           >  1:247554/1‑150 (MQ=255)
                                                           tGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGa                                                                                           >  1:456714/1‑150 (MQ=255)
                                                            gCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGaa                                                                                          <  1:2203556/150‑1 (MQ=255)
                                                              gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                                                                        <  2:102714/150‑1 (MQ=255)
                                                              gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                                                                        >  1:712158/1‑150 (MQ=255)
                                                              gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                                                                        >  1:2825585/1‑150 (MQ=255)
                                                               aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                                                                       <  2:4487417/150‑1 (MQ=255)
                                                               aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                                                                       <  1:2722034/150‑1 (MQ=255)
                                                               aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                                                                       <  1:1661718/150‑1 (MQ=255)
                                                               aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                                                                       <  2:2613441/150‑1 (MQ=255)
                                                                 cGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGt                                                                                     <  1:3159515/150‑1 (MQ=255)
                                                                   gCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGttt                                                                                   <  1:2007817/150‑1 (MQ=255)
                                                                    cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                                                                  >  2:3272305/1‑150 (MQ=255)
                                                                     tCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGGAGAGTGCGGGGAAATTGTTTTc                                                                                 <  2:1225039/150‑1 (MQ=255)
                                                                      ccTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCt                                                                                <  2:861509/150‑1 (MQ=255)
                                                                       cTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTg                                                                               <  2:3131559/150‑1 (MQ=255)
                                                                       cTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTg                                                                               <  1:65598/150‑1 (MQ=255)
                                                                        tGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTgg                                                                              <  1:246254/150‑1 (MQ=255)
                                                                        tGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTgg                                                                              >  2:4208650/1‑150 (MQ=255)
                                                                          cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGAc                                                                            >  2:3892056/1‑150 (MQ=255)
                                                                           gcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACg                                                                           <  2:4856791/150‑1 (MQ=255)
                                                                            ccATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGc                                                                          >  1:1501792/1‑150 (MQ=255)
                                                                                        cgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGt                                                              >  2:3498669/1‑150 (MQ=255)
                                                                                         gaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGtt                                                             >  1:1738022/1‑150 (MQ=255)
                                                                                            ccTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGt                                                          >  1:837093/1‑150 (MQ=255)
                                                                                            ccTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGt                                                          <  1:4968095/150‑1 (MQ=255)
                                                                                            ccTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGt                                                          >  1:2091956/1‑150 (MQ=255)
                                                                                             cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc                                                         <  2:94236/150‑1 (MQ=255)
                                                                                             cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc                                                         <  1:3564802/150‑1 (MQ=255)
                                                                                               aCAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGc                                                       >  1:3263946/1‑150 (MQ=255)
                                                                                               aCAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGc                                                       >  1:1425691/1‑150 (MQ=255)
                                                                                                 aGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCa                                                     <  2:2787325/150‑1 (MQ=255)
                                                                                                  ggAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAg                                                    >  2:3405412/1‑150 (MQ=255)
                                                                                                     aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg                                                 >  1:1167209/1‑150 (MQ=255)
                                                                                                     aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg                                                 >  1:4874840/1‑150 (MQ=255)
                                                                                                       tCGGGTGAAGACCTGGGTGCGCGATGCCGCCGTTCCACGTCTCAATTGTATGGTTCAGCATCTCTTTCGCGTACTCCGGGGGTTACATGACTGATGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCGAGTAGcg                                               <  2:3428867/150‑1 (MQ=255)
                                                                                                        cGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcga                                              >  1:344031/1‑150 (MQ=255)
                                                                                                        cGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcga                                              <  1:1090087/150‑1 (MQ=255)
                                                                                                             tAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACgaga                                         >  1:4696290/1‑150 (MQ=255)
                                                                                                                    tGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGa                                  >  2:1685677/1‑150 (MQ=255)
                                                                                                                      gTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAac                                <  1:2432946/150‑1 (MQ=255)
                                                                                                                       gtGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAaca                               <  2:63850/149‑1 (MQ=255)
                                                                                                                         gcgcGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAAcaca                             <  2:1533998/150‑1 (MQ=255)
                                                                                                                            cgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGct                          >  1:91683/1‑150 (MQ=255)
                                                                                                                              aTGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGctct                        <  1:2900507/150‑1 (MQ=255)
                                                                                                                              aTGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGctct                        <  1:3188408/150‑1 (MQ=255)
                                                                                                                                gCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCa                      <  1:1531521/150‑1 (MQ=255)
                                                                                                                                 ccGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAg                     <  2:4106079/150‑1 (MQ=255)
                                                                                                                                  cgacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAgg                    <  1:4232825/150‑1 (MQ=255)
                                                                                                                                  cgacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAATAGCGACGAGATGCCGGAACACAGCTCTCAgg                    >  2:3802903/1‑150 (MQ=255)
                                                                                                                                   gacgaTCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGt                   >  1:3810296/1‑150 (MQ=255)
                                                                                                                                      gaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAAt                <  1:4407319/150‑1 (MQ=255)
                                                                                                                                       aTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATg               >  2:2882842/1‑150 (MQ=255)
                                                                                                                                        tCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGt              <  1:4839156/150‑1 (MQ=255)
                                                                                                                                         ccACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTa             <  1:1981068/150‑1 (MQ=255)
                                                                                                                                          cACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTACTGTAt            >  2:4294292/1‑150 (MQ=255)
                                                                                                                                          cACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTAt            <  1:3310136/150‑1 (MQ=255)
                                                                                                                                             gTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTAc         >  1:3393925/1‑150 (MQ=255)
                                                                                                                                                  aaTAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgct    >  1:3843631/1‑150 (MQ=255)
                                                                                                                                                  aaTAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTAGTACTCgct    >  1:5071056/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAG  >  CP009273/2604627‑2604924

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: