Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,300,333 T→A 24.4% 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,3330TA24.4% 666.8 / 240.9 312E155D (GAA→GATalsAfused D‑allose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (111/125);  new base A (42/34);  total (153/159)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 316 total aligned reads displayed.

AACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300184‑4300482
                                                                                                                                                     |                                                                                                                                                     
  cTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTc                                                                                                                                                     <  1:4856524/150‑1 (MQ=255)
     aTGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAg                                                                                                                                                  <  2:2403529/150‑1 (MQ=255)
       gCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCa                                                                                                                                                >  2:3806678/1‑150 (MQ=255)
          agaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCt                                                                                                                                             >  1:4693126/1‑150 (MQ=255)
          agaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCt                                                                                                                                             >  1:2497320/1‑150 (MQ=255)
            atatAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGc                                                                                                                                           >  1:1902001/1‑150 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          >  2:1348444/1‑150 (MQ=255)
              ataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCtt                                                                                                                                         >  1:3146209/1‑150 (MQ=255)
              ataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCtt                                                                                                                                         <  2:1609269/150‑1 (MQ=255)
               taGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtg                                                                                                                                        <  2:2981422/150‑1 (MQ=255)
                aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgt                                                                                                                                       <  2:1738321/150‑1 (MQ=255)
                 gACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgtg                                                                                                                                      <  1:2851158/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   <  2:3298687/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   >  1:2736254/1‑150 (MQ=255)
                     aTGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTg                                                                                                                                  <  1:4350165/150‑1 (MQ=255)
                       ggCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgat                                                                                                                                <  1:1044446/150‑1 (MQ=255)
                         ccGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatag                                                                                                                              <  1:2930140/150‑1 (MQ=255)
                          cGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgataga                                                                                                                             >  2:3936543/1‑150 (MQ=255)
                            tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           >  2:3511292/1‑150 (MQ=255)
                              cccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATaaa                                                                                                                         >  1:3241169/1‑150 (MQ=255)
                               ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAt                                                                                                                        >  2:2433661/1‑150 (MQ=255)
                               ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAAt                                                                                                                        <  1:2019094/150‑1 (MQ=255)
                                ctctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAtt                                                                                                                       >  1:1554467/1‑150 (MQ=255)
                                 tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  1:343241/150‑1 (MQ=255)
                                  ctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCg                                                                                                                     >  2:4268383/1‑150 (MQ=255)
                                  ctTTACGCAACTGAGTCAGGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCg                                                                                                                     <  2:3768432/150‑1 (MQ=255)
                                     tACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCa                                                                                                                  >  2:2330041/1‑150 (MQ=255)
                                      aCGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCAc                                                                                                                 >  1:5088373/1‑150 (MQ=255)
                                        gCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt                                                                                                               <  2:4846895/150‑1 (MQ=255)
                                        gCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt                                                                                                               >  2:4012561/1‑150 (MQ=255)
                                         cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              >  2:3701613/1‑150 (MQ=255)
                                          aaCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc                                                                                                             >  2:3157561/1‑150 (MQ=255)
                                            cTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtctc                                                                                                           <  2:4562651/150‑1 (MQ=255)
                                              gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         <  2:5045218/150‑1 (MQ=255)
                                                ttCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCt                                                                                                       >  2:2634313/1‑150 (MQ=255)
                                                   aTGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGa                                                                                                    >  1:4809147/1‑150 (MQ=255)
                                                      aTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGGTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCa                                                                                                 >  2:4961776/1‑150 (MQ=255)
                                                       tCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCaa                                                                                                <  1:1237516/150‑1 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  1:534612/1‑150 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  1:3792540/1‑150 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  1:3854543/1‑150 (MQ=255)
                                                            aaCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTc                                                                                           <  1:977682/150‑1 (MQ=255)
                                                              cAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCaa                                                                                         <  1:3952469/150‑1 (MQ=255)
                                                                gATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGc                                                                                       >  2:4759225/1‑150 (MQ=255)
                                                                 aTAGTCCACCTCTTTAGTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGcc                                                                                      >  1:634482/1‑150 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     <  1:3562562/150‑1 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     <  1:71781/150‑1 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                <  2:4476974/150‑1 (MQ=255)
                                                                        aCCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGt                                                                               <  2:2241033/150‑1 (MQ=255)
                                                                         ccTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGta                                                                              <  2:2900142/150‑1 (MQ=255)
                                                                         ccTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGta                                                                              >  1:1418737/1‑150 (MQ=255)
                                                                            ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           >  1:2395945/1‑150 (MQ=255)
                                                                              ttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAc                                                                         >  1:2541644/1‑150 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     <  2:2945418/150‑1 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     <  2:4391045/150‑1 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     <  2:707604/150‑1 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     >  2:2026499/1‑150 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   >  2:1156634/1‑150 (MQ=255)
                                                                                     tgagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGATTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATgg                                                                  >  1:4432569/1‑150 (MQ=255)
                                                                                          gaggagGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGgcggcg                                                             >  1:4270879/1‑150 (MQ=255)
                                                                                           aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCGGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            <  2:1527720/150‑1 (MQ=255)
                                                                                           aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            >  2:1870822/1‑150 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           <  2:3452084/150‑1 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           <  1:4191200/150‑1 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           <  2:604934/150‑1 (MQ=255)
                                                                                              aggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGa                                                         >  1:167732/1‑150 (MQ=255)
                                                                                                 tGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGacac                                                      >  2:1935209/1‑150 (MQ=255)
                                                                                                   ggTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGc                                                    >  1:58353/1‑150 (MQ=255)
                                                                                                    gTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTAGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCa                                                   >  1:365382/1‑150 (MQ=255)
                                                                                                     ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAt                                                  <  1:4155778/150‑1 (MQ=255)
                                                                                                      tCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 >  2:804924/1‑150 (MQ=255)
                                                                                                      tCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 <  1:3950348/150‑1 (MQ=255)
                                                                                                       cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                >  1:3350954/1‑150 (MQ=255)
                                                                                                       cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                <  1:1630771/150‑1 (MQ=255)
                                                                                                          ccATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtctc                                             <  2:414878/150‑1 (MQ=255)
                                                                                                            atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  2:3566842/1‑150 (MQ=255)
                                                                                                             tgatgatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTCGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGcc                                          <  1:5127405/150‑1 (MQ=255)
                                                                                                             tgatgatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGcc                                          >  2:873333/1‑150 (MQ=255)
                                                                                                                 gatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTc                                      <  2:4732020/150‑1 (MQ=255)
                                                                                                                      cTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGataa                                 <  2:4214699/150‑1 (MQ=255)
                                                                                                                           gCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTg                            <  1:2444886/150‑1 (MQ=255)
                                                                                                                             aTCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGAc                          <  1:930013/150‑1 (MQ=255)
                                                                                                                               cGAGCATCAGCGTTTTGGCAATTTCGAGCATCTGCTTGTGGCTGATAGATAAAGTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGc                        <  2:3235316/150‑1 (MQ=255)
                                                                                                                                gAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGcc                       >  1:2424734/1‑150 (MQ=255)
                                                                                                                                  gCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCcac                     <  1:4639260/150‑1 (MQ=255)
                                                                                                                                   cATCAGCGTTTTGTCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCTTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCctca                    <  2:10123/150‑1 (MQ=255)
                                                                                                                                   cATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCcaca                    <  1:2221235/150‑1 (MQ=255)
                                                                                                                                      cAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAt                 >  1:5068898/1‑150 (MQ=255)
                                                                                                                                       aGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                <  2:2736254/150‑1 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              <  2:3484631/150‑1 (MQ=255)
                                                                                                                                           ttttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttttt            <  1:1080429/150‑1 (MQ=255)
                                                                                                                                            tttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttttttt           >  1:2326608/1‑150 (MQ=255)
                                                                                                                                            tttGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttttttt           <  1:493298/150‑1 (MQ=255)
                                                                                                                                             ttGGGAATTTCGAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCCTGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTc          <  1:1504023/150‑1 (MQ=255)
                                                                                                                                             ttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTc          <  1:3511115/150‑1 (MQ=255)
                                                                                                                                             ttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTc          >  2:703631/1‑150 (MQ=255)
                                                                                                                                              tGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCg         <  2:2383652/150‑1 (MQ=255)
                                                                                                                                                gCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTc       <  1:2111885/150‑1 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      <  2:2498933/150‑1 (MQ=255)
                                                                                                                                                    tttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGAt   <  1:4064899/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: