Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,604,776 A→T 6.9% 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 6.9% 1021.2 / 48.1 333S238C (AGT→TGT) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base A (130/180);  new base T (9/14);  total (139/194)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 341 total aligned reads displayed.

GTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGC  >  CP009273/2604627‑2604925
                                                                                                                                                     |                                                                                                                                                     
gTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATa                                                                                                                                                       >  2:1950845/1‑150 (MQ=255)
gTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATa                                                                                                                                                       >  2:128634/1‑150 (MQ=255)
  aTCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGt                                                                                                                                                     <  1:4737221/150‑1 (MQ=255)
   tCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTa                                                                                                                                                    <  1:4564185/150‑1 (MQ=255)
     gCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATc                                                                                                                                                  <  2:2386790/150‑1 (MQ=255)
     gCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATc                                                                                                                                                  <  1:1547754/150‑1 (MQ=255)
     gCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATc                                                                                                                                                  <  1:693517/150‑1 (MQ=255)
      ccATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCg                                                                                                                                                 >  2:1305217/1‑150 (MQ=255)
      ccATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCg                                                                                                                                                 <  2:4710635/150‑1 (MQ=255)
       cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                                                                >  1:396507/1‑150 (MQ=255)
       cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                                                                >  2:4192915/1‑150 (MQ=255)
       cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                                                                >  1:2149791/1‑150 (MQ=255)
       cATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGt                                                                                                                                                >  1:620573/1‑150 (MQ=255)
          tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                                                             <  2:2813393/150‑1 (MQ=255)
          tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                                                             <  2:3411708/150‑1 (MQ=255)
          tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                                                             <  1:2537531/150‑1 (MQ=255)
               gagaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAc                                                                                                                                        >  2:228606/1‑150 (MQ=255)
                    gcgcGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCttt                                                                                                                                   <  1:987648/150‑1 (MQ=255)
                      gcggcgGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                                                                                 <  2:2125624/150‑1 (MQ=255)
                      gcggcgGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTcg                                                                                                                                 >  2:3284997/1‑150 (MQ=255)
                        ggcggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTcgcg                                                                                                                               >  1:2991045/1‑150 (MQ=255)
                          cggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTa                                                                                                                             <  1:1718716/150‑1 (MQ=255)
                           ggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTAc                                                                                                                            >  2:2692137/1‑150 (MQ=255)
                           ggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTAc                                                                                                                            <  2:4885684/150‑1 (MQ=255)
                            gcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACt                                                                                                                           <  1:1599956/150‑1 (MQ=255)
                              tGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTcc                                                                                                                         >  2:4269883/1‑150 (MQ=255)
                                   ggCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCggggg                                                                                                                    >  2:2518013/1‑150 (MQ=255)
                                    gCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGt                                                                                                                   <  2:1731594/150‑1 (MQ=255)
                                          gTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATg                                                                                                             <  1:3037598/150‑1 (MQ=255)
                                          gTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATg                                                                                                             <  2:4808836/150‑1 (MQ=255)
                                             aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                                                                          >  2:3089258/1‑150 (MQ=255)
                                               gggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGa                                                                                                        <  1:4269275/150‑1 (MQ=255)
                                               gggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGa                                                                                                        <  1:461462/150‑1 (MQ=255)
                                                ggCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGaa                                                                                                       >  2:4451907/1‑150 (MQ=255)
                                                  cGAGAAATGTGCGATCGGCTCCTGCGCCATTTCAGCGCCGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAaga                                                                                                     <  2:3565992/150‑1 (MQ=255)
                                                   gagaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                                                                    <  2:3558620/150‑1 (MQ=255)
                                                   gagaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAagag                                                                                                    >  1:1785887/1‑150 (MQ=255)
                                                    tgaATTTTGCGATCGGCTCCTGCGCCATGTCAGCGTCGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGt                                                                                                   <  1:3006315/149‑1 (MQ=255)
                                                       aaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCg                                                                                                <  2:3656973/150‑1 (MQ=255)
                                                          ttGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCgggg                                                                                             <  2:3854789/150‑1 (MQ=255)
                                                             cGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGaaa                                                                                          >  1:3404364/1‑150 (MQ=255)
                                                                tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTg                                                                                       <  2:1407243/150‑1 (MQ=255)
                                                                tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTg                                                                                       <  1:3331465/150‑1 (MQ=255)
                                                                 cGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGTGGTTACATGACTGAAGAGTGCGGGGAAATTGt                                                                                      >  2:2335275/1‑150 (MQ=255)
                                                                 cGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGt                                                                                      >  2:180916/1‑150 (MQ=255)
                                                                 cGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGt                                                                                      <  1:3539671/150‑1 (MQ=255)
                                                                    cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                                                                   >  1:1979017/1‑150 (MQ=255)
                                                                    cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                                                                   >  1:3831808/1‑150 (MQ=255)
                                                                    cTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGtttt                                                                                   <  1:444194/150‑1 (MQ=255)
                                                                          cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGAc                                                                             <  1:5036180/150‑1 (MQ=255)
                                                                          cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGACATTGTTTTCTGGAc                                                                             >  1:1132325/1‑150 (MQ=255)
                                                                          cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGAc                                                                             <  1:2510786/150‑1 (MQ=255)
                                                                           gcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACg                                                                            <  2:1566076/150‑1 (MQ=255)
                                                                             cATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCt                                                                          <  2:4359036/150‑1 (MQ=255)
                                                                             cATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCt                                                                          <  1:415591/150‑1 (MQ=255)
                                                                                  aaGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaa                                                                     >  2:1133887/1‑150 (MQ=255)
                                                                                    gCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaa                                                                   >  2:1647477/1‑150 (MQ=255)
                                                                                     cgacgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaaa                                                                  >  1:1013726/1‑150 (MQ=255)
                                                                                     cgacgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTATCTGGACGCTGCGaaaaa                                                                  >  1:1284217/1‑150 (MQ=255)
                                                                                       acgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAg                                                                <  2:3260839/150‑1 (MQ=255)
                                                                                         gaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGtt                                                              >  1:268302/1‑150 (MQ=255)
                                                                                           tCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTg                                                            <  1:4118517/150‑1 (MQ=255)
                                                                                             cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc                                                          >  2:4008965/1‑150 (MQ=255)
                                                                                               aCAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGc                                                        >  1:2407436/1‑150 (MQ=255)
                                                                                                   gAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGt                                                    <  1:4196560/150‑1 (MQ=255)
                                                                                                     aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg                                                  >  1:876635/1‑150 (MQ=255)
                                                                                                      aTCGGTTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCATCTCTTTCTCGTACTCCGGGTGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGc                                                 <  1:3598602/150‑1 (MQ=255)
                                                                                                       tCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcg                                                <  2:1875309/150‑1 (MQ=255)
                                                                                                         gggTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTAAATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcgac                                              >  2:2227741/1‑150 (MQ=255)
                                                                                                           gTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTCGcgacga                                            >  2:4744340/1‑150 (MQ=255)
                                                                                                              aaCACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGAt                                         <  2:51415/150‑1 (MQ=255)
                                                                                                               acacCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATg                                        <  2:3987523/150‑1 (MQ=255)
                                                                                                                cacCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGc                                       >  2:184764/1‑150 (MQ=255)
                                                                                                                     ggTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGaa                                  >  1:2219180/1‑150 (MQ=255)
                                                                                                                       ttGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAaca                                <  2:687320/150‑1 (MQ=255)
                                                                                                                       ttGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAaca                                <  1:1497739/150‑1 (MQ=255)
                                                                                                                        tGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAacac                               <  1:1950845/150‑1 (MQ=255)
                                                                                                                          cgcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAg                             >  2:295888/1‑150 (MQ=255)
                                                                                                                          cgcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAg                             <  1:3714921/150‑1 (MQ=255)
                                                                                                                           gcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGc                            <  2:207547/150‑1 (MQ=255)
                                                                                                                           gcgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGc                            >  1:5304184/1‑150 (MQ=255)
                                                                                                                            cgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGct                           <  2:3356044/150‑1 (MQ=255)
                                                                                                                            cgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGct                           <  1:2071663/150‑1 (MQ=255)
                                                                                                                            cgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGct                           <  2:4610401/150‑1 (MQ=255)
                                                                                                                            cgATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGct                           <  1:3284997/150‑1 (MQ=255)
                                                                                                                               tGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCtctc                        <  1:1064212/150‑1 (MQ=255)
                                                                                                                                 ccGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAg                      <  2:2494197/150‑1 (MQ=255)
                                                                                                                                   gacgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGt                    >  1:2230704/1‑150 (MQ=255)
                                                                                                                                     cgaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTaa                  <  2:3208398/150‑1 (MQ=255)
                                                                                                                                      gaTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAAt                 <  1:1305217/150‑1 (MQ=255)
                                                                                                                                          cACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTAt             <  2:2159030/150‑1 (MQ=255)
                                                                                                                                           aCGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTAtt            <  2:5329407/150‑1 (MQ=255)
                                                                                                                                            cGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTa           <  1:4944423/150‑1 (MQ=255)
                                                                                                                                             gTCTCAATTGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTAc          <  2:3503598/150‑1 (MQ=255)
                                                                                                                                                 cAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgc      >  2:4582861/1‑150 (MQ=255)
                                                                                                                                                   aTAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCtcta    >  2:4081770/1‑150 (MQ=255)
                                                                                                                                                   aTAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgcta    <  2:3631508/150‑1 (MQ=255)
                                                                                                                                                     aGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgctagc  >  1:5171033/1‑150 (MQ=255)
                                                                                                                                                     aGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgctagc  >  1:2807261/1‑150 (MQ=255)
                                                                                                                                                     aGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCgctagc  >  2:2447225/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGACGAGATGCCGGAACACAGCTCTCAGGTAATGTATTACTCGCTAGC  >  CP009273/2604627‑2604925

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: