Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,835,413 C→A 11.9% 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,4130CA11.9% 1188.7 / 101.7 352I355I (ATC→ATAascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (156/154);  new base A (26/16);  total (182/170)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 358 total aligned reads displayed.

GTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGA  >  CP009273/2835264‑2835562
                                                                                                                                                     |                                                                                                                                                     
gTCTCTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATc                                                                                                                                                       >  2:760514/1‑150 (MQ=255)
 tctcTTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCa                                                                                                                                                      <  1:3268281/150‑1 (MQ=255)
  ctctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAc                                                                                                                                                     <  1:3345106/150‑1 (MQ=255)
  ctctTTCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCAc                                                                                                                                                     >  2:1673878/1‑150 (MQ=255)
     tttCAGCTATTACGCCTCTCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAt                                                                                                                                                  <  2:3535659/150‑1 (MQ=255)
      ttCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATg                                                                                                                                                 >  1:3045852/1‑150 (MQ=255)
       tCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGa                                                                                                                                                <  2:1823296/150‑1 (MQ=255)
       tCAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGa                                                                                                                                                >  2:4257259/1‑150 (MQ=255)
        cAGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGaa                                                                                                                                               >  2:3384160/1‑150 (MQ=255)
         aGCTATTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAat                                                                                                                                              <  2:3201731/150‑1 (MQ=255)
             aTTACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatga                                                                                                                                          >  2:1377447/1‑150 (MQ=255)
               tACGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATatgatg                                                                                                                                        >  2:2572253/1‑150 (MQ=255)
                aCGCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGt                                                                                                                                       >  1:581340/1‑150 (MQ=255)
                  gCCTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTAc                                                                                                                                     >  1:1442550/1‑150 (MQ=255)
                   ccTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACg                                                                                                                                    >  2:2036209/1‑150 (MQ=255)
                   ccTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACg                                                                                                                                    <  1:212063/150‑1 (MQ=255)
                   ccTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACg                                                                                                                                    <  2:568928/150‑1 (MQ=255)
                   ccTCGCGCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTAGCACCATGACTATGATGTACg                                                                                                                                    >  1:3497478/1‑150 (MQ=255)
                      cgcgCTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGTACGATc                                                                                                                                 >  1:183949/1‑150 (MQ=255)
                        cgcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGt                                                                                                                               <  1:1254338/150‑1 (MQ=255)
                        cgcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGt                                                                                                                               >  2:3074155/1‑150 (MQ=255)
                         gcTGCGCCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGtt                                                                                                                              >  2:1619835/1‑150 (MQ=255)
                             cgcCTCGGCGGAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCa                                                                                                                          <  1:3662612/150‑1 (MQ=255)
                                     cggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGc                                                                                                                  <  2:4889793/150‑1 (MQ=255)
                                     cggAGATGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGc                                                                                                                  <  1:2210035/150‑1 (MQ=255)
                                          aTGAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGttt                                                                                                             >  1:2435388/1‑150 (MQ=255)
                                            gAACGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCt                                                                                                           <  1:4406954/150‑1 (MQ=255)
                                               cGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                        <  1:461270/150‑1 (MQ=255)
                                               cGCCAACAACAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCtggt                                                                                                        <  1:3319887/150‑1 (MQ=255)
                                                   aacaacAGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAAGAATATGATGTACGATCGTTATCAGAAGCCTCTGTTTATGGTGGaa                                                                                                    >  1:3860486/1‑150 (MQ=255)
                                                     caacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaaaa                                                                                                  >  1:5026313/1‑150 (MQ=255)
                                                     caacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCAGGAATATGATGTCCGATCGTTATCAGAAGCCGCTGTTTCTGGTGGaaaa                                                                                                  >  1:5279807/1‑150 (MQ=255)
                                                      aacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                                 >  1:561236/1‑150 (MQ=255)
                                                      aacaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAc                                                                                                 <  1:3542653/150‑1 (MQ=255)
                                                       accgcagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACg                                                                                                <  2:5100430/147‑1 (MQ=255)
                                                       acaGCAGTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACg                                                                                                <  1:470976/150‑1 (MQ=255)
                                                          gcagTGAGGCGAACGTGGTGATATCGCTGCGTAATCCTTATCTACAGGTGAGCTACTGGGGCTTGGGAATGGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGcc                                                                                             <  2:2352575/150‑1 (MQ=255)
                                                           cagTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCt                                                                                            >  2:3030354/1‑150 (MQ=255)
                                                            agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTg                                                                                           >  2:5188289/1‑150 (MQ=255)
                                                            agTGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTg                                                                                           >  2:2119141/1‑150 (MQ=255)
                                                              tGCGGCGAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTggg                                                                                         >  1:5451178/1‑150 (MQ=255)
                                                               gcggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                        <  2:4957995/150‑1 (MQ=255)
                                                               gcggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgc                                                                                        <  1:5156628/150‑1 (MQ=255)
                                                                cggcgAACGTGGTGAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGgcg                                                                                       <  1:823895/150‑1 (MQ=255)
                                                                        gtggtgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGAt                                                                               >  2:536387/1‑150 (MQ=255)
                                                                           gtgAAATCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGaa                                                                            <  1:1170259/150‑1 (MQ=255)
                                                                              aaaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAttt                                                                         >  2:4163100/1‑150 (MQ=255)
                                                                              aaaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAttt                                                                         >  1:3965079/1‑150 (MQ=255)
                                                                               aaTCGCTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTtg                                                                        <  2:5154627/150‑1 (MQ=255)
                                                                                    cTGCGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGcc                                                                   >  2:481460/1‑150 (MQ=255)
                                                                                       cGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                                <  2:4365277/150‑1 (MQ=255)
                                                                                       cGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                                <  2:1302022/150‑1 (MQ=255)
                                                                                       cGTAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGCATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAAt                                                                <  1:255314/150‑1 (MQ=255)
                                                                                         tAATCCGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATgg                                                              <  1:258981/150‑1 (MQ=255)
                                                                                              cGTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCgaga                                                         <  1:2716990/150‑1 (MQ=255)
                                                                                               gTATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAt                                                        >  2:4030078/1‑150 (MQ=255)
                                                                                                tATCTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGAtt                                                       <  1:2996577/150‑1 (MQ=255)
                                                                                                   cTACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAac                                                    >  1:3202926/1‑150 (MQ=255)
                                                                                                    tACAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacg                                                   >  2:3256049/1‑150 (MQ=255)
                                                                                                      cAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac                                                 >  2:1152229/1‑150 (MQ=255)
                                                                                                      cAGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTCTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgac                                                 <  1:1071673/150‑1 (MQ=255)
                                                                                                       aGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                                <  2:5571561/150‑1 (MQ=255)
                                                                                                       aGGTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAacgacg                                                >  1:1347119/1‑150 (MQ=255)
                                                                                                        ggTGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAAcgacga                                               <  1:4108075/150‑1 (MQ=255)
                                                                                                         gTGAGCGACTGGGGCTGGGGAATTGATCAACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCAGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACgacgac                                              >  1:1079818/1‑150 (MQ=255)
                                                                                                          tGAGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACt                                             >  1:3700437/1‑150 (MQ=255)
                                                                                                            aGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                           >  1:5093498/1‑150 (MQ=255)
                                                                                                            aGCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTAt                                           >  1:1818949/1‑150 (MQ=255)
                                                                                                             gCGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATc                                          <  2:3639135/150‑1 (MQ=255)
                                                                                                              cGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg                                         >  2:806853/1‑150 (MQ=255)
                                                                                                              cGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg                                         >  1:2323992/1‑150 (MQ=255)
                                                                                                              cGACTGGGGCTGGGGAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCg                                         >  2:52996/1‑150 (MQ=255)
                                                                                                                  tggggctggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                     >  1:1601088/1‑150 (MQ=255)
                                                                                                                  tggggctggggAATTGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                     >  1:2919711/1‑150 (MQ=255)
                                                                                                                  tggggctggggAAGTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTCCGCTCGTGATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATc                                     <  2:2950710/150‑1 (MQ=255)
                                                                                                                        tggggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTAc                               <  2:2708247/150‑1 (MQ=255)
                                                                                                                          gggAATTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAATAGCTGAATTTGCTGCCAAGGGCGAGCTTAACGACGACTCTCGCCTCAGCTACtt                             >  2:3965988/1‑150 (MQ=255)
                                                                                                                              aTTGATCCACTCGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTAcgc                         <  1:4245140/150‑1 (MQ=255)
                                                                                                                                tGATCCACTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGa                       >  2:4617614/1‑150 (MQ=255)
                                                                                                                                      aCTCGGTCTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACatat                 >  1:4845235/1‑150 (MQ=255)
                                                                                                                                        tCGGTCTTCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATcc               <  2:1919937/150‑1 (MQ=255)
                                                                                                                                         cGGTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcg              >  2:4333202/1‑150 (MQ=255)
                                                                                                                                          ngTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgc             >  2:2650098/2‑150 (MQ=255)
                                                                                                                                          ggTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgc             >  2:705446/1‑150 (MQ=255)
                                                                                                                                          ggTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgc             >  2:1383889/1‑150 (MQ=255)
                                                                                                                                           gTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgcg            >  2:846061/1‑150 (MQ=255)
                                                                                                                                           gTCTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCcgcg            <  2:3539806/150‑1 (MQ=255)
                                                                                                                                             cTGCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCa          >  2:3761704/1‑150 (MQ=255)
                                                                                                                                             cTGCGTATAACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCa          >  1:5662503/1‑150 (MQ=255)
                                                                                                                                               gCGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAAAATGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAAt        >  2:4131465/1‑150 (MQ=255)
                                                                                                                                                cGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg       <  1:4364526/150‑1 (MQ=255)
                                                                                                                                                cGTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg       <  1:536387/150‑1 (MQ=255)
                                                                                                                                                 gTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATgg      >  1:4944254/1‑150 (MQ=255)
                                                                                                                                                 gTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATgg      <  2:2435388/150‑1 (MQ=255)
                                                                                                                                                  tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg     >  2:3911948/1‑150 (MQ=255)
                                                                                                                                                  tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg     >  2:1502397/1‑150 (MQ=255)
                                                                                                                                                  tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg     >  2:4210597/1‑150 (MQ=255)
                                                                                                                                                   aTCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGc    >  1:2820526/1‑150 (MQ=255)
                                                                                                                                                   aTCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGc    >  2:2591502/1‑150 (MQ=255)
                                                                                                                                                     cACCATGAATATGATGTACGATCGTTATCTGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGa  <  2:4487441/150‑1 (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: