Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,300,333 T→A 14.0% 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,3330TA14.0% 881.6 / 118.0 310E155D (GAA→GATalsAfused D‑allose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (139/127);  new base A (19/24);  total (158/152)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 318 total aligned reads displayed.

AACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300184‑4300482
                                                                                                                                                     |                                                                                                                                                     
aaCTTATGCGAGATATAGACGATGGCCGGACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAAtt                                                                                                                                                       <  2:624734/150‑1 (MQ=255)
   ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  1:3455715/150‑1 (MQ=255)
   ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  1:42343/150‑1 (MQ=255)
    tATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTa                                                                                                                                                   <  2:3802194/150‑1 (MQ=255)
    tATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTa                                                                                                                                                   >  1:101974/1‑150 (MQ=255)
      tGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGc                                                                                                                                                 >  2:1118423/1‑150 (MQ=255)
            atatAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGc                                                                                                                                           >  1:2438655/1‑150 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          >  1:2586771/1‑150 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          <  1:2501591/150‑1 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCt                                                                                                                                          <  1:472622/150‑1 (MQ=255)
                aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTtgt                                                                                                                                       <  2:1833764/150‑1 (MQ=255)
                 gACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtgtg                                                                                                                                      >  1:4196030/1‑150 (MQ=255)
                 gACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTtgtg                                                                                                                                      <  2:3075474/150‑1 (MQ=255)
                   cGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGc                                                                                                                                    <  2:801780/150‑1 (MQ=255)
                   cGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGc                                                                                                                                    >  1:3285036/1‑150 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   >  1:4181308/1‑150 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   <  2:3386055/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   >  2:1049320/1‑150 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   <  1:17320/150‑1 (MQ=255)
                    gATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCt                                                                                                                                   >  1:1044547/1‑150 (MQ=255)
                     agggCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTg                                                                                                                                  <  2:2499679/148‑1 (MQ=255)
                     aTGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTg                                                                                                                                  <  1:3879747/150‑1 (MQ=255)
                      tGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTga                                                                                                                                 <  2:837699/150‑1 (MQ=255)
                      tGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGACCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCATCGTTTTGGCAATTTCTATCATCTGCTTGTGGCTga                                                                                                                                 >  1:4525346/1‑150 (MQ=255)
                        gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                                                                                                                               <  1:2390342/150‑1 (MQ=255)
                        gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                                                                                                                               <  1:2272497/150‑1 (MQ=255)
                        gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                                                                                                                               <  1:88427/150‑1 (MQ=255)
                             aCCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaa                                                                                                                          >  2:3972097/1‑150 (MQ=255)
                             aCCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaa                                                                                                                          >  2:3971855/1‑150 (MQ=255)
                                ctctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAtt                                                                                                                       <  2:1285591/150‑1 (MQ=255)
                                 tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      >  1:761024/1‑150 (MQ=255)
                                   tttACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGc                                                                                                                    <  2:1863813/150‑1 (MQ=255)
                                         cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGCAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              >  1:4450564/1‑150 (MQ=255)
                                         cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTCAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              >  1:4475882/1‑150 (MQ=255)
                                          aaCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc                                                                                                             <  2:1166348/150‑1 (MQ=255)
                                           aCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtct                                                                                                            <  2:3664527/150‑1 (MQ=255)
                                             tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          <  1:2132164/150‑1 (MQ=255)
                                               aTTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc                                                                                                        >  2:2515908/1‑150 (MQ=255)
                                                 tCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTa                                                                                                      <  1:1992458/150‑1 (MQ=255)
                                                  cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAg                                                                                                     >  2:2038671/1‑150 (MQ=255)
                                                   aTGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGa                                                                                                    <  2:3520581/150‑1 (MQ=255)
                                                     gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATc                                                                                                  >  2:67219/1‑150 (MQ=255)
                                                       tCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCaa                                                                                                >  1:3194243/1‑150 (MQ=255)
                                                       tCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCaa                                                                                                >  1:2009966/1‑150 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  2:4367583/1‑150 (MQ=255)
                                                        cAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAAc                                                                                               >  1:814451/1‑150 (MQ=255)
                                                         aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACt                                                                                              <  2:2134458/150‑1 (MQ=255)
                                                         aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGAGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCTCTTTCTCATCTAGATCAACt                                                                                              <  2:4383525/150‑1 (MQ=255)
                                                             aCAGATAGTCCACCCCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCa                                                                                          <  1:4412905/150‑1 (MQ=255)
                                                               aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAg                                                                                        >  2:3840226/1‑150 (MQ=255)
                                                                gATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGc                                                                                       >  2:417412/1‑150 (MQ=255)
                                                                    gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCGACTTTCAAGCCCAc                                                                                   <  2:852730/150‑1 (MQ=255)
                                                                    gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAc                                                                                   <  2:3132395/150‑1 (MQ=255)
                                                                      ccACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgc                                                                                 >  2:2176968/1‑150 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                <  1:3148129/150‑1 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                <  1:3874777/150‑1 (MQ=255)
                                                                        aCCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGt                                                                               >  2:1950195/1‑150 (MQ=255)
                                                                            ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           <  1:4006797/150‑1 (MQ=255)
                                                                             tttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaataa                                                                          >  1:3844860/1‑150 (MQ=255)
                                                                              ttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAc                                                                         <  1:1887119/150‑1 (MQ=255)
                                                                                aTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcat                                                                       >  1:1191480/1‑150 (MQ=255)
                                                                                aTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcat                                                                       >  2:3619877/1‑150 (MQ=255)
                                                                                aTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcat                                                                       >  2:1638805/1‑150 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     >  2:221822/1‑150 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   <  2:1609614/150‑1 (MQ=255)
                                                                                      gagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggc                                                                 <  2:4037957/150‑1 (MQ=255)
                                                                                      gagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggc                                                                 >  1:971955/1‑150 (MQ=255)
                                                                                      gagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggc                                                                 <  2:2976215/150‑1 (MQ=255)
                                                                                       agtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcg                                                                <  2:2278823/150‑1 (MQ=255)
                                                                                        gtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcgg                                                               <  2:3100149/150‑1 (MQ=255)
                                                                                         tgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcggc                                                              >  2:2688286/1‑150 (MQ=255)
                                                                                           aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            <  1:3677933/150‑1 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           >  1:4337697/1‑150 (MQ=255)
                                                                                            ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           >  2:1412160/1‑150 (MQ=255)
                                                                                             gaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCgg                                                          >  1:1267189/1‑150 (MQ=255)
                                                                                             gaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCgg                                                          >  1:2946222/1‑150 (MQ=255)
                                                                                               ggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGac                                                        <  1:3508697/150‑1 (MQ=255)
                                                                                                    gTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCa                                                   <  1:853398/150‑1 (MQ=255)
                                                                                                      tCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 >  2:2717010/1‑150 (MQ=255)
                                                                                                           cATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCg                                            >  2:1376271/1‑150 (MQ=255)
                                                                                                            atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  1:1618376/1‑150 (MQ=255)
                                                                                                            atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  1:740569/1‑150 (MQ=255)
                                                                                                            atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  1:98053/1‑150 (MQ=255)
                                                                                                                   tgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGa                                    <  2:732317/150‑1 (MQ=255)
                                                                                                                          ggCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATAtt                             >  1:1522950/1‑150 (MQ=255)
                                                                                                                           gCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTg                            <  2:1264144/150‑1 (MQ=255)
                                                                                                                           gCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTg                            <  2:1263851/150‑1 (MQ=255)
                                                                                                                             aTCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGAc                          <  1:3054605/150‑1 (MQ=255)
                                                                                                                                       aGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                >  2:2739997/1‑150 (MQ=255)
                                                                                                                                       aGCGTTTTGGCAATTTCTAGCATCTGATTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATAAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCACGCCAGTCGATAATATTGACGCCAAAGAtt                >  2:2739739/1‑150 (MQ=255)
                                                                                                                                        gCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttt               <  1:67405/150‑1 (MQ=255)
                                                                                                                                          gTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttttt             >  2:503456/1‑150 (MQ=255)
                                                                                                                                            tttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttttttt           <  1:3518758/150‑1 (MQ=255)
                                                                                                                                            tttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAttttttt           <  2:4488609/150‑1 (MQ=255)
                                                                                                                                             ttGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTc          >  2:1768540/1‑150 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCTAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      <  1:2684698/150‑1 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      <  2:4181391/150‑1 (MQ=255)
                                                                                                                                                  aaTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAg     >  2:4298647/1‑150 (MQ=255)
                                                                                                                                                  aaTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAg     >  2:3508788/1‑150 (MQ=255)
                                                                                                                                                   aTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGa    >  2:1471822/1‑150 (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: