Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,602,418 C→A 19.1% 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,4180CA19.1% 629.0 / 107.0 226I191I (ATC→ATA) hyfGhydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base C (81/101);  new base A (20/24);  total (101/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

GATGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACCG  >  CP009273/2602269‑2602567
                                                                                                                                                     |                                                                                                                                                     
gattccTGGGCGGAATATCTGCATCCGCTGCGCAAAGATGCGTTGGATTATCGACTTCGCCCTGAACCGACGACTGCTTCCGAAACGTATCCGTTTATCATTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATc                                                                                                                                                       <  1:1956616/145‑1 (MQ=255)
  tGACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCAc                                                                                                                                                     <  2:3836225/150‑1 (MQ=255)
   gACTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCAcc                                                                                                                                                    <  1:612434/150‑1 (MQ=255)
     cTGGCCGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTc                                                                                                                                                  <  1:1785526/150‑1 (MQ=255)
         ccGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGAt                                                                                                                                              >  1:976907/1‑150 (MQ=255)
         ccGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGAt                                                                                                                                              >  2:2556099/1‑150 (MQ=255)
          cGGAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCAGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATg                                                                                                                                             <  2:3401257/150‑1 (MQ=255)
           ggAAGATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGa                                                                                                                                            <  2:3262468/150‑1 (MQ=255)
              agagATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAAcc                                                                                                                                         <  2:3935295/150‑1 (MQ=255)
               gATATGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCg                                                                                                                                        <  2:639328/150‑1 (MQ=255)
                atatGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCgg                                                                                                                                       <  2:4451197/150‑1 (MQ=255)
                  atGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGt                                                                                                                                     >  1:833901/1‑150 (MQ=255)
                  atGCATCCGCTGCGCAAAGAAGCGATGGATTATCGACTGCGCACTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGAGCGGGTGATCCCTGTCGGCCCGCTGCAAATCACCTCCGATGAACCGGGt                                                                                                                                     >  1:1454231/1‑150 (MQ=255)
                   tGCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTc                                                                                                                                    <  1:3529588/150‑1 (MQ=255)
                    gCATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCa                                                                                                                                   <  1:2100220/150‑1 (MQ=255)
                     cATCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCAc                                                                                                                                  >  1:4036458/1‑150 (MQ=255)
                       tCCGCTGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACtt                                                                                                                                <  2:2265358/150‑1 (MQ=255)
                            tGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCt                                                                                                                           <  1:3644209/150‑1 (MQ=255)
                            tGCGCAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCt                                                                                                                           <  2:15339/150‑1 (MQ=255)
                             gcgcAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCtt                                                                                                                          <  1:1343680/150‑1 (MQ=255)
                             gcgcAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCtt                                                                                                                          <  2:4813611/150‑1 (MQ=255)
                               gcAAAGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGACCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgt                                                                                                                        >  1:2294239/1‑150 (MQ=255)
                                  aaGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttg                                                                                                                     <  1:1704304/150‑1 (MQ=255)
                                   aGATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCttgtttgt                                                                                                                    <  2:2423295/150‑1 (MQ=255)
                                    gATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTtgtg                                                                                                                   >  2:3467629/1‑150 (MQ=255)
                                    gATGCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTtgtg                                                                                                                   >  2:2663152/1‑150 (MQ=255)
                                       gCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGAt                                                                                                                <  1:5152734/150‑1 (MQ=255)
                                       gCGATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGAt                                                                                                                <  1:1696816/150‑1 (MQ=255)
                                         gATGGATTATCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATgg                                                                                                              <  1:3991076/150‑1 (MQ=255)
                                                 aTCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCaac                                                                                                      >  2:1801298/1‑149 (MQ=255)
                                                  tCGACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAAt                                                                                                     >  2:4510516/1‑150 (MQ=255)
                                                    gACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTg                                                                                                   <  1:542035/150‑1 (MQ=255)
                                                    gACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTg                                                                                                   <  1:4501331/150‑1 (MQ=255)
                                                    gACTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTGTGTGGATGGCGAGCAAATTg                                                                                                   >  1:2249298/1‑150 (MQ=255)
                                                     aCTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGt                                                                                                  >  1:4226081/1‑150 (MQ=255)
                                                     aCTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTGGTTTGTGGATGGCGAGCAAATTGt                                                                                                  >  1:1631600/1‑150 (MQ=255)
                                                      cTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTc                                                                                                 >  2:3414982/1‑150 (MQ=255)
                                                      cTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTc                                                                                                 <  1:1529746/150‑1 (MQ=255)
                                                      cTGCGCCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTc                                                                                                 <  1:2612888/150‑1 (MQ=255)
                                                         cgcCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGAt                                                                                              <  2:3587361/150‑1 (MQ=255)
                                                          gcCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATg                                                                                             >  1:4466998/1‑150 (MQ=255)
                                                          gcCCTGAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATg                                                                                             <  1:17845/150‑1 (MQ=255)
                                                               gAACCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGAt                                                                                        <  1:4446882/150‑1 (MQ=255)
                                                                 aCCGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTa                                                                                      >  2:3335862/1‑150 (MQ=255)
                                                                  ccGACGACTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTAc                                                                                     >  1:2265403/1‑150 (MQ=255)
                                                                    gacgacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCg                                                                                   <  2:1850760/150‑1 (MQ=255)
                                                                       gacTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCt                                                                                >  2:4962528/1‑150 (MQ=255)
                                                                        acTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTg                                                                               >  2:15293/1‑150 (MQ=255)
                                                                         cTGATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGt                                                                              >  2:1635429/1‑150 (MQ=255)
                                                                           gATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGttt                                                                            >  1:1851648/1‑150 (MQ=255)
                                                                           gATTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGttt                                                                            >  2:3585010/1‑150 (MQ=255)
                                                                            aTTCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGtttt                                                                           >  1:4977768/1‑150 (MQ=255)
                                                                              tCCGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTAt                                                                         <  2:4658758/150‑1 (MQ=255)
                                                                                cGAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGt                                                                       <  2:4231317/150‑1 (MQ=255)
                                                                                 gAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTc                                                                      <  2:639877/150‑1 (MQ=255)
                                                                                 gAAACGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTc                                                                      >  2:4149559/1‑150 (MQ=255)
                                                                                   aaCGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCa                                                                    <  1:1222085/150‑1 (MQ=255)
                                                                                     cGTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATc                                                                  >  1:4877831/1‑150 (MQ=255)
                                                                                      gTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATcg                                                                 >  2:51134/1‑150 (MQ=255)
                                                                                      gTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATcg                                                                 <  2:3051214/150‑1 (MQ=255)
                                                                                      gTATCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATcg                                                                 >  2:310174/1‑150 (MQ=255)
                                                                                        aTCCGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATcgcg                                                               >  1:3230233/1‑150 (MQ=255)
                                                                                         tCCTTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCgg                                                              >  1:2218859/1‑150 (MQ=255)
                                                                                           cGTTTATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCa                                                            <  2:4626180/150‑1 (MQ=255)
                                                                                             tttATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTGTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTGTATGTCCATCGCGGCATg                                                          >  1:1993722/1‑150 (MQ=255)
                                                                                              ttATCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATgg                                                         <  1:4384957/150‑1 (MQ=255)
                                                                                                aTCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgag                                                       <  1:3131831/150‑1 (MQ=255)
                                                                                                aTCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgag                                                       >  2:3715940/1‑150 (MQ=255)
                                                                                                aTCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgag                                                       >  2:1117349/1‑150 (MQ=255)
                                                                                                 tCAATGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGgaga                                                      <  2:698633/150‑1 (MQ=255)
                                                                                                   aaTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGaaa                                                    >  2:2086596/1‑150 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCACATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   >  2:4141854/1‑150 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   >  2:280896/1‑150 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   <  2:751742/150‑1 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   >  2:2067944/1‑150 (MQ=255)
                                                                                                    aTGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAAc                                                   >  2:1065658/1‑150 (MQ=255)
                                                                                                     tGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACt                                                  <  1:3406169/150‑1 (MQ=255)
                                                                                                     tGAGGGCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACt                                                  <  1:4057805/150‑1 (MQ=255)
                                                                                                          gCAACAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAg                                             <  2:5695/150‑1 (MQ=255)
                                                                                                            aaCAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGaa                                           >  1:2851860/1‑150 (MQ=255)
                                                                                                              cAGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAc                                         <  2:4409146/150‑1 (MQ=255)
                                                                                                               aGCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAcg                                        >  1:4744105/1‑150 (MQ=255)
                                                                                                                gCGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAcgc                                       <  1:2663152/150‑1 (MQ=255)
                                                                                                                 cGATGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAAcgcg                                      <  1:2570943/150‑1 (MQ=255)
                                                                                                                   aTGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGa                                    >  1:4168791/1‑150 (MQ=255)
                                                                                                                    tGCGCGGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGAt                                   >  1:3821347/1‑150 (MQ=255)
                                                                                                                       gcgGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGACTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATggg                                <  1:3299913/150‑1 (MQ=255)
                                                                                                                        cgGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGc                               >  1:34223/1‑150 (MQ=255)
                                                                                                                        cgGGTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGc                               <  1:3600892/150‑1 (MQ=255)
                                                                                                                         gggTGATCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCt                              <  1:255370/150‑1 (MQ=255)
                                                                                                                               tCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACg                        <  2:833901/150‑1 (MQ=255)
                                                                                                                               tCCCTGTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACg                        <  2:1345362/150‑1 (MQ=255)
                                                                                                                                    gTCGGCCCGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTg                   >  2:31134/1‑150 (MQ=255)
                                                                                                                                    gTCGGCCCGCTGCATATAACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTg                   >  1:3751407/1‑150 (MQ=255)
                                                                                                                                         cccGCTGCATATCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACctt              <  1:1065818/150‑1 (MQ=255)
                                                                                                                                         cccGCTGCATATCACCTACGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACctt              <  1:1031019/150‑1 (MQ=255)
                                                                                                                                                  tataACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGa     >  2:3566372/1‑150 (MQ=255)
                                                                                                                                                  tataACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGa     >  2:3291259/1‑150 (MQ=255)
                                                                                                                                                   atCACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGAc    >  2:2579569/1‑150 (MQ=255)
                                                                                                                                                     cACCTCCGATGAACCGGGTCACTTCCGCTTGTTTGTGGATGGCGAGCAAATTGTCGATGCTGATTACCGCCTGTTTTATGTCCATCGCGGCATGGAGAAACTGGCAGAAACGCGGATGGGCTACAACGAAGTGACCTTCTTATCGGACcg  >  2:4196310/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: ‑