Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,413 C→A 17.6% 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,4130CA17.6% 955.0 / 152.3 332I355I (ATC→ATAascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (150/123);  new base A (25/34);  total (175/157)
Fisher's exact test for biased strand distribution p-value = 8.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

GTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGA  >  CP009273/2835264‑2835562
                                                                                                                                                     |                                                                                                                                                     
gTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATc                                                                                                                                                       <  2:4405967/150‑1 (MQ=255)
gTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATc                                                                                                                                                       >  2:3699310/1‑150 (MQ=255)
        cAGCTATTACGCCTCGCGCGGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGaa                                                                                                                                               <  2:351755/150‑1 (MQ=255)
         aGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAat                                                                                                                                              >  1:3902219/1‑150 (MQ=255)
         aGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAat                                                                                                                                              >  1:4928462/1‑150 (MQ=255)
           cTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAatat                                                                                                                                            >  1:3602439/1‑150 (MQ=255)
           cTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAatat                                                                                                                                            >  2:4085452/1‑150 (MQ=255)
           cTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAatat                                                                                                                                            <  1:110112/150‑1 (MQ=255)
             aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatga                                                                                                                                          >  2:174917/1‑150 (MQ=255)
                aCGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGt                                                                                                                                       >  1:3033304/1‑150 (MQ=255)
                aCGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGt                                                                                                                                       >  1:2355167/1‑150 (MQ=255)
                  gCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGTGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTAc                                                                                                                                     <  1:2801066/150‑1 (MQ=255)
                    cTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGa                                                                                                                                   <  1:4357626/150‑1 (MQ=255)
                       gcgcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCg                                                                                                                                <  2:5003075/150‑1 (MQ=255)
                           tGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTAt                                                                                                                            >  1:1872462/1‑150 (MQ=255)
                           tGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTAt                                                                                                                            >  2:2937743/1‑150 (MQ=255)
                                 tCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAg                                                                                                                      >  2:1359602/1‑150 (MQ=255)
                                   ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGcc                                                                                                                    <  1:840481/150‑1 (MQ=255)
                                   ggcggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGcc                                                                                                                    >  1:1755716/1‑150 (MQ=255)
                                       gagaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTg                                                                                                                >  1:2151282/1‑150 (MQ=255)
                                        agaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGt                                                                                                               <  2:834883/150‑1 (MQ=255)
                                        agaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGt                                                                                                               >  2:5107023/1‑150 (MQ=255)
                                        agaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGt                                                                                                               <  1:2853266/150‑1 (MQ=255)
                                         gaTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGtt                                                                                                              <  1:1709027/150‑1 (MQ=255)
                                             aaCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtg                                                                                                          <  1:2982974/150‑1 (MQ=255)
                                             aaCGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtg                                                                                                          <  1:4587868/150‑1 (MQ=255)
                                               cGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                        <  1:3110018/150‑1 (MQ=255)
                                               cGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                        <  1:4338769/150‑1 (MQ=255)
                                               cGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                        >  1:2294152/1‑150 (MQ=255)
                                                        cagcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACgg                                                                                               >  2:3067378/1‑150 (MQ=255)
                                                         agcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGc                                                                                              >  1:1421936/1‑150 (MQ=255)
                                                         agcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGc                                                                                              >  2:292615/1‑150 (MQ=255)
                                                          gcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                                                             <  2:3254634/150‑1 (MQ=255)
                                                           cagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGTCTGTGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCt                                                                                            <  2:4332997/150‑1 (MQ=255)
                                                            agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTg                                                                                           <  1:3421958/150‑1 (MQ=255)
                                                            agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTg                                                                                           >  1:1055479/1‑150 (MQ=255)
                                                              tGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTggg                                                                                         >  1:781236/1‑150 (MQ=255)
                                                               gcggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGACTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                        >  2:2460319/1‑150 (MQ=255)
                                                               gcggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                        <  1:1666730/150‑1 (MQ=255)
                                                                 ggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGACTGGgcgc                                                                                      >  2:2548600/1‑150 (MQ=255)
                                                                  gcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCa                                                                                     <  2:566499/150‑1 (MQ=255)
                                                                   cgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGACGCTGTTTCTGGTGGAAAACGGCCTGGGCGCac                                                                                    >  1:3258876/1‑149 (MQ=255)
                                                                      aCGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAg                                                                                 <  2:3581655/150‑1 (MQ=255)
                                                                      aCGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAg                                                                                 <  2:3375298/150‑1 (MQ=255)
                                                                       cGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGa                                                                                >  1:4041781/1‑150 (MQ=255)
                                                                       cGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTCGGCGCAAAAGa                                                                                >  1:5104432/1‑150 (MQ=255)
                                                                       cGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGa                                                                                <  1:1121034/150‑1 (MQ=255)
                                                                        gtggtgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGAt                                                                               <  1:952752/150‑1 (MQ=255)
                                                                          ggtgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGa                                                                             >  2:1862936/1‑150 (MQ=255)
                                                                              aaaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAttt                                                                         >  1:555737/1‑150 (MQ=255)
                                                                               aaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtg                                                                        >  1:2654783/1‑150 (MQ=255)
                                                                               aaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGATATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtg                                                                        >  2:1425413/1‑150 (MQ=255)
                                                                                  cGCTGCGTTATCCGTATCTACTGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtgctg                                                                     <  2:2763240/150‑1 (MQ=255)
                                                                                  cGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtgctg                                                                     >  2:3149632/1‑150 (MQ=255)
                                                                                    cTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                                                   >  2:457091/1‑150 (MQ=255)
                                                                                     tGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCa                                                                  <  2:689431/150‑1 (MQ=255)
                                                                                       cGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                                >  2:4954794/1‑150 (MQ=255)
                                                                                         tAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATgg                                                              >  1:3641014/1‑150 (MQ=255)
                                                                                          aagccGTATCTACAGGTGAGCGACTGGGGGTGGGGAATTGCTCCACTCGGTCTGCGTATAACCATGAATATGAGGTACGATCGTGATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGACTTTGCTGCCAATGGc                                                             <  2:1934656/147‑1 (MQ=255)
                                                                                           aTCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCg                                                            >  1:3455891/1‑150 (MQ=255)
                                                                                            tCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCga                                                           <  2:1508100/150‑1 (MQ=255)
                                                                                              cGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgaga                                                         <  1:2240217/150‑1 (MQ=255)
                                                                                                tATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                                                       <  1:3632060/150‑1 (MQ=255)
                                                                                                 aTCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTa                                                      >  2:2164669/1‑150 (MQ=255)
                                                                                                 aTCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTa                                                      <  2:934818/150‑1 (MQ=255)
                                                                                                 aTCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTa                                                      >  2:3537694/1‑150 (MQ=255)
                                                                                                  tCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTaa                                                     >  2:1128393/1‑150 (MQ=255)
                                                                                                   cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGTAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                    <  2:2585522/150‑1 (MQ=255)
                                                                                                   cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                    <  1:3704048/150‑1 (MQ=255)
                                                                                                   cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                    <  1:2486167/150‑1 (MQ=255)
                                                                                                       aGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                                >  1:1531538/1‑150 (MQ=255)
                                                                                                       aGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                                >  1:1250802/1‑150 (MQ=255)
                                                                                                         gTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACgacgac                                              <  1:2199560/150‑1 (MQ=255)
                                                                                                          tGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt                                             >  2:2370620/1‑150 (MQ=255)
                                                                                                            aGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                           <  1:4168279/150‑1 (MQ=255)
                                                                                                            aGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                           <  2:4162982/150‑1 (MQ=255)
                                                                                                            aGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                           >  2:3935402/1‑150 (MQ=255)
                                                                                                               gACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGc                                        >  1:3720678/1‑150 (MQ=255)
                                                                                                                 ctggggctggggAATTGATCCACTCGGTCTTCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGATTTTGCTGCCAATGGCGAGATTAACGACGACTATCGCAt                                      <  1:2039854/150‑1 (MQ=255)
                                                                                                                 ctggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCAt                                      <  1:2549398/150‑1 (MQ=255)
                                                                                                                  tggggctggggAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                     <  2:841649/150‑1 (MQ=255)
                                                                                                                   ggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCa                                    >  1:4700849/1‑150 (MQ=255)
                                                                                                                        ttgtgatttGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGAAAAAGAAGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTAc                               <  2:4779067/143‑1 (MQ=255)
                                                                                                                         ggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACt                              >  2:4239209/1‑150 (MQ=255)
                                                                                                                             aaTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTAcg                          >  1:597826/1‑150 (MQ=255)
                                                                                                                             aaTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTAcg                          >  1:4424846/1‑150 (MQ=255)
                                                                                                                                tGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGa                       <  2:3853160/150‑1 (MQ=255)
                                                                                                                                  aTCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAAc                     >  1:4081306/1‑150 (MQ=255)
                                                                                                                                  aTCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAAc                     >  1:3499771/1‑150 (MQ=255)
                                                                                                                                    ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACat                   >  1:4729686/1‑150 (MQ=255)
                                                                                                                                    ccACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACat                   <  1:862050/150‑1 (MQ=255)
                                                                                                                                    ccACTCGGTCTGCGTATCAACATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACat                   <  2:3370376/150‑1 (MQ=255)
                                                                                                                                      aCTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat                 <  1:4868930/150‑1 (MQ=255)
                                                                                                                                      aCTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat                 >  1:1036502/1‑150 (MQ=255)
                                                                                                                                        tCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATcc               <  1:1331206/150‑1 (MQ=255)
                                                                                                                                        tCGGTCTGCGTATCACCATGAATATGATCTACGATAGTTATCAGAAGCCGATGTTTCTGGTGGAAAAAGTCCTGGGCGCAAAATATGAATATGCTGCAAATGGCGAGATTAACTACGACTATCGAATCAGATACTTACGAGAACATATca               >  2:407219/1‑149 (MQ=255)
                                                                                                                                          ggTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgc             <  1:1986776/150‑1 (MQ=255)
                                                                                                                                                cGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg       >  1:4253288/1‑150 (MQ=255)
                                                                                                                                                  tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg     >  1:4491477/1‑150 (MQ=255)
                                                                                                                                                     cACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGa  >  1:2453890/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGA  >  CP009273/2835264‑2835562

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: