Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,602,418 C→A 8.8% I191I (ATC→ATA)  hyfG → hydrogenase 4, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,602,4180CA 8.8% 1050.5 / 55.4 284I191I (ATC→ATA) hyfGhydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base C (114/145);  new base A (12/13);  total (126/158)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 287 total aligned reads displayed.

GATGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACCG  >  CP009273/2602269‑2602567
                                                                                                                                                     |                                                                                                                                                     
 aTGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCa                                                                                                                                                      <  1:3956706/150‑1 (MQ=255)
   gACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCAcc                                                                                                                                                    <  1:776832/150‑1 (MQ=255)
     cTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTc                                                                                                                                                  >  2:4515321/1‑150 (MQ=255)
      tGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTcc                                                                                                                                                 <  2:787063/150‑1 (MQ=255)
      tGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTcc                                                                                                                                                 >  2:1460566/1‑150 (MQ=255)
      tGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTcc                                                                                                                                                 <  1:3273859/150‑1 (MQ=255)
        gCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGa                                                                                                                                               <  2:3530898/150‑1 (MQ=255)
         ccGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGAt                                                                                                                                              <  1:3240731/150‑1 (MQ=255)
          cGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATg                                                                                                                                             >  1:306501/1‑150 (MQ=255)
           ggAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGa                                                                                                                                            >  1:1373525/1‑150 (MQ=255)
              aGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAAcc                                                                                                                                         <  1:4107694/150‑1 (MQ=255)
               gATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCg                                                                                                                                        <  2:4320002/150‑1 (MQ=255)
                  atGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGt                                                                                                                                     <  2:3114035/150‑1 (MQ=255)
                   tGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTc                                                                                                                                    <  1:4326805/150‑1 (MQ=255)
                   tGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTc                                                                                                                                    <  2:1249446/150‑1 (MQ=255)
                     cATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCAc                                                                                                                                  >  1:2685194/1‑150 (MQ=255)
                       tCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACtt                                                                                                                                <  2:4257588/150‑1 (MQ=255)
                        ccGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTc                                                                                                                               <  2:1889252/150‑1 (MQ=255)
                            tGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCt                                                                                                                           >  1:3535479/1‑150 (MQ=255)
                              cgcAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGAAGACCCGGGTCACTTCCACttg                                                                                                                         >  2:1773313/1‑150 (MQ=255)
                                 aaaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgttt                                                                                                                      >  2:2720856/1‑150 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCACTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     <  1:2984851/150‑1 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     <  2:699648/150‑1 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     <  1:3234605/150‑1 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     <  2:2821291/150‑1 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     >  2:2570967/1‑150 (MQ=255)
                                       gCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGAt                                                                                                                <  1:725458/150‑1 (MQ=255)
                                        cGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATg                                                                                                               <  2:1090292/150‑1 (MQ=255)
                                         gATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATgg                                                                                                              <  1:192556/150‑1 (MQ=255)
                                         gATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATgg                                                                                                              <  2:3937318/150‑1 (MQ=255)
                                         gATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATgg                                                                                                              <  1:1681009/150‑1 (MQ=255)
                                            ggATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGa                                                                                                           <  1:3796339/150‑1 (MQ=255)
                                            ggATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGa                                                                                                           >  1:3484881/1‑150 (MQ=255)
                                              aTTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGc                                                                                                         >  1:3384054/1‑150 (MQ=255)
                                                  tCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAt                                                                                                     >  1:985781/1‑150 (MQ=255)
                                                  tCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAt                                                                                                     >  1:4047150/1‑150 (MQ=255)
                                                   cGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACGGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAtt                                                                                                    <  1:938664/150‑1 (MQ=255)
                                                   cGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAtt                                                                                                    <  1:1687296/150‑1 (MQ=255)
                                                   cGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAtt                                                                                                    <  2:3831982/150‑1 (MQ=255)
                                                   cGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAtt                                                                                                    <  2:18841/150‑1 (MQ=255)
                                                    gACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTg                                                                                                   >  2:2239233/1‑150 (MQ=255)
                                                         cgcCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGAt                                                                                              <  1:656111/150‑1 (MQ=255)
                                                          gcCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATg                                                                                             <  2:1373525/150‑1 (MQ=255)
                                                          gcCCTGAACCGACGACTGATTCCGAAACATATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATg                                                                                             >  1:4611996/1‑150 (MQ=255)
                                                           cccTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTGTGTGGATGGCGAGCAAATTGTCGATGc                                                                                            >  1:3697113/1‑150 (MQ=255)
                                                            ccTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCt                                                                                           <  2:2093512/150‑1 (MQ=255)
                                                              tGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGa                                                                                         <  2:2750438/150‑1 (MQ=255)
                                                                aaCCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGAtt                                                                                       <  2:3822892/150‑1 (MQ=255)
                                                                aaCCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGAtt                                                                                       <  2:1632620/150‑1 (MQ=255)
                                                                 aCCTACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTa                                                                                      <  1:2222516/150‑1 (MQ=255)
                                                                 aCCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTa                                                                                      <  2:4387629/150‑1 (MQ=255)
                                                                   cgacgaCTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTAcc                                                                                    <  2:4110239/150‑1 (MQ=255)
                                                                     acgacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGc                                                                                  <  2:2983378/150‑1 (MQ=255)
                                                                      cgacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCACATTGTCGATGCTGATTACCGcc                                                                                 >  2:2614972/1‑150 (MQ=255)
                                                                      cgacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGcc                                                                                 >  1:3777070/1‑150 (MQ=255)
                                                                       gacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCt                                                                                <  2:4057410/150‑1 (MQ=255)
                                                                        acTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTg                                                                               <  2:2245221/150‑1 (MQ=255)
                                                                         cTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGt                                                                              >  1:1745559/1‑150 (MQ=255)
                                                                          tGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGtt                                                                             <  2:1313739/150‑1 (MQ=255)
                                                                          tGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGtt                                                                             >  2:1918607/1‑150 (MQ=255)
                                                                             ttCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTa                                                                          <  1:3037968/150‑1 (MQ=255)
                                                                                cGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGt                                                                       <  1:154813/150‑1 (MQ=255)
                                                                                 gAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTc                                                                      >  1:533774/1‑150 (MQ=255)
                                                                                  aaaCGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTcc                                                                     <  2:459800/150‑1 (MQ=255)
                                                                                  aaaCGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGGTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTcc                                                                     >  1:531574/1‑150 (MQ=255)
                                                                                        aTCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATcgcg                                                               <  2:3683785/150‑1 (MQ=255)
                                                                                           cGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCa                                                            >  2:4627117/1‑150 (MQ=255)
                                                                                            gTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCAt                                                           >  1:1312972/1‑150 (MQ=255)
                                                                                            gTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCAt                                                           <  1:2281731/150‑1 (MQ=255)
                                                                                             tttATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATg                                                          <  1:3561231/150‑1 (MQ=255)
                                                                                              ttATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATgg                                                         >  1:4246852/1‑150 (MQ=255)
                                                                                                aTCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgag                                                       >  1:467256/1‑150 (MQ=255)
                                                                                                 tCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgaga                                                      <  1:1043473/150‑1 (MQ=255)
                                                                                                 tCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgaga                                                      <  1:1043412/150‑1 (MQ=255)
                                                                                                  cAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGaa                                                     <  1:1499529/150‑1 (MQ=255)
                                                                                                   aaTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGaaa                                                    >  2:2006038/1‑150 (MQ=255)
                                                                                                   aaTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGaaa                                                    >  2:1683115/1‑150 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   >  1:2125238/1‑150 (MQ=255)
                                                                                                       aGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTgg                                                >  1:3672684/1‑150 (MQ=255)
                                                                                                        gggCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGc                                               >  1:1857887/1‑150 (MQ=255)
                                                                                                         ggCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCa                                              <  2:1558354/150‑1 (MQ=255)
                                                                                                         ggCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCAGATCACCTCCGATGAACCGGGTCACTTCCGCTTGATTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCa                                              <  2:2396834/150‑1 (MQ=255)
                                                                                                           cAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGa                                            <  1:2845203/150‑1 (MQ=255)
                                                                                                            aaCAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGaa                                           >  1:2762515/1‑150 (MQ=255)
                                                                                                            aaCAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGaa                                           >  2:2249354/1‑150 (MQ=255)
                                                                                                               aGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAcg                                        >  2:3180371/1‑150 (MQ=255)
                                                                                                                 cGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAcgcg                                      >  1:4485489/1‑150 (MQ=255)
                                                                                                                         gggTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCt                              <  1:4145090/150‑1 (MQ=255)
                                                                                                                          ggTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTa                             <  2:985674/150‑1 (MQ=255)
                                                                                                                            tGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACa                           >  2:2414437/1‑150 (MQ=255)
                                                                                                                            tGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACa                           <  2:3089205/150‑1 (MQ=255)
                                                                                                                            tGATCCCAGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGAGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACa                           >  2:3790822/1‑150 (MQ=255)
                                                                                                                                     tCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGa                  >  2:4694761/1‑150 (MQ=255)
                                                                                                                                       ggCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGAcc                >  1:864514/1‑150 (MQ=255)
                                                                                                                                            gCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACcttctt           <  2:519573/150‑1 (MQ=255)
                                                                                                                                             cTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGCTGGGCGACAACGAAGTGACCTTCTTa          >  1:2803050/1‑150 (MQ=255)
                                                                                                                                             cTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTa          >  2:2961549/1‑150 (MQ=255)
                                                                                                                                                    tCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGAcc   <  1:1662570/150‑1 (MQ=255)
                                                                                                                                                     cACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACcg  >  1:3006064/1‑150 (MQ=255)
                                                                                                                                                     cACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACcg  >  1:633878/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACCG  >  CP009273/2602269‑2602567

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: ‑