Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,620,961 T→G 100% Q226P (CAA→CCA)  dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,620,9610TG97.0% 906.5 / 6.7 229Q226P (CAA→CCA) dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (99/123);  total (103/126)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 234 total aligned reads displayed.

AGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCC  >  CP009273/1620814‑1621107
                                                                                                                                                   |                                                                                                                                                  
      gCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTa                                                                                                                                            <  1:1064055/150‑1 (MQ=255)
       ccTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTAt                                                                                                                                           >  1:2410874/1‑150 (MQ=255)
        cTGGTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTAtt                                                                                                                                          <  2:3314612/150‑1 (MQ=255)
         tgttgtGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTc                                                                                                                                         <  1:124182/150‑1 (MQ=255)
          gttgtGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTca                                                                                                                                        <  1:786570/150‑1 (MQ=255)
           ttgtGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTcag                                                                                                                                       >  2:2168165/1‑150 (MQ=255)
            tgtgCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTcagc                                                                                                                                      >  1:740086/1‑150 (MQ=255)
             gtgCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTcagca                                                                                                                                     <  1:4114426/150‑1 (MQ=255)
                  cGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGa                                                                                                                                <  1:1993459/150‑1 (MQ=255)
                   ggATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGaa                                                                                                                               <  2:2870337/150‑1 (MQ=255)
                    gATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAAt                                                                                                                              <  2:527476/150‑1 (MQ=255)
                     aTCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATc                                                                                                                             <  1:493911/150‑1 (MQ=255)
                       ctcCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAg                                                                                                                           <  2:1710950/150‑1 (MQ=255)
                         ccaccaGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGcc                                                                                                                         <  2:1248877/150‑1 (MQ=255)
                           accaGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATACACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAt                                                                                                                       <  2:2995685/150‑1 (MQ=255)
                            ccaGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCACGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAtt                                                                                                                      >  2:3664387/1‑150 (MQ=255)
                             caGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAttt                                                                                                                     <  1:2367769/150‑1 (MQ=255)
                             caGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAttt                                                                                                                     <  2:4397971/150‑1 (MQ=255)
                             caGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAgtt                                                                                                                     >  2:1476561/1‑150 (MQ=255)
                              aGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTg                                                                                                                    <  1:1886764/150‑1 (MQ=255)
                               gACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGt                                                                                                                   <  2:615044/150‑1 (MQ=255)
                                 cGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTa                                                                                                                 <  1:221352/150‑1 (MQ=255)
                                  gTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTAt                                                                                                                >  1:3522223/1‑150 (MQ=255)
                                     gAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCa                                                                                                             <  1:1120525/150‑1 (MQ=255)
                                        tGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGAc                                                                                                          >  1:2621932/1‑150 (MQ=255)
                                         gATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGAcc                                                                                                         >  1:4183624/1‑150 (MQ=255)
                                          aTAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCt                                                                                                        <  2:3638775/150‑1 (MQ=255)
                                           tAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCtt                                                                                                       >  2:1885899/1‑150 (MQ=255)
                                           tAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCtt                                                                                                       >  1:2788709/1‑150 (MQ=255)
                                           tAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCtt                                                                                                       >  2:1372979/1‑150 (MQ=255)
                                               gcgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGNAATCAGCCATTTGTTATCCAGACCTTtctc                                                                                                   <  1:4281212/150‑1 (MQ=255)
                                               gcgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTtctc                                                                                                   <  1:2766796/150‑1 (MQ=255)
                                                cgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTcg                                                                                                  >  1:889897/1‑150 (MQ=255)
                                                  ggTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTcgcg                                                                                                >  2:4088864/1‑150 (MQ=255)
                                                    tATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAg                                                                                              <  2:4419465/150‑1 (MQ=255)
                                                      tCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAgcc                                                                                            >  2:1594982/1‑150 (MQ=255)
                                                      tCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAgcc                                                                                            <  1:4289252/150‑1 (MQ=255)
                                                      tCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAgcc                                                                                            <  1:3790081/150‑1 (MQ=255)
                                                           ggCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCctc                                                                                       <  1:4535101/150‑1 (MQ=255)
                                                           ggCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCctc                                                                                       >  1:275816/1‑150 (MQ=255)
                                                             catcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTg                                                                                     <  2:2098061/150‑1 (MQ=255)
                                                             catcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTg                                                                                     >  2:2897091/1‑150 (MQ=255)
                                                               tcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgcc                                                                                   <  2:4462946/150‑1 (MQ=255)
                                                               tcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgcc                                                                                   <  1:636020/150‑1 (MQ=255)
                                                               tcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgcc                                                                                   <  1:4262214/150‑1 (MQ=255)
                                                                 atTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgccgc                                                                                 >  2:3499187/1‑150 (MQ=255)
                                                                    tttttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAg                                                                              >  2:1441435/1‑150 (MQ=255)
                                                                    tttttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAg                                                                              <  1:2944421/150‑1 (MQ=255)
                                                                       ttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAgcgc                                                                           >  2:4399866/1‑150 (MQ=255)
                                                                       ttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAgcgc                                                                           >  1:1725234/1‑150 (MQ=255)
                                                                       ttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAgcgc                                                                           >  1:3001428/1‑150 (MQ=255)
                                                                          cGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAt                                                                        <  2:2383461/150‑1 (MQ=255)
                                                                           gCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCATACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAtc                                                                       >  1:1194321/1‑150 (MQ=255)
                                                                             tCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAtctc                                                                     >  1:3612552/1‑150 (MQ=255)
                                                                              cgcgTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATctct                                                                    >  1:2654977/1‑150 (MQ=255)
                                                                               gcgTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAtctcgt                                                                   <  2:1609764/150‑3 (MQ=255)
                                                                               gcgTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCtt                                                                   <  2:174272/150‑1 (MQ=255)
                                                                                 gTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGc                                                                 >  1:4522315/1‑150 (MQ=255)
                                                                                  tCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGct                                                                >  2:3541564/1‑150 (MQ=255)
                                                                                  tCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGct                                                                <  1:3371268/150‑1 (MQ=255)
                                                                                    cAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctca                                                              >  2:1069634/1‑150 (MQ=255)
                                                                                     aGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcac                                                             <  1:956547/150‑1 (MQ=255)
                                                                                       cccGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcactc                                                           <  2:3801492/150‑1 (MQ=255)
                                                                                         cGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAt                                                         <  1:1490506/150‑1 (MQ=255)
                                                                                         cGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAt                                                         >  2:3066310/1‑150 (MQ=255)
                                                                                              tAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTg                                                    <  1:3687837/150‑1 (MQ=255)
                                                                                                  cAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAg                                                <  2:3944972/150‑1 (MQ=255)
                                                                                                        ttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTgcgc                                          <  1:2125924/150‑1 (MQ=255)
                                                                                                        ttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTgcgc                                          >  2:3746883/1‑150 (MQ=255)
                                                                                                         tCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGcgcg                                         <  2:4328580/150‑1 (MQ=255)
                                                                                                         tCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGcgcg                                         >  2:3280372/1‑150 (MQ=255)
                                                                                                            cgcgTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGata                                      >  2:3262546/1‑150 (MQ=255)
                                                                                                              cgTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATc                                    >  1:3291879/1‑150 (MQ=255)
                                                                                                                 cGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGtt                                 >  2:3080946/1‑150 (MQ=255)
                                                                                                                  gCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTc                                <  1:2327761/150‑1 (MQ=255)
                                                                                                                   cACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTca                               >  1:921560/1‑150 (MQ=255)
                                                                                                                   cACGATCGCGCATTTCGGCAAGCGCCGCTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATATGTTATCCAGACCTTACTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTAGCTCACTCAGTCCTGCCAGCTGCGCGATATCGATca                               >  2:3329674/1‑150 (MQ=255)
                                                                                                                    aCGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTcac                              <  1:4394500/150‑1 (MQ=255)
                                                                                                                    aCGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTcac                              >  2:424411/1‑150 (MQ=255)
                                                                                                                    aCGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTcac                              >  2:424065/1‑150 (MQ=255)
                                                                                                                      gATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACaa                            >  2:2108891/1‑150 (MQ=255)
                                                                                                                      gATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACaa                            >  2:822036/1‑150 (MQ=255)
                                                                                                                      gATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACaa                            >  2:4490019/1‑150 (MQ=255)
                                                                                                                       aTCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAAc                           <  2:3176313/150‑1 (MQ=255)
                                                                                                                         cgcgCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCa                         <  1:2464458/150‑1 (MQ=255)
                                                                                                                            gcATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGAc                      >  2:4264358/1‑150 (MQ=255)
                                                                                                                              aTTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGAccg                    >  2:4569905/1‑150 (MQ=255)
                                                                                                                              aTTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGAccg                    <  2:566832/150‑1 (MQ=255)
                                                                                                                                ttCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGAccgcc                  <  2:113438/150‑1 (MQ=255)
                                                                                                                                   ggCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGa               <  2:3916176/150‑1 (MQ=255)
                                                                                                                                   ggCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGa               <  1:2848312/150‑1 (MQ=255)
                                                                                                                                    gCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGat              <  2:2731688/150‑1 (MQ=255)
                                                                                                                                     cAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGatt             >  1:1326043/1‑150 (MQ=255)
                                                                                                                                       aGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGattta           >  2:735812/1‑150 (MQ=255)
                                                                                                                                         cgcCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGatttatt         >  2:4622202/1‑150 (MQ=255)
                                                                                                                                          gcCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGatttattt        <  2:3604904/150‑1 (MQ=255)
                                                                                                                                           ccGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtg       >  2:2732259/1‑150 (MQ=255)
                                                                                                                                           ccGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtg       <  1:2091908/150‑1 (MQ=255)
                                                                                                                                            cGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtgc      <  2:894215/150‑1 (MQ=255)
                                                                                                                                                tgtggGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCATACCGCCGGATTTATTTGCTGcc  >  2:2964837/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCC  >  CP009273/1620814‑1621107

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: