Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,334 C→A 25.0% I141I (ATC→ATA yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,3340CA25.0% 1208.0 / inf 507I141I (ATC→ATAyagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (177/202);  new base A (68/60);  total (245/262)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

ATCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTGCC  >  CP009273/279185‑279483
                                                                                                                                                     |                                                                                                                                                     
   ggCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCACGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGATCGATTCCCTGCCGTACCGCAACGACGCGGCGATCgtg                                                                                                                                                    >  1:981358/1‑150 (MQ=255)
    gCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCgtgt                                                                                                                                                   <  1:3511098/150‑1 (MQ=255)
     cATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGTGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGtt                                                                                                                                                  >  2:4281258/1‑150 (MQ=255)
     cATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGtt                                                                                                                                                  >  2:1017402/1‑150 (MQ=255)
      atgcagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACAGCAACGACGCGGCGATCGTGttt                                                                                                                                                 >  1:3762406/1‑150 (MQ=255)
       tgcagatgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTc                                                                                                                                                <  2:433797/150‑1 (MQ=255)
            atgcagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTcgccgc                                                                                                                                           >  1:4309011/1‑150 (MQ=255)
               cagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTg                                                                                                                                        >  2:1153990/1‑150 (MQ=255)
               cagGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCACTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTg                                                                                                                                        >  1:3104825/1‑150 (MQ=255)
                agGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGa                                                                                                                                       >  2:2805574/1‑150 (MQ=255)
                   cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATcc                                                                                                                                    <  2:1431674/150‑1 (MQ=255)
                   cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATcc                                                                                                                                    <  1:3260272/150‑1 (MQ=255)
                   cggcggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATcc                                                                                                                                    <  2:2652301/150‑1 (MQ=255)
                      cggcgAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCt                                                                                                                                 >  1:506783/1‑150 (MQ=255)
                          gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                             <  2:4538178/150‑1 (MQ=255)
                          gAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCt                                                                                                                             >  1:2561083/1‑150 (MQ=255)
                           aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTg                                                                                                                            <  1:4666560/150‑1 (MQ=255)
                           aaGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTg                                                                                                                            <  1:2264939/150‑1 (MQ=255)
                            aGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGc                                                                                                                           >  2:3774688/1‑150 (MQ=255)
                             ggAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGcc                                                                                                                          >  1:4781188/1‑150 (MQ=255)
                             ggAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGcc                                                                                                                          <  2:881267/150‑1 (MQ=255)
                              gagaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCg                                                                                                                         >  2:1264877/1‑150 (MQ=255)
                               agaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGa                                                                                                                        >  2:4715144/1‑150 (MQ=255)
                               agaTCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGa                                                                                                                        >  2:2067978/1‑150 (MQ=255)
                                   cACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGAcgcg                                                                                                                    <  1:4609761/150‑1 (MQ=255)
                                    aCCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcgg                                                                                                                   >  1:3196752/1‑150 (MQ=255)
                                     cccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggc                                                                                                                  <  2:3349961/150‑1 (MQ=255)
                                      ccGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACgcggcg                                                                                                                 <  1:267470/150‑1 (MQ=255)
                                       cGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGcggcgg                                                                                                                >  1:1823653/1‑150 (MQ=255)
                                           aTGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGgcgggcgg                                                                                                            >  2:3910926/1‑150 (MQ=255)
                                              gCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGa                                                                                                         <  2:457389/150‑1 (MQ=255)
                                                gggATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATc                                                                                                       >  2:3032521/1‑150 (MQ=255)
                                                 ggATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCg                                                                                                      <  2:448148/150‑1 (MQ=255)
                                                    tCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCg                                                                                                   >  1:4358244/1‑150 (MQ=255)
                                                     cccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGt                                                                                                  >  1:3975603/1‑150 (MQ=255)
                                                     cccGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGt                                                                                                  <  2:1872309/150‑1 (MQ=255)
                                                        attCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGc                                                                                               <  1:794041/149‑1 (MQ=255)
                                                         ttCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCg                                                                                              <  1:3032530/150‑1 (MQ=255)
                                                         ttCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCg                                                                                              <  2:3714473/150‑1 (MQ=255)
                                                           cgcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAc                                                                                            <  2:3525561/150‑1 (MQ=255)
                                                           cgcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAc                                                                                            >  2:1598280/1‑150 (MQ=255)
                                                            gcgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGAcc                                                                                           >  1:283770/1‑150 (MQ=255)
                                                             cgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCt                                                                                          >  2:1075614/1‑150 (MQ=255)
                                                             cgGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGAGTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCt                                                                                          >  1:384520/1‑150 (MQ=255)
                                                              ggCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTg                                                                                         <  2:3115798/150‑1 (MQ=255)
                                                                ccTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCg                                                                                       >  2:1751404/1‑150 (MQ=255)
                                                                     gTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaaa                                                                                  >  1:1700844/1‑150 (MQ=255)
                                                                     gTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATaaa                                                                                  >  1:1551771/1‑150 (MQ=255)
                                                                        aGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAggg                                                                               >  1:922563/1‑150 (MQ=255)
                                                                          cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCt                                                                             >  2:1512197/1‑150 (MQ=255)
                                                                          cGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCt                                                                             <  1:4087647/150‑1 (MQ=255)
                                                                           gATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTg                                                                            <  1:1738779/150‑1 (MQ=255)
                                                                           gATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTg                                                                            <  2:703981/150‑1 (MQ=255)
                                                                            aTCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGc                                                                           <  1:3380797/150‑1 (MQ=255)
                                                                                  gCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGcca                                                                     <  1:2908402/150‑1 (MQ=255)
                                                                                  gCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGcca                                                                     <  1:11999/150‑1 (MQ=255)
                                                                                  gCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGcca                                                                     <  1:8653/150‑1 (MQ=255)
                                                                                   cGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccac                                                                    <  1:1486462/150‑1 (MQ=255)
                                                                                   cGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccac                                                                    >  2:1820592/1‑150 (MQ=255)
                                                                                     aCGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGccacca                                                                  <  2:504329/150‑1 (MQ=255)
                                                                                       gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                                >  1:2461681/1‑150 (MQ=255)
                                                                                       gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                                >  2:2278330/1‑150 (MQ=255)
                                                                                       gggCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATg                                                                >  2:2645782/1‑150 (MQ=255)
                                                                                            cTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgc                                                           <  2:982641/150‑1 (MQ=255)
                                                                                             tCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcg                                                          <  1:704974/150‑1 (MQ=255)
                                                                                             tCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcg                                                          <  1:4252803/150‑1 (MQ=255)
                                                                                             tCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTgcg                                                          >  1:2807353/1‑150 (MQ=255)
                                                                                               gCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCt                                                        >  2:777113/1‑150 (MQ=255)
                                                                                               gCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCt                                                        <  2:2480412/150‑1 (MQ=255)
                                                                                                 aGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTgg                                                      >  2:3839691/1‑150 (MQ=255)
                                                                                                 aGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTgg                                                      <  1:3686038/150‑1 (MQ=255)
                                                                                                    gcacgcacGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCcg                                                   <  1:3800629/150‑1 (MQ=255)
                                                                                                         cacgGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATg                                              <  1:1912602/150‑1 (MQ=255)
                                                                                                              tATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGa                                         >  2:2027175/1‑150 (MQ=255)
                                                                                                                tGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGAcc                                       <  1:4155593/150‑1 (MQ=255)
                                                                                                                   tCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGc                                    <  1:240842/150‑1 (MQ=255)
                                                                                                                     gATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCg                                  >  2:2725131/1‑150 (MQ=255)
                                                                                                                     gATTCCCTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCg                                  >  2:646675/1‑150 (MQ=255)
                                                                                                                      aTTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGa                                 >  1:3429093/1‑150 (MQ=255)
                                                                                                                       ttCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGAc                                >  1:3339986/1‑150 (MQ=255)
                                                                                                                       ttCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGAc                                >  1:253224/1‑150 (MQ=255)
                                                                                                                        tCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACt                               >  1:4065532/1‑150 (MQ=255)
                                                                                                                        tCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCACTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACt                               <  1:2025255/150‑1 (MQ=255)
                                                                                                                           cTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAtt                            <  2:1101943/150‑1 (MQ=255)
                                                                                                                           cTGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTAtt                            >  1:4511747/1‑150 (MQ=255)
                                                                                                                            tGCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTc                           >  2:4140121/1‑150 (MQ=255)
                                                                                                                             gCCGTACCGCAACGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCt                          <  1:1354082/150‑1 (MQ=255)
                                                                                                                                  aCCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGc                     >  2:3979399/1‑150 (MQ=255)
                                                                                                                                      cAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcggg                 <  2:3403828/150‑1 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc                >  1:1650747/1‑150 (MQ=255)
                                                                                                                                       aaCGACGCGGCGATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggc                <  2:2380799/150‑1 (MQ=255)
                                                                                                                                        acgacgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggcg               <  1:1017402/150‑1 (MQ=255)
                                                                                                                                          gacgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCcgggcggg             <  2:283770/150‑1 (MQ=255)
                                                                                                                                           acgCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCgggcgggg            <  2:1969328/150‑1 (MQ=255)
                                                                                                                                               ggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGAc        <  1:3133510/150‑1 (MQ=255)
                                                                                                                                               ggcgATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGAGAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGAc        >  2:3240948/1‑150 (MQ=255)
                                                                                                                                               ggcgATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGAc        >  1:2861761/1‑150 (MQ=255)
                                                                                                                                                 cgATAGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGc      >  2:2189663/1‑150 (MQ=255)
                                                                                                                                                     cGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTgcc  <  1:2609801/150‑1 (MQ=255)
                                                                                                                                                     cGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTgcc  <  1:3000137/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGGCATGCAGATGCAGGCGGCGGCGAAGGAGATCACCCGCAATGGCGGGATCCCGTTCGCGGCCTTCGTCAGCGATCCGTGCGACGGGCGCTCGCAGGGCACGCACGGTATGTTCGATTCCCTGCCGTACCGCAACGACGCGGCGATCGTGTTTCGCCGCCTGATCCGCTCCCTGCCGACGCGGCGGGCGGTGATCGGCGTAGCGACCTGCGATAAAGGGCTGCCCGCCACCATGATTGCGCTGGCCGCGATGCACGACCTGCCGACTATTCTGGTGCCGGGCGGGGCGACGCTGCC  >  CP009273/279185‑279483

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: