Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,300,333 T→A 7.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,3330TA 7.1% 860.4 / 43.1 252E155D (GAA→GATalsAfused D‑allose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (109/124);  new base A (10/8);  total (120/132)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 256 total aligned reads displayed.

AACTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300184‑4300482
                                                                                                                                                     |                                                                                                                                                     
   ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  1:1631049/150‑1 (MQ=255)
   ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  2:2966545/150‑1 (MQ=255)
   tgaggCGAGATATAGACGAGGGCCGGACCCTCTGTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  1:3908902/146‑1 (MQ=255)
     aTGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAg                                                                                                                                                  >  1:343491/1‑150 (MQ=255)
       gCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCa                                                                                                                                                >  1:4500638/1‑150 (MQ=255)
         gagaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATc                                                                                                                                              <  1:959798/150‑1 (MQ=255)
         gagaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATc                                                                                                                                              >  2:2856947/1‑150 (MQ=255)
          agaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCt                                                                                                                                             >  1:1022081/1‑150 (MQ=255)
            atatAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGAATCGAGCATCATCGTTTTGGCAATTTCTAGCaaatgc                                                                                                                                           >  2:123940/1‑145 (MQ=255)
             tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          >  1:2183362/1‑150 (MQ=255)
               taGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTtg                                                                                                                                        <  2:114772/150‑1 (MQ=255)
                aGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTtgt                                                                                                                                       <  1:4220131/150‑1 (MQ=255)
                       ggCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                                                                                                                                >  1:4451397/1‑150 (MQ=255)
                          cGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgataga                                                                                                                             <  1:2415620/150‑1 (MQ=255)
                          cGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgataga                                                                                                                             <  1:788347/150‑1 (MQ=255)
                          cGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgataga                                                                                                                             >  1:3998661/1‑150 (MQ=255)
                           gTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagat                                                                                                                            <  2:1517651/150‑1 (MQ=255)
                              cccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaaa                                                                                                                         >  2:2094010/1‑150 (MQ=255)
                              cccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaaa                                                                                                                         >  2:3804139/1‑150 (MQ=255)
                                 tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      >  1:3451339/1‑150 (MQ=255)
                                 tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  1:2683575/150‑1 (MQ=255)
                                       cGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACt                                                                                                                <  2:2922699/150‑1 (MQ=255)
                                        gCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt                                                                                                               >  1:3597663/1‑150 (MQ=255)
                                          aaCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc                                                                                                             <  1:839629/150‑1 (MQ=255)
                                           aCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtct                                                                                                            <  2:3198572/150‑1 (MQ=255)
                                             tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          <  2:3277771/150‑1 (MQ=255)
                                              gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         >  1:1198335/1‑150 (MQ=255)
                                              gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         >  1:1198194/1‑150 (MQ=255)
                                               aTTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc                                                                                                        >  2:918985/1‑150 (MQ=255)
                                               aTTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc                                                                                                        <  1:2435112/150‑1 (MQ=255)
                                                 tCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTa                                                                                                      >  2:2398071/1‑150 (MQ=255)
                                                  cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCAGATAGATAAATTCGCCACTTTCTCATCTAg                                                                                                     <  1:1739346/150‑1 (MQ=255)
                                                   aTGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGa                                                                                                    >  1:1928249/1‑150 (MQ=255)
                                                     gATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATc                                                                                                  >  1:3219609/1‑150 (MQ=255)
                                                      aTCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATGTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTGTCTCATCTAGATCa                                                                                                 >  2:4501712/1‑150 (MQ=255)
                                                         aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACt                                                                                              <  2:3150903/150‑1 (MQ=255)
                                                          gAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACtt                                                                                             >  1:2343254/1‑150 (MQ=255)
                                                             aCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCa                                                                                          <  2:3200801/150‑1 (MQ=255)
                                                              cAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCaa                                                                                         <  2:1072594/150‑1 (MQ=255)
                                                               aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCatg                                                                                        >  1:3385582/1‑148 (MQ=255)
                                                               aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAg                                                                                        >  1:2624647/1‑150 (MQ=255)
                                                               aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAg                                                                                        >  2:3687897/1‑150 (MQ=255)
                                                                gATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGc                                                                                       >  2:1599594/1‑150 (MQ=255)
                                                                 aTAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGcc                                                                                      >  1:3405219/1‑150 (MQ=255)
                                                                  tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGAATCTTCTTGTGGCTGATAGATAAATTCGCCACTTTCTAATCTAGATCAACTTACAAcccc                                                                                     >  1:1621128/1‑150 (MQ=255)
                                                                      ccACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgc                                                                                 >  2:1807745/1‑150 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                <  1:2904082/150‑1 (MQ=255)
                                                                       cACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcgcg                                                                                <  1:641825/150‑1 (MQ=255)
                                                                        aCCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGAAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGt                                                                               >  1:732458/1‑150 (MQ=255)
                                                                          ctctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaa                                                                             >  1:4457903/1‑150 (MQ=255)
                                                                          ctctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaa                                                                             <  2:2470542/150‑1 (MQ=255)
                                                                          ctctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaa                                                                             >  2:4374213/1‑150 (MQ=255)
                                                                            ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           <  2:449911/150‑1 (MQ=255)
                                                                               aaTTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTTGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTACAAGACCACGCGTAATAAca                                                                        >  2:255290/2‑150 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     <  1:1331616/150‑1 (MQ=255)
                                                                                  tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     <  1:2122285/150‑1 (MQ=255)
                                                                                    gtgagtgaGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATg                                                                   <  2:4181248/150‑1 (MQ=255)
                                                                                     tgagtgagGAGGTGGTTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATgg                                                                  <  2:697020/150‑1 (MQ=255)
                                                                                       agtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcg                                                                >  1:83601/1‑150 (MQ=255)
                                                                                                     ttCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAt                                                  <  1:3422828/150‑1 (MQ=255)
                                                                                                       cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                <  2:3264785/150‑1 (MQ=255)
                                                                                                        aTCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtc                                               <  2:1045950/150‑1 (MQ=255)
                                                                                                        aTCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtc                                               >  2:4447518/1‑150 (MQ=255)
                                                                                                        aTCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtc                                               >  1:3336453/1‑150 (MQ=255)
                                                                                                          ccATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtctc                                             >  2:3714373/1‑150 (MQ=255)
                                                                                                          ccATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtctc                                             >  1:907445/1‑150 (MQ=255)
                                                                                                             tgatgatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGcc                                          >  1:2508067/1‑150 (MQ=255)
                                                                                                             tgatgatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGcc                                          >  2:3329997/1‑150 (MQ=255)
                                                                                                                 gatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTc                                      <  1:1673884/150‑1 (MQ=255)
                                                                                                                 gatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTc                                      <  2:3524960/150‑1 (MQ=255)
                                                                                                                    gaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGat                                   <  2:2229390/150‑1 (MQ=255)
                                                                                                                     aCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGata                                  <  1:2551342/150‑1 (MQ=255)
                                                                                                                         tGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAatat                              >  2:1922927/1‑150 (MQ=255)
                                                                                                                         tGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAatat                              >  1:4390761/1‑150 (MQ=255)
                                                                                                                          ggCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATAtt                             <  2:2307070/150‑1 (MQ=255)
                                                                                                                          ggCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATAtt                             <  1:3711222/150‑1 (MQ=255)
                                                                                                                             aTCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGAc                          <  2:1489083/150‑1 (MQ=255)
                                                                                                                             aTCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGAc                          >  1:4191321/1‑150 (MQ=255)
                                                                                                                             aTCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGAc                          >  1:223194/1‑150 (MQ=255)
                                                                                                                                gAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAGTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGcc                       <  2:1686693/150‑1 (MQ=255)
                                                                                                                                gAGCATCAGCGTTTTGGAAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGcc                       <  1:2388334/150‑1 (MQ=255)
                                                                                                                                 aGCATCAGCGTTTTTGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCca                      <  2:237046/150‑1 (MQ=255)
                                                                                                                                 aGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCca                      >  1:3879374/1‑150 (MQ=255)
                                                                                                                                    aTCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAg                   >  2:2435582/1‑150 (MQ=255)
                                                                                                                                      cAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCCTCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAt                 <  2:3992167/150‑1 (MQ=255)
                                                                                                                                      cAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAt                 >  1:2027591/1‑150 (MQ=255)
                                                                                                                                       aGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                >  2:1748428/1‑150 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              >  1:4243658/1‑150 (MQ=255)
                                                                                                                                         cGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttt              >  2:1938450/1‑150 (MQ=255)
                                                                                                                                           ttttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtttttt            <  1:533119/150‑1 (MQ=255)
                                                                                                                                               ggCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGt        <  1:3290289/150‑1 (MQ=255)
                                                                                                                                               ggCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGt        <  1:2175659/150‑1 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      >  2:408384/1‑150 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      >  1:4475641/1‑150 (MQ=255)
                                                                                                                                                 cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      >  1:1785850/1‑150 (MQ=255)
                                                                                                                                                  aaTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAg     <  1:3669406/150‑1 (MQ=255)
                                                                                                                                                  aaTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAg     >  2:642655/1‑150 (MQ=255)
                                                                                                                                                    tttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGAt   >  2:3857203/1‑150 (MQ=255)
                                                                                                                                                    tttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGAt   <  1:1066681/150‑1 (MQ=255)
                                                                                                                                                     ttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCCTCTAGATCAACTTTCACGCCCACGCGTAATAACATCAGGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATg  <  2:189487/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: