Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,896,189 G→T 6.8% T291T (ACC→ACA bglB ← cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,896,1890GT 6.8% 859.5 / 28.3 227T291T (ACC→ACAbglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (138/73);  new base T (10/6);  total (148/79)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

CAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAA  >  CP009273/3896041‑3896338
                                                                                                                                                     |                                                                                                                                                     
cAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACgg                                                                                                                                                      >  1:739995/1‑150 (MQ=255)
cAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACgg                                                                                                                                                      >  1:993357/1‑150 (MQ=255)
    cGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATg                                                                                                                                                  >  1:3165823/1‑150 (MQ=255)
     ggATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGt                                                                                                                                                 >  1:2345774/1‑150 (MQ=255)
     ggATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGt                                                                                                                                                 >  1:2103685/1‑150 (MQ=255)
         cAATTTGCCACCCCCACTCTGAACTTTTCAGTTGCGGATTGGGGATAATATTCTGTATGTTGCCCTGCGAATTTTTATTACTGCTTTCGTCGTGGGAAACAAAACCAGTCATGTAATAACTAAAAGAGATGAACTCGACTGTATGTTTTa                                                                                                                                             <  2:2350877/150‑1 (MQ=255)
          aaTTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTaa                                                                                                                                            >  2:2827001/1‑150 (MQ=255)
           aTTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGCTCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTaaa                                                                                                                                           >  2:1155642/1‑150 (MQ=255)
           aTTTCCCACCCCCACTCTGAACTTTTCATATGCGGTTTGGTTATAATATTCATTTTTTTGCCCTGCTATTTTTTATTAATTCTGTCGTCGTGGGAAACACACCCTGTCATGTAATAAATAAAAGAGATTAAATCGACGGTATGTTTTaaa                                                                                                                                           <  2:147777/150‑1 (MQ=255)
           aTTGGCCACCCCCACTCTTAACTTTTCAGTTGCGGTTTGGGTATCATATTCAGTATGTTGCCCTGCGAATTTTTATTAATGCTTTCGTGGTGGTAAACACAACCAGTCATTTAATAACTAAAAGAGATCAACTCGACGGTATGTTTTaaa                                                                                                                                           <  2:145792/150‑1 (MQ=255)
                ccACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttc                                                                                                                                      >  1:585148/1‑150 (MQ=255)
                 cACCCCCACTCTTAACTTTTCAGATGCGGTTTGGGGATCATATTCATTATGTTGCCCTGCGAATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATTTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttct                                                                                                                                     <  2:1934569/150‑1 (MQ=255)
                 cACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttct                                                                                                                                     >  1:1615557/1‑150 (MQ=255)
                 cACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttct                                                                                                                                     >  1:4351605/1‑150 (MQ=255)
                    ccccACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCg                                                                                                                                  >  2:3719318/1‑150 (MQ=255)
                           ctGAACTTTTCTGTTGGGTTTTGGGGATCATTTTCTGTTTGTTGCCCGGCGTTTTTTTTTTCTTGTTTTGGGCGTGGTAAACACACCCCGTCTTTTACTAACTAACAGAGTTTAACTCGACGGTATGTTTTAAATCTTCTGCGTCACttt                                                                                                                           <  2:26961/150‑1 (MQ=255)
                                 ttttCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCa                                                                                                                     >  2:1657320/1‑150 (MQ=255)
                                  tttCTGATGCGGTTTGGGTATAATTTTCATTTTGTTGCCCTGCGAATTTTTATTATTGCTTTCGTGGTGGTAAACACAACCAGTCATTTAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCAt                                                                                                                    <  2:1989263/150‑1 (MQ=255)
                                       gttGCGGTTTGGGTATCATTTTCTGTTTGTTGCCCTGCGATTTTTTATTATGGCTTTCGTCGTGGTAAACACAACCAGTCATGTAATAACTAAAAGAGAGGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCaa                                                                                                               <  2:4344337/148‑1 (MQ=255)
                                           cGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGt                                                                                                           >  1:1011103/1‑150 (MQ=255)
                                             ttttGGGTATCATATTCTTTTTTTTGCCGTGCTATTTTTTATTATTGTTTTCGTGGTGGTAAAAACAACCAGAAATTTAATAAATAAATTTTTTTAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGa                                                                                                         <  1:257674/148‑1 (MQ=255)
                                              aTTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        >  2:86968/1‑150 (MQ=255)
                                              aTTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        >  1:395391/1‑150 (MQ=255)
                                               ttGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGata                                                                                                       >  2:180244/1‑150 (MQ=255)
                                                  ggcataatatTCAGTATTTTCCCCTGCGAATTTTTATTATTCCTTTCGTCGTGGTAAACAAAACAATAAATTTAATAAATAAAAGATATCAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATtg                                                                                                    <  1:1948486/144‑1 (MQ=255)
                                                         tatTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTcgc                                                                                             >  2:2873870/1‑150 (MQ=255)
                                                         tatTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTcgc                                                                                             >  1:2833963/1‑150 (MQ=255)
                                                           ttCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCgg                                                                                           >  1:2700778/1‑150 (MQ=255)
                                                             cAGTTTTTTGCCCTGCGAATTTTTATAAATGCTTGCGGCGTGGTAAACACAACCAGGCATGTAAAAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaa                                                                                         <  2:4108478/150‑1 (MQ=255)
                                                             cAGTAGGGTGCCCTGCGAATTTTTATTAATGCTTGCGTCGTGGGAAACACAACCAGTCATTTAATAACTAAAAGAGAGGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaa                                                                                         <  2:1746384/150‑1 (MQ=255)
                                                                tttggtGCCGTGCGATTTTTTTTAATGGTTTGCGTGTTGGTAAACACACCAAGACTTGAAATAACAAAAACTTATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaagaa                                                                                      <  1:2955363/145‑1 (MQ=255)
                                                                  ttttGCCCTGCTATTTTTTATTATTGCTTTCGTCGTGGTAAACACAACCAGTCATGTAATAAATAAAAGAGATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAAc                                                                                    <  2:2483802/148‑1 (MQ=255)
                                                                  tGTTGCCCTGCTAATTTTTATTAATGCTTTCGGCGTGGTAAACACAACCAGTCATTAAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAAc                                                                                    <  2:2955262/150‑1 (MQ=255)
                                                                  tGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAAc                                                                                    >  2:2046748/1‑150 (MQ=255)
                                                                  tGTTGCCCTGCGATTTTTTATTAAGGCTTTCGTCGTGGGAAACACAACCAGTCATTTAAAAAATAAAAGATATTAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAAc                                                                                    <  1:1155642/150‑1 (MQ=255)
                                                                   tttGCCCTGCTAATTTTTATTATTGCTTTGGTGGTGGAAAACACAACCAGCATTTTAATAACAAAACCATATCAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACg                                                                                   <  1:3733873/149‑1 (MQ=255)
                                                                   gTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACg                                                                                   >  2:2127202/1‑150 (MQ=255)
                                                                   gTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACg                                                                                   >  2:1074335/1‑150 (MQ=255)
                                                                     tGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGACGAAACGCt                                                                                 >  1:1882810/1‑150 (MQ=255)
                                                                        ccTGCGTTTTTTTATTATGGCTTGCGGGGTGGCAAACACACCCCGTCATTAAAAAACAAAAAGACAACAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCa                                                                              <  1:4006743/150‑1 (MQ=255)
                                                                        ccTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCa                                                                              >  2:424597/1‑150 (MQ=255)
                                                                          tGCGCATTTTTATTAATGCTTTCGTCGTGGTAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCata                                                                            <  2:789472/150‑1 (MQ=255)
                                                                           gcgcATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCatat                                                                           >  1:2884416/1‑150 (MQ=255)
                                                                            cgcATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGACGAAACGCTGCAtata                                                                          >  2:3065707/1‑150 (MQ=255)
                                                                              atttttttttAATGCTTTCTTCGTGGTAAACACAACCAGTAATGTAATAACTAAAAGATATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGc                                                                        <  2:699151/148‑1 (MQ=255)
                                                                               aTTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGcc                                                                       >  1:4242209/1‑150 (MQ=255)
                                                                               aTTTTTATTAATGCTTACGTCGTGTTAAACACAACCAGACATTTAATAACAAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGcc                                                                       <  2:2594718/150‑1 (MQ=255)
                                                                                tttttATTAATTCTTTCGTCGTGTGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGAATGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGTTGATATTGTGGTCTCGGAAGAAACGCTGCATATAGCCg                                                                      >  1:3588928/1‑150 (MQ=255)
                                                                                tttttATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCg                                                                      >  1:883970/1‑150 (MQ=255)
                                                                                 ttttATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCgg                                                                     >  2:436941/1‑150 (MQ=255)
                                                                                    tattaATGCTTTCGTCGTGGGAAACACAACCCGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGAt                                                                  <  2:3143893/150‑1 (MQ=255)
                                                                                    tattaATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGAt                                                                  >  1:3678952/1‑150 (MQ=255)
                                                                                     attaATGCTTTCGTCGTGGTAAACACAACCAGTCATTTAATAACTAAAACAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATa                                                                 <  2:4614138/150‑1 (MQ=255)
                                                                                         aTGCGTTCGTCGGGGGAAAAACAACCAGACATTTAATAACTAAAATCTATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTgg                                                             <  2:1367904/150‑1 (MQ=255)
                                                                                           gCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGcc                                                           >  1:541308/1‑150 (MQ=255)
                                                                                              ttCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcg                                                        >  2:1492382/1‑150 (MQ=255)
                                                                                              ttCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcg                                                        >  1:2073410/1‑150 (MQ=255)
                                                                                              ttCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcg                                                        >  1:3552872/1‑150 (MQ=255)
                                                                                               tcgtcgTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcgc                                                       >  1:2776709/1‑150 (MQ=255)
                                                                                                tcgtgttGGTAAACACAACCTGTCATGTAATAACTAAAAGATATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcgc                                                       <  2:4763174/144‑1 (MQ=255)
                                                                                                      tGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGa                                                 >  2:3688055/1‑150 (MQ=255)
                                                                                                       gggAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGaa                                                >  2:4469586/1‑150 (MQ=255)
                                                                                                        ggAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAAc                                               >  2:250128/1‑150 (MQ=255)
                                                                                                         taaaCACAACCAGTCATTTAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACa                                              <  2:2068814/149‑1 (MQ=255)
                                                                                                          aaaaaaaaCCAGTCATTTAATAACAAAAAGACATCAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACAt                                             <  1:4612702/144‑1 (MQ=255)
                                                                                                          aaaCACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACAt                                             >  1:3203436/1‑150 (MQ=255)
                                                                                                               caacaaGTAATTTAAAAAAAAAAAGATATCAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCa                                        <  1:3011656/150‑1 (MQ=255)
                                                                                                                aaCCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaa                                       >  2:409207/1‑150 (MQ=255)
                                                                                                                 ccctgTCTTGTAATAAAAAAATTAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaaa                                      <  1:4774323/146‑1 (MQ=255)
                                                                                                                 aCCTGTCTTGTAATAAATAAATGATTTTAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaaa                                      <  1:1475797/150‑1 (MQ=255)
                                                                                                                  ccAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAg                                     >  2:4116295/1‑150 (MQ=255)
                                                                                                                   cAGTCTTGTATTAAATAAAAGTTATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGa                                    <  2:739931/150‑1 (MQ=255)
                                                                                                                   cAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGTCGGGATACTGGCCACGCGCCTGAACATCACCAAAGa                                    >  2:2266477/1‑150 (MQ=255)
                                                                                                                        ttttaataaCTAAAAGATATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATc                               <  2:3552872/147‑1 (MQ=255)
                                                                                                                        aTGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATc                               >  2:705677/1‑150 (MQ=255)
                                                                                                                         tGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATcc                              >  1:1910824/1‑150 (MQ=255)
                                                                                                                           taataaCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAg                            >  1:3209962/1‑150 (MQ=255)
                                                                                                                           taataaATAAATGATATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAg                            <  1:2541080/150‑1 (MQ=255)
                                                                                                                             ataaCTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAgcg                          >  2:4658434/1‑150 (MQ=255)
                                                                                                                               aaCTAAAATTGTTGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGcc                        <  2:3498195/150‑1 (MQ=255)
                                                                                                                                   aaaaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGtt                    >  1:3798984/1‑150 (MQ=255)
                                                                                                                                     aaGAGATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTtct                  <  2:311038/150‑1 (MQ=255)
                                                                                                                                     aaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAg                            >  1:1841550/1‑140 (MQ=255)
                                                                                                                                      agagATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTtctc                 >  1:978498/1‑150 (MQ=255)
                                                                                                                                       gagaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTctct                >  1:3384295/1‑150 (MQ=255)
                                                                                                                                        tgaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               <  1:808381/149‑1 (MQ=255)
                                                                                                                                        atattaaaTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               <  1:807987/145‑1 (MQ=255)
                                                                                                                                        agaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               >  2:274449/1‑150 (MQ=255)
                                                                                                                                        agaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               >  2:3313735/1‑150 (MQ=255)
                                                                                                                                        agaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               >  2:3021797/1‑150 (MQ=255)
                                                                                                                                          aTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTcc             >  2:1303257/1‑150 (MQ=255)
                                                                                                                                           tGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGACACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCa            >  2:4543219/1‑150 (MQ=255)
                                                                                                                                            gAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCAt           <  2:585148/150‑1 (MQ=255)
                                                                                                                                             aaaTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATg          >  2:1112149/1‑150 (MQ=255)
                                                                                                                                               aTCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGc        >  1:671455/1‑150 (MQ=255)
                                                                                                                                                tCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGcc       <  1:4644918/150‑1 (MQ=255)
                                                                                                                                                 cGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCt      >  2:4029906/1‑150 (MQ=255)
                                                                                                                                                 cGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCt      <  2:3362552/150‑1 (MQ=255)
                                                                                                                                                    cGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCa   >  2:1828038/1‑150 (MQ=255)
                                                                                                                                                     ggTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCaa  <  2:1223178/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAA  >  CP009273/3896041‑3896338

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