Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,413 C→A 13.4% I355I (ATC→ATA ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,835,4130CA13.4% 918.6 / 96.2 276I355I (ATC→ATAascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (101/138);  new base A (12/25);  total (113/163)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 278 total aligned reads displayed.

TCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGC  >  CP009273/2835265‑2835560
                                                                                                                                                    |                                                                                                                                                   
  tctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAc                                                                                                                                                   >  2:848290/1‑149 (MQ=255)
   ctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAcc                                                                                                                                                  <  1:1090598/149‑1 (MQ=255)
     ttCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAt                                                                                                                                                >  1:2520971/1‑149 (MQ=255)
      tCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATg                                                                                                                                               >  1:3737186/1‑149 (MQ=255)
      tCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATg                                                                                                                                               >  1:4061528/1‑149 (MQ=255)
        aGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGaa                                                                                                                                             >  2:380458/1‑149 (MQ=255)
           tATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAatat                                                                                                                                          <  2:1021975/149‑1 (MQ=255)
           tATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAatat                                                                                                                                          >  1:70542/1‑149 (MQ=255)
            aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatg                                                                                                                                         >  2:775899/1‑149 (MQ=255)
            aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatg                                                                                                                                         >  2:1759251/1‑149 (MQ=255)
              tACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatgat                                                                                                                                       <  2:2496842/149‑1 (MQ=255)
              tACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATatgat                                                                                                                                       <  1:528192/149‑1 (MQ=255)
                 gCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTa                                                                                                                                    <  1:3699910/149‑1 (MQ=255)
                   cTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACg                                                                                                                                  <  2:2283071/149‑1 (MQ=255)
                     cgcgCTGCGCCTCGGCGGAGATGAACGCCAACAACCGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGAt                                                                                                                                <  1:2093678/149‑1 (MQ=255)
                         cTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGtt                                                                                                                            <  1:3245153/149‑1 (MQ=255)
                            cgcCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATc                                                                                                                         <  2:4301470/149‑1 (MQ=255)
                              ccTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAg                                                                                                                       >  1:453263/1‑149 (MQ=255)
                               cTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACCGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGa                                                                                                                      >  2:1454866/1‑149 (MQ=255)
                                tCGGCGGAGATGAACGCCTACAACAGCAGTGCGGCGAACGTGGTGTAATCGCTGCGTAATCCGTATCTACAGGTGTGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGGATCACCATGAATATGATGTACGATCGTTATCAGaa                                                                                                                     <  2:664248/149‑1 (MQ=255)
                                     ggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGc                                                                                                                <  2:3394764/149‑1 (MQ=255)
                                      gagaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCt                                                                                                               >  1:3459398/1‑149 (MQ=255)
                                       agaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTg                                                                                                              >  2:2385683/1‑149 (MQ=255)
                                       agaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTg                                                                                                              <  2:241772/149‑1 (MQ=255)
                                          tGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGttt                                                                                                           <  1:1399559/149‑1 (MQ=255)
                                           gAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTc                                                                                                          <  2:4473677/149‑1 (MQ=255)
                                            aaCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCt                                                                                                         <  2:1105376/149‑1 (MQ=255)
                                               gCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                      >  2:622106/1‑149 (MQ=255)
                                                ccAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggtg                                                                                                     <  2:4655463/149‑1 (MQ=255)
                                                ccAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggtg                                                                                                     >  2:954314/1‑149 (MQ=255)
                                                  aacaacAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGa                                                                                                   <  1:1142579/149‑1 (MQ=255)
                                                     aacaGCAGTGCGGCGAACGTGGTTAAATCGCTGCTTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaaaa                                                                                                <  1:1570060/149‑1 (MQ=255)
                                                      acaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                               <  2:531466/149‑1 (MQ=255)
                                                          cagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                                                           <  2:1650073/149‑1 (MQ=255)
                                                           agTGCGGCGAACGTGGTTAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCt                                                                                          <  1:1368244/149‑1 (MQ=255)
                                                           agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAATCGGCCt                                                                                          >  2:2433425/1‑149 (MQ=255)
                                                              gcggcgAACGTGGTGAAATCGCTGCTTAATCCGTATCTACAGGGGCGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTggg                                                                                       <  1:3578726/149‑1 (MQ=255)
                                                               cggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                      <  1:2251115/149‑1 (MQ=255)
                                                                 gcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgcgc                                                                                    <  2:4231177/149‑1 (MQ=255)
                                                                     aCGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCaaaa                                                                                <  1:1905606/149‑1 (MQ=255)
                                                                           ttaaaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                                                          <  1:3000397/147‑1 (MQ=255)
                                                                                tCGCTGCGTAATCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtgc                                                                     <  1:2004668/149‑1 (MQ=255)
                                                                                  gCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtgctg                                                                   >  1:1542041/1‑149 (MQ=255)
                                                                                    tGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                                                 >  1:4657850/1‑149 (MQ=255)
                                                                                     gCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCa                                                                >  2:1860214/1‑149 (MQ=255)
                                                                                       gTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                              >  1:1316351/1‑149 (MQ=255)
                                                                                       gTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTACGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                              <  1:1536263/149‑1 (MQ=255)
                                                                                         aaTCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATgg                                                            <  1:2830311/149‑1 (MQ=255)
                                                                                         aaTCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATgg                                                            <  1:1716868/149‑1 (MQ=255)
                                                                                         aaTCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATCATGTACGATCGTTATCAGAAGCCGCTGTTTCGGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATgg                                                            >  1:2924688/1‑149 (MQ=255)
                                                                                           tCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCg                                                          <  2:1290716/149‑1 (MQ=255)
                                                                                             cGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgag                                                        <  1:4437569/149‑1 (MQ=255)
                                                                                               tATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTGATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAt                                                      <  2:2858153/149‑1 (MQ=255)
                                                                                                aTCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                                                     >  1:3993990/1‑149 (MQ=255)
                                                                                                  cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTaa                                                   <  1:2305513/149‑1 (MQ=255)
                                                                                                   tACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                  >  1:4383457/1‑149 (MQ=255)
                                                                                                   tACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                  <  1:4457029/149‑1 (MQ=255)
                                                                                                     cAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacga                                                <  1:347043/149‑1 (MQ=255)
                                                                                                     cAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacga                                                >  2:1205147/1‑149 (MQ=255)
                                                                                                      tggTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac                                               <  1:2959850/148‑1 (MQ=255)
                                                                                                        gtgtGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAAcgacga                                             <  1:1934405/149‑1 (MQ=255)
                                                                                                        gTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAAcgacga                                             >  2:3188506/1‑149 (MQ=255)
                                                                                                        gTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAAcgacga                                             <  2:2719809/149‑1 (MQ=255)
                                                                                                        gTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAAcgacga                                             <  1:2730369/149‑1 (MQ=255)
                                                                                                         tGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACgacgac                                            <  2:2089596/149‑1 (MQ=255)
                                                                                                            gCGACTGGGGGTGGGGAATTGATCCTCTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                         <  1:1740850/149‑1 (MQ=255)
                                                                                                            gCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                         <  1:4448847/149‑1 (MQ=255)
                                                                                                               aCTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGc                                      >  1:3857692/1‑149 (MQ=255)
                                                                                                               aCTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGc                                      >  2:3280201/1‑149 (MQ=255)
                                                                                                                 tggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCAt                                    <  1:3870120/149‑1 (MQ=255)
                                                                                                                   gggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCa                                  <  2:109144/149‑1 (MQ=255)
                                                                                                                       tggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTa                              >  2:26010/1‑149 (MQ=255)
                                                                                                                        ggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTAc                             >  1:3400487/1‑149 (MQ=255)
                                                                                                                         gggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACt                            >  2:2084567/1‑149 (MQ=255)
                                                                                                                         gggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACt                            <  2:3768/149‑1 (MQ=255)
                                                                                                                          ggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACtt                           <  2:3709268/149‑1 (MQ=255)
                                                                                                                          ggAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACtt                           <  2:2019310/149‑1 (MQ=255)
                                                                                                                           gAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTa                          <  2:3617720/149‑1 (MQ=255)
                                                                                                                             aTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAGGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTAcg                        >  1:10469/1‑149 (MQ=255)
                                                                                                                                gATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGTCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGa                     >  1:202015/1‑149 (MQ=255)
                                                                                                                                   ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACa                  >  2:3897517/1‑149 (MQ=255)
                                                                                                                                   ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACa                  >  2:4154615/1‑149 (MQ=255)
                                                                                                                                   ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACa                  <  2:2686229/149‑1 (MQ=255)
                                                                                                                                   ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACa                  >  1:2142665/1‑149 (MQ=255)
                                                                                                                                   ccACTCGGTCTGCGTATCACCATGAATAGGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACa                  <  2:3576232/149‑1 (MQ=255)
                                                                                                                                     aCTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACata                >  2:3146925/1‑149 (MQ=255)
                                                                                                                                      cTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat               <  1:3074740/149‑1 (MQ=255)
                                                                                                                                      cTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat               >  2:4052170/1‑149 (MQ=255)
                                                                                                                                      cTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat               >  2:2168434/1‑149 (MQ=255)
                                                                                                                                        cGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATcc             <  1:2056932/149‑1 (MQ=255)
                                                                                                                                         ggTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCACTGGCGAGATTAACGACGACTCTCGCATCAGCTACTTACGCGAACATATCcg            >  2:2919518/1‑149 (MQ=255)
                                                                                                                                         ggTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcg            >  2:1828882/1‑149 (MQ=255)
                                                                                                                                          gTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgc           <  1:2609771/149‑1 (MQ=255)
                                                                                                                                              gCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCaa       >  2:4297231/1‑149 (MQ=255)
                                                                                                                                              gCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCaa       <  2:224635/149‑1 (MQ=255)
                                                                                                                                                gTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg     >  1:2383473/1‑149 (MQ=255)
                                                                                                                                                gTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg     <  2:3235294/149‑1 (MQ=255)
                                                                                                                                                gTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg     <  2:1490631/149‑1 (MQ=255)
                                                                                                                                                   tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGc  >  2:1545158/1‑149 (MQ=255)
                                                                                                                                                   tCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGc  >  1:828395/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGC  >  CP009273/2835265‑2835560

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: