Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,300,333 T→A 37.1% E155D (GAA→GAT alsA ← fused D‑allose transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,300,3330TA37.1% 227.4 / 283.2 224E155D (GAA→GATalsAfused D‑allose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (67/74);  new base A (34/49);  total (101/123)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 229 total aligned reads displayed.

AACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300184‑4300482
                                                                                                                                                     |                                                                                                                                                     
 aCTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAttt                                                                                                                                                      <  2:2501912/150‑1 (MQ=255)
  cTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTc                                                                                                                                                     <  2:2883281/150‑1 (MQ=255)
  cTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTc                                                                                                                                                     <  1:1975664/150‑1 (MQ=255)
   ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  1:2603010/150‑1 (MQ=255)
    tATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTa                                                                                                                                                   <  2:3910614/150‑1 (MQ=255)
    tATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTa                                                                                                                                                   <  2:3911013/150‑1 (MQ=255)
             tataGACGATGGTCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCt                                                                                                                                          <  1:1564365/150‑1 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          >  1:3146365/1‑150 (MQ=255)
                 gACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgtg                                                                                                                                      <  1:2690416/150‑1 (MQ=255)
                  aCGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTgg                                                                                                                                     >  2:1301825/1‑150 (MQ=255)
                   cGATGGCCATACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGc                                                                                                                                    <  2:3899394/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   >  1:2023672/1‑150 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   <  2:2329157/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   <  2:2329619/150‑1 (MQ=255)
                       ggCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                                                                                                                                <  1:290913/150‑1 (MQ=255)
                        gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgata                                                                                                                               >  2:112402/1‑150 (MQ=255)
                          cGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgataga                                                                                                                             <  1:325890/150‑1 (MQ=255)
                           gTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagat                                                                                                                            <  2:3712945/150‑1 (MQ=255)
                            tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           >  1:1797465/1‑150 (MQ=255)
                            tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           <  1:1054984/150‑1 (MQ=255)
                            tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           <  2:20209/150‑1 (MQ=255)
                            tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           <  2:83244/150‑1 (MQ=255)
                                    ttACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGcc                                                                                                                   <  1:3696534/150‑1 (MQ=255)
                                      aCGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCAc                                                                                                                 >  2:1282771/1‑150 (MQ=255)
                                       cGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACt                                                                                                                >  1:1822949/1‑150 (MQ=255)
                                        gCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCCAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt                                                                                                               <  2:3313477/150‑1 (MQ=255)
                                         cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              >  1:155801/1‑150 (MQ=255)
                                            cTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtctc                                                                                                           >  1:885972/1‑150 (MQ=255)
                                             tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          <  1:360039/150‑1 (MQ=255)
                                              gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         <  1:2258627/150‑1 (MQ=255)
                                                ttCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCt                                                                                                       <  2:1765421/150‑1 (MQ=255)
                                                    tGATCAGAAACAGATGGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGAt                                                                                                   <  2:696670/150‑1 (MQ=255)
                                                     gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATc                                                                                                  >  1:3135108/1‑150 (MQ=255)
                                                     gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATc                                                                                                  >  1:51442/1‑150 (MQ=255)
                                                      aTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCa                                                                                                 <  2:4736999/150‑1 (MQ=255)
                                                       tCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCaa                                                                                                >  2:3017785/1‑150 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               <  1:2556431/150‑1 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  2:1597497/1‑150 (MQ=255)
                                                            aaCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTc                                                                                           >  1:1486590/1‑150 (MQ=255)
                                                             aCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCa                                                                                          <  1:467004/150‑1 (MQ=255)
                                                              cAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCaa                                                                                         <  1:295819/150‑1 (MQ=255)
                                                               aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAg                                                                                        <  1:2359521/150‑1 (MQ=255)
                                                                gATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGc                                                                                       <  2:433894/150‑1 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     >  2:4708200/1‑150 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     <  1:1597526/150‑1 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     <  2:150916/150‑1 (MQ=255)
                                                                   aGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCa                                                                                    >  2:2196534/1‑150 (MQ=255)
                                                                   aGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCa                                                                                    >  2:4597238/1‑150 (MQ=255)
                                                                   aGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCa                                                                                    >  1:214885/1‑150 (MQ=255)
                                                                    gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCGCTTTCTCATCTAGATCAACTTTCAAGCCCAc                                                                                   >  2:4785264/1‑150 (MQ=255)
                                                                    gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAc                                                                                   <  1:4462003/150‑1 (MQ=255)
                                                                     tCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcg                                                                                  >  2:2704808/1‑150 (MQ=255)
                                                                      ccACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgc                                                                                 <  1:3314449/150‑1 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                >  2:762794/1‑150 (MQ=255)
                                                                        aCCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGt                                                                               >  2:3054249/1‑150 (MQ=255)
                                                                            ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           <  1:2477794/150‑1 (MQ=255)
                                                                             tttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaataa                                                                          <  1:1696474/150‑1 (MQ=255)
                                                                               tATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAca                                                                        >  2:2577226/1‑150 (MQ=255)
                                                                               tATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAca                                                                        >  1:4362642/1‑150 (MQ=255)
                                                                                aTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcat                                                                       <  1:4691035/150‑1 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   >  2:2975269/1‑150 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   <  2:4305585/150‑1 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   <  1:2503486/150‑1 (MQ=255)
                                                                                        gtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcgg                                                               <  2:2280304/150‑1 (MQ=255)
                                                                                         tgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcggc                                                              >  2:3740640/1‑150 (MQ=255)
                                                                                         tgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcggc                                                              >  1:1184839/1‑150 (MQ=255)
                                                                                         tgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcggc                                                              >  2:3487225/1‑150 (MQ=255)
                                                                                           aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            <  2:1109244/150‑1 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           <  1:2811001/150‑1 (MQ=255)
                                                                                             gaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCgg                                                          <  2:3852831/150‑1 (MQ=255)
                                                                                               ggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGac                                                        <  2:3008444/150‑1 (MQ=255)
                                                                                                     ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAt                                                  <  1:1624402/150‑1 (MQ=255)
                                                                                                     ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAt                                                  <  1:2973933/150‑1 (MQ=255)
                                                                                                     ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAt                                                  <  1:2973046/150‑1 (MQ=255)
                                                                                                        aTCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtc                                               >  2:2441790/1‑150 (MQ=255)
                                                                                                        aTCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtc                                               >  2:2440192/1‑150 (MQ=255)
                                                                                                         tCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTCGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtct                                              <  2:45666/150‑1 (MQ=255)
                                                                                                          ccATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtctc                                             <  1:3472011/150‑1 (MQ=255)
                                                                                                            atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           <  2:1797465/150‑1 (MQ=255)
                                                                                                                  atgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCg                                     <  1:2354962/150‑1 (MQ=255)
                                                                                                                    gaCTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGat                                   >  2:1067282/1‑150 (MQ=255)
                                                                                                                     aCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGata                                  <  2:1177498/150‑1 (MQ=255)
                                                                                                                      cTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGataa                                 <  1:112402/150‑1 (MQ=255)
                                                                                                                        ttGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGataata                               >  2:1955327/1‑150 (MQ=255)
                                                                                                                        ttGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGataata                               <  2:4148753/150‑1 (MQ=255)
                                                                                                                         tGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAatat                              >  1:2943339/1‑150 (MQ=255)
                                                                                                                          ggCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATAtt                             >  2:1330872/1‑150 (MQ=255)
                                                                                                                           gCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTg                            <  2:1565642/150‑1 (MQ=255)
                                                                                                                            cATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGa                           <  1:2858971/150‑1 (MQ=255)
                                                                                                                            cATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGa                           <  1:434745/150‑1 (MQ=255)
                                                                                                                              tCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACg                         <  1:4123901/150‑1 (MQ=255)
                                                                                                                                 aGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCca                      <  1:3391746/150‑1 (MQ=255)
                                                                                                                                 aGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCca                      <  1:2337/150‑1 (MQ=255)
                                                                                                                                  gCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCcac                     >  1:1303338/1‑150 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              >  1:396195/1‑150 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              >  2:3825823/1‑150 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              >  2:2024903/1‑150 (MQ=255)
                                                                                                                                              tGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCg         <  2:2513166/150‑1 (MQ=255)
                                                                                                                                                gCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTc       <  1:347791/150‑1 (MQ=255)
                                                                                                                                                     ttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATg  <  1:2577226/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300184‑4300482

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: