Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 207,709 C→A 5.5% A515E (GCA→GAA)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273207,7090CA 5.5% 987.1 / 22.6 256A515E (GCA→GAA) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (122/120);  new base A (10/4);  total (132/124)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 258 total aligned reads displayed.

CCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGC  >  CP009273/207561‑207858
                                                                                                                                                    |                                                                                                                                                     
cccGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCa                                                                                                                                                      >  1:4423102/1‑150 (MQ=255)
cccGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCa                                                                                                                                                      <  1:2333666/150‑1 (MQ=255)
   gTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCaaaa                                                                                                                                                   >  1:80816/1‑150 (MQ=255)
    ttGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAAt                                                                                                                                                  <  1:1245589/150‑1 (MQ=255)
    ttGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAAt                                                                                                                                                  >  1:1492559/1‑150 (MQ=255)
    ttGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAAt                                                                                                                                                  <  1:3058088/150‑1 (MQ=255)
           tGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTcga                                                                                                                                           >  1:4551865/1‑150 (MQ=255)
                 aCAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACg                                                                                                                                     <  1:699534/150‑1 (MQ=255)
                    gTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTgg                                                                                                                                  <  2:4109772/150‑1 (MQ=255)
                       gCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGAt                                                                                                                               <  1:139656/150‑1 (MQ=255)
                        cACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATc                                                                                                                              <  2:3289442/150‑1 (MQ=255)
                         aCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                                                                                             <  1:3518252/150‑1 (MQ=255)
                         aCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                                                                                             <  1:2436228/150‑1 (MQ=255)
                         aCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                                                                                             <  1:1519883/150‑1 (MQ=255)
                         aCGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                                                                                             <  1:1519830/150‑1 (MQ=255)
                            gCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag                                                                                                                          >  1:384242/1‑150 (MQ=255)
                            gCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag                                                                                                                          >  1:2818119/1‑150 (MQ=255)
                                 aaCGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGAAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                                                                                                     <  1:1163461/150‑1 (MQ=255)
                                    gATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                                                                                                  <  1:3855774/150‑1 (MQ=255)
                                       gCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                                                                               <  2:3838036/150‑1 (MQ=255)
                                        cGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCg                                                                                                              >  1:4401773/1‑150 (MQ=255)
                                             gCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCt                                                                                                         >  1:2171744/1‑150 (MQ=255)
                                              ccGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                                                                                                        <  2:153968/150‑1 (MQ=255)
                                              ccGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                                                                                                        >  1:474420/1‑150 (MQ=255)
                                              ccGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCtt                                                                                                        >  2:219124/1‑150 (MQ=255)
                                               cGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                                                                       >  2:2064967/1‑150 (MQ=255)
                                               cGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                                                                                       <  2:42181/150‑1 (MQ=255)
                                                   cATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                                                                                   >  2:4574741/1‑150 (MQ=255)
                                                         tttCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgctg                                                                                             <  2:3106212/150‑1 (MQ=255)
                                                          ttCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGt                                                                                            >  1:2525968/1‑150 (MQ=255)
                                                          ttCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGt                                                                                            <  1:1734976/150‑1 (MQ=255)
                                                           tctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGtt                                                                                           <  2:3994609/150‑1 (MQ=255)
                                                           tctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGtt                                                                                           >  2:367233/1‑150 (MQ=255)
                                                           tctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGtt                                                                                           >  1:986745/1‑150 (MQ=255)
                                                           tctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGtt                                                                                           <  2:156788/150‑1 (MQ=255)
                                                            ctcGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGttt                                                                                          >  1:173323/1‑150 (MQ=255)
                                                              cGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtct                                                                                        >  2:630842/1‑150 (MQ=255)
                                                               gATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                                       >  1:1000981/1‑150 (MQ=255)
                                                               gATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGAAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                                                                                       >  1:3428168/1‑150 (MQ=255)
                                                                aTCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTctct                                                                                      >  1:101524/1‑150 (MQ=255)
                                                                  ccGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                                    >  1:754071/1‑150 (MQ=255)
                                                                  ccGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCttt                                                                                    <  2:4388729/150‑1 (MQ=255)
                                                                   cGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa                                                                                   <  2:4087806/150‑1 (MQ=255)
                                                                   cGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa                                                                                   <  2:2808333/150‑1 (MQ=255)
                                                                    ggTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                                                                                  <  1:431278/150‑1 (MQ=255)
                                                                    ggTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                                                                                  <  2:1226/150‑1 (MQ=255)
                                                                    ggTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                                                                                  <  2:2881917/150‑1 (MQ=255)
                                                                       tAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                                                                               >  2:4443436/1‑150 (MQ=255)
                                                                       tAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                                                                               <  1:2127000/150‑1 (MQ=255)
                                                                       tAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                                                                               >  2:1454088/1‑150 (MQ=255)
                                                                          aGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTgg                                                                            <  2:587191/150‑1 (MQ=255)
                                                                           gTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGc                                                                           <  1:2693040/150‑1 (MQ=255)
                                                                            tCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa                                                                          <  2:2371800/150‑1 (MQ=255)
                                                                             cACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt                                                                         >  2:1169179/1‑150 (MQ=255)
                                                                              aCTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc                                                                        >  1:2857734/1‑150 (MQ=255)
                                                                               cTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCg                                                                       <  1:3019461/150‑1 (MQ=255)
                                                                                tATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGa                                                                      >  1:1830364/1‑150 (MQ=255)
                                                                                tATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGa                                                                      >  2:3417859/1‑150 (MQ=255)
                                                                                  ttttGACACCGGGGATGGACGAGCAGGGCAGTATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATa                                                                    >  1:2378752/1‑150 (MQ=255)
                                                                                       acacCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAcc                                                               >  2:3594771/1‑150 (MQ=255)
                                                                                       acacCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAcc                                                               >  1:1961097/1‑150 (MQ=255)
                                                                                        cacCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                              >  1:4592040/1‑150 (MQ=255)
                                                                                        cacCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa                                                              >  1:1647955/1‑150 (MQ=255)
                                                                                         acCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaa                                                             >  1:2047227/1‑150 (MQ=255)
                                                                                         acCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCaa                                                             <  2:2599304/150‑1 (MQ=255)
                                                                                           cGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAg                                                           <  2:2390400/150‑1 (MQ=255)
                                                                                            ggggATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGc                                                          >  2:843363/1‑150 (MQ=255)
                                                                                             gggATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa                                                         >  1:130236/1‑150 (MQ=255)
                                                                                                  ggACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATggg                                                    >  2:1715488/1‑150 (MQ=255)
                                                                                                   gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGAAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGa                                                   >  1:224656/1‑150 (MQ=255)
                                                                                                    aCGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat                                                  <  1:547130/150‑1 (MQ=255)
                                                                                                       aGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattat                                               >  1:4151096/1‑150 (MQ=255)
                                                                                                       aGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattat                                               >  2:2993445/1‑150 (MQ=255)
                                                                                                          aGGGCAATATGAGCGAGGAGGGGATCCCGGAGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGTGATTATTGc                                            >  1:2379514/1‑150 (MQ=255)
                                                                                                               aaTATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                       >  1:2092988/1‑150 (MQ=255)
                                                                                                               aaTATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                       <  2:1820757/150‑1 (MQ=255)
                                                                                                               aaTATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTAGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTggg                                       >  2:1373410/1‑150 (MQ=255)
                                                                                                                 tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGtt                                     <  2:17536/150‑1 (MQ=255)
                                                                                                                 tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGtt                                     <  2:2517531/150‑1 (MQ=255)
                                                                                                                  atGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTg                                    <  2:2838704/150‑1 (MQ=255)
                                                                                                                     aGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACg                                 >  2:4626323/1‑150 (MQ=255)
                                                                                                                     aGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGGTGACg                                 >  2:770672/1‑150 (MQ=255)
                                                                                                                       cGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGa                               >  2:3113319/1‑150 (MQ=255)
                                                                                                                       cGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGa                               >  2:372659/1‑150 (MQ=255)
                                                                                                                        gaggagGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGaa                              <  1:589094/150‑1 (MQ=255)
                                                                                                                          ggaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAAtt                            >  1:2632092/1‑150 (MQ=255)
                                                                                                                          ggaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAAtt                            <  1:3707211/150‑1 (MQ=255)
                                                                                                                           gaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTc                           <  2:277816/150‑1 (MQ=255)
                                                                                                                           gaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTc                           <  2:1326798/150‑1 (MQ=255)
                                                                                                                           gaggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTc                           <  1:3533711/150‑1 (MQ=255)
                                                                                                                            aggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCa                          >  1:2720671/1‑150 (MQ=255)
                                                                                                                              gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCaaa                        >  2:2578312/1‑150 (MQ=255)
                                                                                                                               ggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAAc                       <  2:2000308/150‑1 (MQ=255)
                                                                                                                                     cggcggCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCtt                 <  1:4469641/150‑1 (MQ=255)
                                                                                                                                         ggcgCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGa             >  1:1608472/1‑150 (MQ=255)
                                                                                                                                          gcgcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAt            <  2:3910515/150‑1 (MQ=255)
                                                                                                                                             cTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGAtctc         >  1:4224511/1‑150 (MQ=255)
                                                                                                                                                  aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCt    <  2:80816/150‑1 (MQ=255)
                                                                                                                                                  aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCt    <  1:4266111/150‑1 (MQ=255)
                                                                                                                                                    cAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGc  <  2:297924/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCGTTGCGCCTGGCGAACAGTGGCACGGCTTTAACGATGCGGATGCCGATCATATGTTTCTCGATCCGGTTAAAGTCACTATTTTGACACCGGGGATGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTGGGTTGACGGAATTCAAACGCTCTTACGATCTCAACCTGC  >  CP009273/207561‑207858

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: