Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,604,731 C→T 15.8% R223W (CGG→TGG)  hyfI → hydrogenase 4, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,604,7310CT15.8% 953.7 / 137.0 337R223W (CGG→TGG) hyfIhydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (148/135);  new base T (19/34);  total (167/170)
Fisher's exact test for biased strand distribution p-value = 3.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

GGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGA  >  CP009273/2604583‑2604880
                                                                                                                                                    |                                                                                                                                                     
 gTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTgg                                                                                                                                                     <  2:4435174/150‑1 (MQ=255)
   gACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGt                                                                                                                                                   >  2:4573743/1‑150 (MQ=255)
   gACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGt                                                                                                                                                   >  2:1300474/1‑150 (MQ=255)
    aCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGtt                                                                                                                                                  >  1:4766567/1‑150 (MQ=255)
      tATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTaa                                                                                                                                                >  2:3233910/1‑150 (MQ=255)
      tATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTaa                                                                                                                                                >  2:2023102/1‑150 (MQ=255)
      tATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTaa                                                                                                                                                >  2:4198725/1‑150 (MQ=255)
       aTGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAac                                                                                                                                               <  2:2358267/150‑1 (MQ=255)
        tGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAaca                                                                                                                                              <  2:1441488/150‑1 (MQ=255)
         gCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATGGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAacac                                                                                                                                             <  1:4605848/150‑1 (MQ=255)
            cccGTTGTGGCCGCATATCCCGCCAAAACCGCGTATCGCCATTTAGCGAGATGCGCGGGGGCTGGCGGGCTATCTTCAGGGGCGTGAAATTTGCTATGGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTg                                                                                                                                          <  1:4451226/149‑1 (MQ=255)
            aCCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTg                                                                                                                                          >  2:1781201/1‑150 (MQ=255)
            aCCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTg                                                                                                                                          <  1:303273/150‑1 (MQ=255)
            aCCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTg                                                                                                                                          <  1:1426284/150‑1 (MQ=255)
               gTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGtt                                                                                                                                       <  2:1767480/150‑1 (MQ=255)
                ttGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTg                                                                                                                                      <  2:875146/150‑1 (MQ=255)
                 tgtgGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTgc                                                                                                                                     >  1:3574831/1‑150 (MQ=255)
                      ccGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGAt                                                                                                                                <  2:194580/150‑1 (MQ=255)
                       cGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATg                                                                                                                               >  1:528099/1‑150 (MQ=255)
                        gCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGc                                                                                                                              <  1:4681376/150‑1 (MQ=255)
                         cAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGcc                                                                                                                             >  2:3517381/1‑150 (MQ=255)
                          aGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcg                                                                                                                            >  2:719816/1‑150 (MQ=255)
                          aGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcg                                                                                                                            >  2:2818810/1‑150 (MQ=255)
                          aGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGGTGCGCGATGCcg                                                                                                                            >  2:3813448/1‑150 (MQ=255)
                           gATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcga                                                                                                                           >  1:6614/1‑150 (MQ=255)
                           gATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcga                                                                                                                           <  2:3461765/150‑1 (MQ=255)
                             tCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCcgacg                                                                                                                         >  2:1561623/1‑150 (MQ=255)
                               ccGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGAt                                                                                                                       >  1:1250574/1‑150 (MQ=255)
                                   cATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCAc                                                                                                                   >  1:4379732/1‑150 (MQ=255)
                                     tCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGt                                                                                                                 >  2:3125090/1‑150 (MQ=255)
                                      cacaGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGtc                                                                                                                <  1:4259722/150‑1 (MQ=255)
                                      cacaGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGtc                                                                                                                >  1:1003845/1‑150 (MQ=255)
                                            gTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATa                                                                                                          <  1:797545/150‑1 (MQ=255)
                                             tATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAg                                                                                                         <  1:1450638/150‑1 (MQ=255)
                                              aTCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGt                                                                                                        <  1:1125623/150‑1 (MQ=255)
                                                cGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTAt                                                                                                      <  2:1918407/150‑1 (MQ=255)
                                                cGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTAt                                                                                                      <  1:4303698/150‑1 (MQ=255)
                                                  ccATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCg                                                                                                    >  1:2707226/1‑150 (MQ=255)
                                                    aTTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                                  >  2:1352512/1‑150 (MQ=255)
                                                    aTTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                                  >  1:3247780/1‑150 (MQ=255)
                                                    aTTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGtt                                                                                                  >  1:2209276/1‑150 (MQ=255)
                                                      tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCTGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                <  1:2562498/150‑1 (MQ=255)
                                                      tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                >  2:133095/1‑150 (MQ=255)
                                                      tGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCa                                                                                                >  2:2494279/1‑150 (MQ=255)
                                                           gagaAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAAc                                                                                           >  1:2136387/1‑150 (MQ=255)
                                                               aGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCtt                                                                                       <  1:2769217/150‑1 (MQ=255)
                                                                gcgcGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCttt                                                                                      <  2:1463176/150‑1 (MQ=255)
                                                                      cggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTa                                                                                <  1:1846044/150‑1 (MQ=255)
                                                                       ggcTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTAc                                                                               >  1:3464441/1‑150 (MQ=255)
                                                                         cTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTc                                                                             >  2:4288454/1‑150 (MQ=255)
                                                                           ggCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCg                                                                           <  2:2593814/150‑1 (MQ=255)
                                                                                  tATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTa                                                                    <  2:72808/150‑1 (MQ=255)
                                                                                   aTCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTAc                                                                   <  1:1793665/150‑1 (MQ=255)
                                                                                      gTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATg                                                                >  1:164444/1‑150 (MQ=255)
                                                                                       tCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGa                                                               <  1:4184258/150‑1 (MQ=255)
                                                                                        cAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGAc                                                              >  1:3914776/1‑150 (MQ=255)
                                                                                        cAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGAc                                                              >  2:845793/1‑150 (MQ=255)
                                                                                         aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCCATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                             >  2:2718669/1‑150 (MQ=255)
                                                                                         aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                             <  1:1703289/150‑1 (MQ=255)
                                                                                         aGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACt                                                             <  2:2707226/150‑1 (MQ=255)
                                                                                                agaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGt                                                      <  1:3291476/150‑1 (MQ=255)
                                                                                                 gaAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGta                                                     <  2:4053108/150‑2 (MQ=255)
                                                                                                  aaaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTGCTCCGGGGGTTACATGACTGAAGAGTGc                                                    <  1:2602554/150‑1 (MQ=255)
                                                                                                  aaaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGc                                                    <  1:2500692/150‑1 (MQ=255)
                                                                                                   aaTTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCg                                                   >  1:2341114/1‑150 (MQ=255)
                                                                                                     tttGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCggg                                                 >  1:315524/1‑150 (MQ=255)
                                                                                                         cGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGaaa                                             >  2:3299609/1‑150 (MQ=255)
                                                                                                          gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATTGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                            <  1:132549/150‑1 (MQ=255)
                                                                                                          gATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAt                                            >  2:4317050/1‑150 (MQ=255)
                                                                                                           aTCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCCGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAAtt                                           <  1:149600/150‑1 (MQ=255)
                                                                                                            tCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTg                                          <  1:3092762/150‑1 (MQ=255)
                                                                                                                 tCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTc                                     >  2:168969/1‑150 (MQ=255)
                                                                                                                  ccTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCt                                    >  2:4744818/1‑150 (MQ=255)
                                                                                                                     gcgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGa                                 <  2:219647/150‑1 (MQ=255)
                                                                                                                      cgcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGAc                                >  2:4151997/1‑150 (MQ=255)
                                                                                                                       gcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACg                               <  2:927961/150‑1 (MQ=255)
                                                                                                                       gcCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACg                               >  1:4755748/1‑150 (MQ=255)
                                                                                                                          aTTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGACGAGTGCGGGGAAATTGTTGTCTGGCCGCTg                            >  2:2178852/1‑150 (MQ=255)
                                                                                                                          aTTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTg                            >  2:1520393/1‑150 (MQ=255)
                                                                                                                            ttAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCg                          <  2:3902583/150‑1 (MQ=255)
                                                                                                                             tAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGa                         >  1:2244987/1‑150 (MQ=255)
                                                                                                                                  gacgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGaaaaaa                    >  2:3286058/1‑150 (MQ=255)
                                                                                                                                    cgaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGt                  >  1:978002/1‑150 (MQ=255)
                                                                                                                                     gaTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGtt                 >  1:939493/1‑150 (MQ=255)
                                                                                                                                      aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt                >  1:1114081/1‑150 (MQ=255)
                                                                                                                                      aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt                <  1:3775777/150‑1 (MQ=255)
                                                                                                                                      aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt                >  2:3487636/1‑150 (MQ=255)
                                                                                                                                      aTCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGttt                >  2:2918764/1‑150 (MQ=255)
                                                                                                                                       tCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGACATTGTTTTCTGGACGCTGCGAAAAAAGTTTg               >  2:3100632/1‑150 (MQ=255)
                                                                                                                                         cTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGCTCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTc             <  2:3697007/150‑1 (MQ=255)
                                                                                                                                          tACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCg            <  1:3846200/150‑1 (MQ=255)
                                                                                                                                          tACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGAGCCACGTCTCAAGAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCg            <  1:2178852/150‑1 (MQ=255)
                                                                                                                                            cAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGcc          >  2:3535652/1‑150 (MQ=255)
                                                                                                                                            cAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGcc          >  2:2448390/1‑150 (MQ=255)
                                                                                                                                                aaaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTa      <  2:3401057/150‑1 (MQ=255)
                                                                                                                                                aaaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTa      >  1:1604552/1‑150 (MQ=255)
                                                                                                                                                 aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg     <  2:1003845/150‑1 (MQ=255)
                                                                                                                                                 aaTCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGAACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAg     >  2:1180799/1‑150 (MQ=255)
                                                                                                                                                   tCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcg   <  2:1708908/150‑1 (MQ=255)
                                                                                                                                                   tCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGcg   <  1:1920363/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCGGCTCCTGCGCCATTTAAGCGACGATCCTACAGGAAATCGGGTTAACACCTGGTTGCGCGATGCCGACGATCCACGTCTCAATAGTATCGTTCAGCAACTCTTTCGCGTACTCCGGGGGTTACATGACTGAAGAGTGCGGGGAAATTGTTTTCTGGACGCTGCGAAAAAAGTTTGTCGCCAGTAGCGA  >  CP009273/2604583‑2604880

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: