Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,300,333 T→A 19.9% 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,3330TA19.9% 635.3 / 153.0 256E155D (GAA→GATalsAfused D‑allose transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base T (101/103);  new base A (21/29);  total (123/133)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 259 total aligned reads displayed.

CTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300186‑4300482
                                                                                                                                                   |                                                                                                                                                     
cTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATc                                                                                                                                                     >  1:1020812/1‑150 (MQ=255)
 ttATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCt                                                                                                                                                    <  2:2636141/150‑1 (MQ=255)
  tATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTa                                                                                                                                                   <  2:1836898/150‑1 (MQ=255)
   aTGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAg                                                                                                                                                  <  2:741473/150‑1 (MQ=255)
     gCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCa                                                                                                                                                >  1:2830325/1‑150 (MQ=255)
     gCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCa                                                                                                                                                >  1:1868451/1‑150 (MQ=255)
      cGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCAt                                                                                                                                               >  2:4235527/1‑150 (MQ=255)
         gaTATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTg                                                                                                                                            >  1:3199604/1‑150 (MQ=255)
           tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCt                                                                                                                                          <  1:775790/150‑1 (MQ=255)
           tataGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCt                                                                                                                                          <  1:1264926/150‑1 (MQ=255)
                   aTGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTg                                                                                                                                  >  2:431688/1‑150 (MQ=255)
                     ggCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgat                                                                                                                                <  1:2616313/150‑1 (MQ=255)
                      gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                                                                                                                               <  1:1162279/150‑1 (MQ=255)
                      gCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgata                                                                                                                               >  1:2414464/1‑150 (MQ=255)
                       ccGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatag                                                                                                                              <  2:3482569/150‑1 (MQ=255)
                       ccGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatag                                                                                                                              >  1:1834555/1‑150 (MQ=255)
                         gTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagat                                                                                                                            >  1:3011074/1‑150 (MQ=255)
                         gTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTgatagat                                                                                                                            >  1:3864323/1‑150 (MQ=255)
                          tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATTATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTTTGGCTgatagata                                                                                                                           >  1:2373957/1‑150 (MQ=255)
                          tACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTgatagata                                                                                                                           <  1:2920690/150‑1 (MQ=255)
                           aCCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATaa                                                                                                                          >  1:4085025/1‑150 (MQ=255)
                             ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAt                                                                                                                        <  2:2292164/150‑1 (MQ=255)
                             ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAt                                                                                                                        <  1:2890607/150‑1 (MQ=255)
                             ccTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAAt                                                                                                                        <  2:175030/150‑1 (MQ=255)
                               tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  1:95809/150‑1 (MQ=255)
                               tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      >  2:2935089/1‑150 (MQ=255)
                               tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  2:1913241/150‑1 (MQ=255)
                               tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  1:2312666/150‑1 (MQ=255)
                               tctTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTc                                                                                                                      <  1:341482/150‑1 (MQ=255)
                                  ttACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGcc                                                                                                                   >  2:848654/1‑150 (MQ=255)
                                  ttACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGcc                                                                                                                   <  1:3027626/150‑1 (MQ=255)
                                    aCGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCAc                                                                                                                 <  1:766128/150‑1 (MQ=255)
                                      gCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACtt                                                                                                               >  2:3272004/1‑150 (MQ=255)
                                       cAACTGATTCATGATCAGAAACAGTTAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              <  1:3565185/150‑1 (MQ=255)
                                       cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              <  1:1153685/150‑1 (MQ=255)
                                       cAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACttt                                                                                                              <  2:2134296/150‑1 (MQ=255)
                                        aaCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc                                                                                                             <  2:34136/150‑1 (MQ=255)
                                        aaCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtc                                                                                                             >  1:2529019/1‑150 (MQ=255)
                                         aCTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtct                                                                                                            <  2:1141059/150‑1 (MQ=255)
                                          cTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTtctc                                                                                                           >  1:289496/1‑150 (MQ=255)
                                           tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          >  2:3925283/1‑150 (MQ=255)
                                           tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          >  1:2225334/1‑150 (MQ=255)
                                           tGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCa                                                                                                          >  2:493543/1‑150 (MQ=255)
                                            gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         <  1:1032636/150‑1 (MQ=255)
                                            gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         <  1:4346139/150‑1 (MQ=255)
                                            gATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCAt                                                                                                         <  1:4131463/150‑1 (MQ=255)
                                             aTTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATc                                                                                                        <  2:2798445/150‑1 (MQ=255)
                                                cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAg                                                                                                     >  1:2784427/1‑150 (MQ=255)
                                                cATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAg                                                                                                     <  2:4325901/150‑1 (MQ=255)
                                                       aGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACt                                                                                              >  1:2082286/1‑150 (MQ=255)
                                                        gAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACtt                                                                                             >  1:3172479/1‑150 (MQ=255)
                                                          aaCAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTc                                                                                           <  1:2425918/150‑1 (MQ=255)
                                                            cAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCaa                                                                                         >  1:1725682/1‑150 (MQ=255)
                                                             aGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAg                                                                                        <  2:2913959/150‑1 (MQ=255)
                                                                tAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGccc                                                                                     >  2:1804642/1‑150 (MQ=255)
                                                                 aGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCa                                                                                    >  1:1989360/1‑150 (MQ=255)
                                                                  gTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAc                                                                                   <  2:2201142/150‑1 (MQ=255)
                                                                   tCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCAcg                                                                                  <  1:338444/150‑1 (MQ=255)
                                                                      aCCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGt                                                                               >  2:194923/1‑150 (MQ=255)
                                                                          ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           >  1:1620679/1‑150 (MQ=255)
                                                                          ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           <  1:1411024/150‑1 (MQ=255)
                                                                          ctTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGtaata                                                                           <  1:3975681/150‑1 (MQ=255)
                                                                            ttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAc                                                                         >  2:3170959/1‑150 (MQ=255)
                                                                            ttATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAc                                                                         >  1:3917190/1‑150 (MQ=255)
                                                                                tGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAAcatca                                                                     >  1:99243/1‑150 (MQ=255)
                                                                                   tgagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATgg                                                                  >  2:567176/1‑150 (MQ=255)
                                                                                   tgagtgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATATAGATCAACTTTCAAGCCCACGCGTAATAACATCATgg                                                                  <  1:901455/150‑1 (MQ=255)
                                                                                       tgagGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATggcggc                                                              >  2:928878/1‑150 (MQ=255)
                                                                                        gaggagGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGgcggcg                                                             >  2:2778151/1‑150 (MQ=255)
                                                                                         aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            >  1:745817/1‑150 (MQ=255)
                                                                                         aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            >  1:1460122/1‑150 (MQ=255)
                                                                                         aggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGgcgc                                                            >  2:2559079/1‑150 (MQ=255)
                                                                                          ggaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGcgcg                                                           <  2:3040596/150‑1 (MQ=255)
                                                                                           gaggTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCgg                                                          >  1:2609832/1‑150 (MQ=255)
                                                                                              gTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGaca                                                       >  2:4154959/1‑150 (MQ=255)
                                                                                              gTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGaca                                                       >  2:3800561/1‑150 (MQ=255)
                                                                                                    tCATCCCTGATGATGACTTTGGCCTCGAGCATCAGCGTTTTGGCAATGTCTAGCATCTGCTTGTGGCTGATAGATAAAGTCGCCACTTTCTCATCTAGATCAACTTTCCAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 <  1:848659/150‑1 (MQ=255)
                                                                                                    tCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAGGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 >  2:1243902/1‑150 (MQ=255)
                                                                                                    tCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAtt                                                 <  2:1105338/150‑1 (MQ=255)
                                                                                                     cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATTGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                <  2:3830049/150‑1 (MQ=255)
                                                                                                     cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                <  2:646585/150‑1 (MQ=255)
                                                                                                     cATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGACAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCAttt                                                >  1:3825264/1‑150 (MQ=255)
                                                                                                        ccATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTtctc                                             >  1:1911226/1‑150 (MQ=255)
                                                                                                         cATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCg                                            >  1:1481720/1‑150 (MQ=255)
                                                                                                          atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  2:2022181/1‑150 (MQ=255)
                                                                                                          atgatgatgACTTTGGCATCGAGCATCAGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGc                                           >  2:1557570/1‑150 (MQ=255)
                                                                                                              tgatgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGt                                       >  1:2870291/1‑150 (MQ=255)
                                                                                                                atgaCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCg                                     <  2:2441581/150‑1 (MQ=255)
                                                                                                                   aCTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGata                                  >  2:1468739/1‑150 (MQ=255)
                                                                                                                     tttGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTCCTCGCCAGTCGataat                                <  1:3700294/150‑1 (MQ=255)
                                                                                                                                 cATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCcaca                    >  2:1212259/1‑150 (MQ=255)
                                                                                                                                 cATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCcaca                    >  2:4049040/1‑150 (MQ=255)
                                                                                                                                     aGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                <  1:2545802/150‑1 (MQ=255)
                                                                                                                                     aGCGTTTTGGCAATTTCTAGCATCTGATTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                <  2:3017147/150‑1 (MQ=255)
                                                                                                                                     aGCGTTTTGGCAATATCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGAtt                >  2:3556807/1‑150 (MQ=255)
                                                                                                                                           ttGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTc          <  2:4181426/150‑1 (MQ=255)
                                                                                                                                             ggCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGt        >  1:62584/1‑150 (MQ=255)
                                                                                                                                               cAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCa      <  1:2588757/150‑1 (MQ=255)
                                                                                                                                                aaTTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAg     >  1:4331774/1‑150 (MQ=255)
                                                                                                                                                  tttCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGAt   <  1:1183706/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTTATGCGAGATATAGACGATGGCCGTACCCTCTTTACGCAACTGATTCATGATCAGAAACAGATAGTCCACCTCTTTATTGGTGAGTGAGGAGGTGGGTTCATCCATGATGATGACTTTGGCATCGAGCATCAGCGTTTTGGCAATTTCTAGCATCTGCTTGTGGCTGATAGATAAATTCGCCACTTTCTCATCTAGATCAACTTTCAAGCCCACGCGTAATAACATCATGGCGGCGCGGACACGCATTTCTCGCCAGTCGATAATATTGACGCCACAGATTTTTTTCGTCAGATG  >  CP009273/4300186‑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: