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,9610TG99.5% 761.8 / ‑5.3 185Q226P (CAA→CCA) dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (87/97);  total (88/97)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 192 total aligned reads displayed.

TGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTG  >  CP009273/1620816‑1621105
                                                                                                                                                 |                                                                                                                                                
tgtTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTggg                                                                                                                                                >  2:1842199/1‑148 (MQ=255)
 gtTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTgggg                                                                                                                                               <  2:1805357/148‑1 (MQ=255)
  ttGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGggt                                                                                                                                              >  1:3399838/1‑148 (MQ=255)
    gCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGggtgg                                                                                                                                            <  1:3050883/148‑1 (MQ=255)
        gttgtGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTAtt                                                                                                                                        >  2:404621/1‑148 (MQ=255)
         ttgtGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTc                                                                                                                                       <  1:2698570/148‑1 (MQ=255)
          tgtgCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTca                                                                                                                                      <  2:2235106/148‑1 (MQ=255)
              cACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATAATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTcagcag                                                                                                                                  <  1:2567948/148‑1 (MQ=255)
                cGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCg                                                                                                                                >  2:1123883/1‑148 (MQ=255)
                  gATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGa                                                                                                                              >  1:2307204/1‑148 (MQ=255)
                   aTCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGaa                                                                                                                             >  1:2835574/1‑148 (MQ=255)
                   aTCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGaa                                                                                                                             >  1:2022878/1‑148 (MQ=255)
                   aTCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGaa                                                                                                                             <  2:2760047/148‑1 (MQ=255)
                   aTCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTAACGCGTCGCACGATCGCGCATGTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGaa                                                                                                                             >  2:2402650/1‑148 (MQ=255)
                     ctcCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATc                                                                                                                           >  1:3236991/1‑148 (MQ=255)
                         accaGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGcc                                                                                                                       <  2:647690/148‑1 (MQ=255)
                           caGCCGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAt                                                                                                                     <  2:3111436/148‑1 (MQ=255)
                           caGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAt                                                                                                                     <  1:636805/148‑1 (MQ=255)
                             gACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAttt                                                                                                                   >  1:2616489/1‑148 (MQ=255)
                              aCGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTg                                                                                                                  >  2:1553161/1‑148 (MQ=255)
                              aCGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTg                                                                                                                  >  2:2093756/1‑148 (MQ=255)
                              aCGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTg                                                                                                                  >  2:1372564/1‑148 (MQ=255)
                               cGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGt                                                                                                                 <  1:3145906/148‑1 (MQ=255)
                                gTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGtt                                                                                                                >  2:1948215/1‑148 (MQ=255)
                                gTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGtt                                                                                                                <  2:2980143/148‑1 (MQ=255)
                                 ttGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTa                                                                                                               <  2:464142/148‑1 (MQ=255)
                                  tGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTAt                                                                                                              <  1:2702431/148‑1 (MQ=255)
                                  tGAATGATAGCGCGGGTAACATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTAt                                                                                                              <  2:1944402/148‑1 (MQ=255)
                                   gAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATc                                                                                                             <  2:3446950/148‑1 (MQ=255)
                                   gAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATc                                                                                                             >  1:1126461/1‑148 (MQ=255)
                                    aaTGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATcc                                                                                                            <  1:2149690/148‑1 (MQ=255)
                                      tGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAg                                                                                                          <  2:2412756/148‑1 (MQ=255)
                                      tGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAg                                                                                                          <  2:2415236/148‑1 (MQ=255)
                                       gATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGa                                                                                                         >  2:1697614/1‑148 (MQ=255)
                                       gATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGa                                                                                                         >  2:1696988/1‑148 (MQ=255)
                                       gATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGATTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGa                                                                                                         >  2:2129701/1‑148 (MQ=255)
                                         tAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGAcc                                                                                                       >  1:2744343/1‑148 (MQ=255)
                                            cgcgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCTCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCttt                                                                                                    <  2:471246/148‑1 (MQ=255)
                                            cgcgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCttt                                                                                                    <  2:3152012/148‑1 (MQ=255)
                                              cgGGTATCATTGGCCTCATTTTTGTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTGTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTtct                                                                                                  <  2:2357898/148‑1 (MQ=255)
                                              cgGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCAGTTGTTATCCAGACCTTtct                                                                                                  >  2:571920/1‑148 (MQ=255)
                                                ggTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTcg                                                                                                >  1:1627903/1‑148 (MQ=255)
                                                  tATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTcgcg                                                                                              >  1:1481443/1‑148 (MQ=255)
                                                  tATCATTGGCATCATTTTGTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTcgcg                                                                                              <  2:413356/148‑1 (MQ=255)
                                                          gCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCct                                                                                      >  1:556601/1‑148 (MQ=255)
                                                           catcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCctc                                                                                     >  2:186645/1‑148 (MQ=255)
                                                           catcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCctc                                                                                     >  2:4438/1‑148 (MQ=255)
                                                             tcatTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTg                                                                                   >  1:580957/1‑148 (MQ=255)
                                                                 ttttttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgccgc                                                                               >  2:2713657/1‑148 (MQ=255)
                                                                 ttttttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgccgc                                                                               <  1:1557385/148‑1 (MQ=255)
                                                                 ttttttCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTgccgc                                                                               <  1:2937015/148‑1 (MQ=255)
                                                                           tCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAtc                                                                     >  1:425551/1‑148 (MQ=255)
                                                                             gcgTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAAtctc                                                                   <  1:3365229/148‑1 (MQ=255)
                                                                                tCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCACCAGCGCAATCTCTTg                                                                >  2:311055/1‑148 (MQ=255)
                                                                                 ccAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGc                                                               >  2:3256259/1‑148 (MQ=255)
                                                                                     cccGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcac                                                           >  1:39773/1‑148 (MQ=255)
                                                                                       cGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcactc                                                         <  1:787877/148‑1 (MQ=255)
                                                                                        gCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcactca                                                        >  1:701661/1‑148 (MQ=255)
                                                                                        gCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcactca                                                        >  1:487965/1‑148 (MQ=255)
                                                                                        gCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGctcactca                                                        <  2:2835574/148‑1 (MQ=255)
                                                                                         cAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAt                                                       <  1:736015/148‑1 (MQ=255)
                                                                                         cAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAt                                                       >  2:2648546/1‑148 (MQ=255)
                                                                                         cAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAt                                                       >  2:1639286/1‑148 (MQ=255)
                                                                                           aTAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTc                                                     <  1:1927082/148‑1 (MQ=255)
                                                                                           aTAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTc                                                     <  1:1875286/148‑1 (MQ=255)
                                                                                           aTAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTc                                                     >  2:1101851/1‑148 (MQ=255)
                                                                                            tAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTcc                                                    >  2:1973383/1‑148 (MQ=255)
                                                                                            tAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTcc                                                    >  2:809291/1‑148 (MQ=255)
                                                                                            tAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTcc                                                    >  2:2515997/1‑148 (MQ=255)
                                                                                            tAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCAGTcc                                                    >  2:874520/1‑148 (MQ=255)
                                                                                             aaaCAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCt                                                   >  1:1650035/1‑148 (MQ=255)
                                                                                             aaaCAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCt                                                   >  2:2567897/1‑148 (MQ=255)
                                                                                              aaCAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTg                                                  <  1:1450052/148‑1 (MQ=255)
                                                                                               aCAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGc                                                 <  1:516169/148‑1 (MQ=255)
                                                                                                cAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGcc                                                >  2:3399509/1‑148 (MQ=255)
                                                                                                  gTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAg                                              <  1:170328/148‑1 (MQ=255)
                                                                                                  gTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAg                                              <  1:1009048/148‑1 (MQ=255)
                                                                                                   tttttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGc                                             <  2:3500370/148‑1 (MQ=255)
                                                                                                   tttttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGc                                             >  1:1651428/1‑148 (MQ=255)
                                                                                                   tttttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGc                                             >  1:1253752/1‑148 (MQ=255)
                                                                                                    ttttCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCt                                            >  1:1120872/1‑148 (MQ=255)
                                                                                                       tCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTgcg                                         <  2:504604/148‑1 (MQ=255)
                                                                                                            cgTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGata                                    <  1:950921/148‑1 (MQ=255)
                                                                                                             gTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGatat                                   <  1:2908213/148‑1 (MQ=255)
                                                                                                             gTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGatat                                   <  1:2524434/148‑1 (MQ=255)
                                                                                                                   cGATCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTca                             >  1:1165940/1‑148 (MQ=255)
                                                                                                                     aTCGCGCATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTcaca                           >  2:679773/1‑148 (MQ=255)
                                                                                                                           cATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGa                     <  1:1777702/148‑1 (MQ=255)
                                                                                                                           cATTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGa                     <  1:1329868/148‑1 (MQ=255)
                                                                                                                            aTTTCGGCAAGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGAc                    <  2:1310093/148‑1 (MQ=255)
                                                                                                                                    aaGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGat            <  2:3055133/148‑1 (MQ=255)
                                                                                                                                     aGCGCCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGatt           >  2:88736/1‑148 (MQ=255)
                                                                                                                                      gcgcCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGattt          <  2:792190/148‑1 (MQ=255)
                                                                                                                                      gcgcCGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGattt          >  1:2685987/1‑148 (MQ=255)
                                                                                                                                         ccGGTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGatttatt       <  2:1636774/148‑1 (MQ=255)
                                                                                                                                           ggTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtg     <  2:1196055/148‑1 (MQ=255)
                                                                                                                                           ggTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtg     <  2:1196025/148‑1 (MQ=255)
                                                                                                                                           ggTTGTGGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtg     >  1:941246/1‑148 (MQ=255)
                                                                                                                                            gttgtgGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtgc    <  1:1314426/148‑1 (MQ=255)
                                                                                                                                             ttgtggGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTtgct   >  1:2052535/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGTTGCCTGTTGTGCACGGATCTCCACCAGACGTTGAATGATAGCGCGGGTATCATTGGCATCATTTTTTTCCGCTCGCGTCCAGCCCGCAATAAACAGTTTTTCACGCGTCGCACGATCGCGCATTTCGGCAAGCGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTG  >  CP009273/1620816‑1621105

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: