Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,896,189 G→T 14.9% 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,1890GT14.9% 879.7 / 93.8 274T291T (ACC→ACAbglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (113/120);  new base T (27/14);  total (140/134)
Fisher's exact test for biased strand distribution p-value = 4.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

CCAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAA  >  CP009273/3896040‑3896338
                                                                                                                                                     |                                                                                                                                                     
ccAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACg                                                                                                                                                       <  1:3372622/150‑1 (MQ=255)
   aCCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTa                                                                                                                                                    >  1:2365796/1‑150 (MQ=255)
    ccGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTAt                                                                                                                                                   >  2:4062543/1‑150 (MQ=255)
    ccGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTAt                                                                                                                                                   >  2:1520819/1‑150 (MQ=255)
          cAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTa                                                                                                                                             >  2:273774/1‑150 (MQ=255)
          cAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTa                                                                                                                                             >  1:3713384/1‑150 (MQ=255)
          cAATTGGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATTTTCATTATGTTGCCCTGCGCTTTTTTATTAATGATTTCGTCGTGGTAAACACAACCAGTCATGTAATAAATAAAAGAGATTAAATCGACTGTATGTTTTa                                                                                                                                             <  1:4250315/150‑1 (MQ=255)
           attttGCCACCCCCACTGGGAACTTTTCAGATGCGGATTGGGGATCATATTCTGTATGTGGCCCTGCGATTTTTTATTAATGCTTTGGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTaa                                                                                                                                            <  1:1235815/148‑1 (MQ=255)
           aaTTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTaa                                                                                                                                            >  1:2951691/1‑150 (MQ=255)
           aaTTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTaa                                                                                                                                            >  1:1216837/1‑150 (MQ=255)
           aaTTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTaa                                                                                                                                            >  2:400519/1‑150 (MQ=255)
                gCCACCCCCACTCTGAACTTTTCAGATGCGGTTTGGGTATCATATTCATTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATTTAATAACTAAAACAGATGAAATCGACGGTATGTTTTAAAtctt                                                                                                                                       <  2:266659/150‑1 (MQ=255)
                gCCACCCCCACACTGAACTTTTCAGATGCGGATTGGGGATAATATTCAGTATTTTGCCCTGCGCATTTTTATTAATGATTTCGTCGTGGTAAACAAAACCAGTCATGTAATAACTAAAAGAGATCAACTCGACGGTATGTTTTAAAtctt                                                                                                                                       <  1:5146129/150‑1 (MQ=255)
                 ccACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttc                                                                                                                                      >  1:2075862/1‑150 (MQ=255)
                  cACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAAtcttct                                                                                                                                     <  2:3435851/150‑1 (MQ=255)
                       ccACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTc                                                                                                                                >  2:4089793/1‑150 (MQ=255)
                          ctctTAACTTTTCTGATGCGTATTGGGGATCATTTTCAGTATGTGGCCCTGCGAATTTTTTTTAATGCTTTCGTCGTGGGAAACACAACCAGTCATTTAATAACTAAAAGAGATTAACTCGACTGTATGTTTTAAATCTTCTGCGTCACt                                                                                                                             <  1:4558659/150‑1 (MQ=255)
                           tctGAACATTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACtt                                                                                                                            >  2:3358735/1‑150 (MQ=255)
                            ctcaaCTTTTCTGTTGGGGTTTGGGTATAATATTAATTATGTTGCCCTGCGAATTTTTATTAATGCTTTCGTCGTGTGAAACAAAACCAGTCATGTAATAAATAAAAGATATCAAATCGACGGTATGTTTTAAATCTTCTGCGTCACttt                                                                                                                           <  1:2741323/147‑1 (MQ=255)
                            ctGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACttt                                                                                                                           >  1:2725368/1‑150 (MQ=255)
                             tGAACTTTTCAGATGCGGAGTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTc                                                                                                                          >  1:82344/1‑150 (MQ=255)
                                 cTTTTCTGAGGCGTATTGGGTATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCCGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTc                                                                                                                      <  2:2333066/150‑1 (MQ=255)
                                  ttttCAGATGCGGATTGGGGATCATTTTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCCGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCa                                                                                                                     <  2:320094/150‑1 (MQ=255)
                                  ttttCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCa                                                                                                                     >  2:3274499/1‑150 (MQ=255)
                                    ttCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCAtc                                                                                                                   >  1:2237868/1‑150 (MQ=255)
                                       aGATGCGTTTTGGGTATAATATTCTGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGTAAACACAACCAGTCATGTAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCa                                                                                                                <  1:399991/150‑1 (MQ=255)
                                        gATGCGGATTGGGGAACATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCTTCGTGGGAAACACAACCCGTCATGTAATAACTAAAACACATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCaa                                                                                                               <  2:4286021/150‑1 (MQ=255)
                                            cgtttggggtATCATATTCTGTATGTTGCCCGGCGCATTTTTATTAAGGCTTTGGCGGTGGTAAACACAACCAGTCATTTAATAAATAAAACAGATCAACTCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGt                                                                                                           <  1:4341613/144‑1 (MQ=255)
                                             ggATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTg                                                                                                          >  1:2337185/1‑150 (MQ=255)
                                             ggATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTg                                                                                                          >  1:1340071/1‑150 (MQ=255)
                                              gATTGGGGATCATATTCTGTATGTTGCCCTGCGAATTTTTATTAATGCTTTCGTGGTGGGAAACACAACCAGTCATGTAATAACTAAAAGATATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGa                                                                                                         <  1:3652692/150‑1 (MQ=255)
                                               tttGGGTATAATATTCTGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        <  1:3062160/149‑1 (MQ=255)
                                               tttGGGGATCATATTCCGTATGTTGCCCGGCGCATTTTTATTAATGCTTTCGTCGTGGTAAACACAACCAGTCTTGTAATAACAAAAATAGATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        <  1:3542641/149‑1 (MQ=255)
                                               aTTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        >  2:92849/1‑150 (MQ=255)
                                               aTTGGGGAAAATATTCTGTATGTTGCCCTGCGAATTTTTATTAATGCTTTCGTCGTGGTAAACAAAACCAGTCATGTAATAACTAAAAGACATCAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGat                                                                                                        <  1:2808650/150‑1 (MQ=255)
                                                ttGGGTATCATATTCAGTATGTTGCCCTGCGAATTTTTATTAATGATTTCGTCGTGGGAAACACAACCAGTCATGTAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGata                                                                                                       <  1:374234/150‑1 (MQ=255)
                                                  ggggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGatagt                                                                                                     >  2:4707400/1‑148 (MQ=255)
                                                  ggggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATAtt                                                                                                     >  1:1371092/1‑150 (MQ=255)
                                                   gggATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATtg                                                                                                    >  2:2140720/1‑150 (MQ=255)
                                                      aTCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTgg                                                                                                 >  2:4786910/1‑150 (MQ=255)
                                                              cAGTTTGTTGCCCTGCGCATTTTTATAAATGCTTTCGTCGTGGTAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaa                                                                                         <  1:1459042/150‑1 (MQ=255)
                                                               aGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaag                                                                                        >  1:3241582/1‑150 (MQ=255)
                                                                gTATGTTGCCCTGCGCTTTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaaga                                                                                       <  1:4133165/150‑1 (MQ=255)
                                                                 tagtttGCCCGGCGCATTTTTATTAAGGATTTCGTCGTGGGAAACACAACCAGTCATTTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGgaagaa                                                                                      <  2:1030617/146‑1 (MQ=255)
                                                                   tGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAAc                                                                                    >  2:545643/1‑150 (MQ=255)
                                                                       gCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTg                                                                                <  2:390803/150‑1 (MQ=255)
                                                                        cccTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGc                                                                               >  2:1960187/1‑150 (MQ=255)
                                                                           tGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCata                                                                            >  2:2951432/1‑150 (MQ=255)
                                                                               cATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGc                                                                        >  2:412733/1‑150 (MQ=255)
                                                                               cATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGc                                                                        >  2:2263931/1‑150 (MQ=255)
                                                                               cATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGc                                                                        >  1:646045/1‑150 (MQ=255)
                                                                                  ttttATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCgg                                                                     >  2:4984136/1‑150 (MQ=255)
                                                                                   tttATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCggg                                                                    >  1:4797385/1‑150 (MQ=255)
                                                                                         aaTGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTg                                                              >  2:721577/1‑150 (MQ=255)
                                                                                            gCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGcc                                                           >  1:3495731/1‑150 (MQ=255)
                                                                                              tttGGTGGTGGGAAACACAACCTGTCTTGTAATAACTAAATGTGTTGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAc                                                         <  2:3118331/150‑1 (MQ=255)
                                                                                              tttCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAc                                                         <  2:48349/150‑1 (MQ=255)
                                                                                                tcgtcgTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCAcgc                                                       <  1:1992495/150‑1 (MQ=255)
                                                                                                  gtcgtGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACgcgc                                                     >  1:4953291/1‑150 (MQ=255)
                                                                                                  gtcgtGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACgcgc                                                     >  1:1094982/1‑150 (MQ=255)
                                                                                                  gtcgtGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACgcgc                                                     >  2:4921283/1‑150 (MQ=255)
                                                                                                   tcgtGGTAAACACAACCAGTCTTTAAATAAATAAAAGAGATTAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGcc                                                    <  1:437364/150‑1 (MQ=255)
                                                                                                    tgtGGGAAACACACCACGTCTTTTAATAACTAAAAGAGTTGATATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAATAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCt                                                   <  2:3548434/149‑1 (MQ=255)
                                                                                                    cgtGGGAAACACAACCAGTCATGTAATAACTAAATGAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCt                                                   <  2:2996650/150‑1 (MQ=255)
                                                                                                     ttGGTAAACACAACCCGTCATTTAATAACTAAAAGAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTg                                                  <  2:983606/149‑1 (MQ=255)
                                                                                                      tGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGa                                                 >  2:3781906/1‑150 (MQ=255)
                                                                                                        ggAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAAc                                               >  2:4612975/1‑150 (MQ=255)
                                                                                                              acaACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCAcc                                         <  2:1130903/150‑1 (MQ=255)
                                                                                                                aaCCAGTCATGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaa                                       >  2:3841929/1‑150 (MQ=255)
                                                                                                                aaCCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCaa                                       >  1:1162993/1‑150 (MQ=255)
                                                                                                                  ccAGTCTTGTAATAACTAAAAGAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAg                                     <  1:1563969/150‑1 (MQ=255)
                                                                                                                  ccAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAg                                     <  1:2099725/150‑1 (MQ=255)
                                                                                                                    aGTCTTGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGaa                                   <  2:4218193/150‑1 (MQ=255)
                                                                                                                      tCATGTAATAAATAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACa                                 <  2:469685/150‑1 (MQ=255)
                                                                                                                        aTGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATc                               >  1:4030877/1‑150 (MQ=255)
                                                                                                                         tttaataaCTAAAAGTGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATcc                              <  1:1848196/148‑1 (MQ=255)
                                                                                                                         tGTAATAACTAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATcc                              >  1:3899672/1‑150 (MQ=255)
                                                                                                                         tGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATcc                              >  2:1911708/1‑150 (MQ=255)
                                                                                                                                 atAAAAGAGATGAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCgg                      <  1:4489337/149‑1 (MQ=255)
                                                                                                                                  tAAAAGAGATGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGt                     >  1:1012475/1‑150 (MQ=255)
                                                                                                                                   aaaaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGtt                    >  1:503658/1‑150 (MQ=255)
                                                                                                                                   aaaaGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGtt                    >  1:963956/1‑150 (MQ=255)
                                                                                                                                      agagATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTtctc                 >  2:5117288/1‑150 (MQ=255)
                                                                                                                                       gagaTTAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTctct                <  2:4485677/150‑1 (MQ=255)
                                                                                                                                       gagaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTctct                >  1:719724/1‑150 (MQ=255)
                                                                                                                                        agaTGAAATCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCtt               >  1:2490713/1‑150 (MQ=255)
                                                                                                                                         gaTTAACTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTc              <  1:3082364/150‑1 (MQ=255)
                                                                                                                                         gaTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTc              <  2:2431741/150‑1 (MQ=255)
                                                                                                                                          aTGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTcc             <  1:3120501/150‑1 (MQ=255)
                                                                                                                                           tGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCa            >  1:4841873/1‑150 (MQ=255)
                                                                                                                                            gAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCAt           >  2:2791217/1‑150 (MQ=255)
                                                                                                                                            gAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCAt           <  1:1160366/150‑1 (MQ=255)
                                                                                                                                              aaTCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATgg         <  1:338559/150‑1 (MQ=255)
                                                                                                                                               aTCGACTGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGc        >  1:2684319/1‑150 (MQ=255)
                                                                                                                                                 cGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCt      <  1:161742/150‑1 (MQ=255)
                                                                                                                                                  gACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTg     >  1:2373944/1‑150 (MQ=255)
                                                                                                                                                  gACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTg     >  1:2047901/1‑150 (MQ=255)
                                                                                                                                                  gACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTg     >  2:1625871/1‑150 (MQ=255)
                                                                                                                                                    cGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCa   >  2:21136/1‑150 (MQ=255)
                                                                                                                                                     ggTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCaa  >  1:4141897/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAACCGGATCAATTTGCCACCCCCACTCTGAACTTTTCAGATGCGGATTGGGGATCATATTCAGTATGTTGCCCTGCGCATTTTTATTAATGCTTTCGTCGTGGGAAACACAACCAGTCATGTAATAACTAAAAGAGATGAAATCGACGGTATGTTTTAAATCTTCTGCGTCACTTTCAGTCATCTCAATGGTGATATTGTGGTCGCGGAAGAAACGCTGCATATAGCCGGGATACTGGCCACGCGCCTGAACATCACCAAAGAACATCCAGCGCCGGTTCTCTTCCATGGCCTGCAA  >  CP009273/3896040‑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: