Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,833,263 C→A 28.3% V127V (GTC→GTA ascF → fused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,833,2630CA28.3% 652.1 / 209.5 288V127V (GTC→GTAascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
Reads supporting (aligned to +/- strand):  ref base C (113/93);  new base A (37/45);  total (150/138)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

CGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACC  >  CP009273/2833116‑2833407
                                                                                                                                                   |                                                                                                                                                
cgtcAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTaaa                                                                                                                                                >  2:3448177/1‑150 (MQ=255)
cgtcAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCaa                                                                                                                                                >  1:3172744/1‑150 (MQ=255)
  tcAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAAc                                                                                                                                              >  2:1238592/1‑150 (MQ=255)
   cAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAAct                                                                                                                                             <  2:2545933/150‑1 (MQ=255)
     gCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAActgc                                                                                                                                           >  1:3156811/1‑150 (MQ=255)
           tgcCGGGAGATATGCAGCCAGCACAGCCCGTGGGTAAACCCAAACACACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCaa                                                                                                                                     >  2:2007924/1‑150 (MQ=255)
             ccGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATg                                                                                                                                   >  1:1172489/1‑150 (MQ=255)
             ccGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATg                                                                                                                                   <  2:3545786/150‑1 (MQ=255)
             ccGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATg                                                                                                                                   <  2:2412869/150‑1 (MQ=255)
               gggAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGAt                                                                                                                                 >  1:1303359/1‑150 (MQ=255)
                  agaTATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCt                                                                                                                              <  1:3238310/150‑1 (MQ=255)
                  agaTATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCt                                                                                                                              <  2:3742236/150‑1 (MQ=255)
                  agaTATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCt                                                                                                                              <  2:51431/150‑1 (MQ=255)
                  agaTATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCt                                                                                                                              <  2:1664711/150‑1 (MQ=255)
                      atGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCgag                                                                                                                          <  1:1169839/150‑1 (MQ=255)
                       tGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCgaga                                                                                                                         >  1:2450848/1‑150 (MQ=255)
                       tGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCgaga                                                                                                                         >  1:1691081/1‑150 (MQ=255)
                        gCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTCTCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGAt                                                                                                                        >  2:4570011/1‑150 (MQ=255)
                          aGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGa                                                                                                                      <  2:1857732/150‑1 (MQ=255)
                          aGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGa                                                                                                                      >  1:3493480/1‑150 (MQ=255)
                           gCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAg                                                                                                                     <  2:1019764/150‑1 (MQ=255)
                              cGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAgcgg                                                                                                                  <  2:4363246/150‑1 (MQ=255)
                                   aGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACGGCTGGCAATGATCCTCGAGATGAGCGGCGTGc                                                                                                             >  2:1129397/1‑150 (MQ=255)
                                   aGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGc                                                                                                             >  1:3227218/1‑150 (MQ=255)
                                    gCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCt                                                                                                            <  1:2191228/150‑1 (MQ=255)
                                        gTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACa                                                                                                        <  1:3292284/150‑1 (MQ=255)
                                          gtgTAAACCCATACTCTCGCTTCGTCGCTTTGGTGCGGGGATCCTCGATGCGCTGTTCGGCTCCTTGTCACAGCTGATCCCGGAGATTATCGGCGGATCGATGGTCAAACTGCTGGCTATGATCCTCGAGATGAGCGGAGTTCTGACaaa                                                                                                      <  2:909056/148‑1 (MQ=255)
                                          gggTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACaaa                                                                                                      >  1:202169/1‑150 (MQ=255)
                                          gggTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACaaa                                                                                                      >  2:3530163/1‑150 (MQ=255)
                                              aaaCCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                                                                                  >  2:3256372/1‑150 (MQ=255)
                                              aaaCCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                                                                                  >  2:3574388/1‑150 (MQ=255)
                                               aaCCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGAt                                                                                                 >  2:2948914/1‑150 (MQ=255)
                                               aaCCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGAt                                                                                                 >  2:141240/1‑150 (MQ=255)
                                               aaCCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGAt                                                                                                 <  2:1214770/150‑1 (MQ=255)
                                                  ccAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                                                                              >  1:2095042/1‑150 (MQ=255)
                                                   cAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                                                                             >  1:4420285/1‑150 (MQ=255)
                                                      aCTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGAc                                                                                          <  1:4626920/150‑1 (MQ=255)
                                                        tCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCt                                                                                        >  1:603965/1‑150 (MQ=255)
                                                         cACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCtt                                                                                       <  1:371070/150‑1 (MQ=255)
                                                         cACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCtt                                                                                       >  1:1249139/1‑150 (MQ=255)
                                                            gCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAc                                                                                    <  2:4250127/150‑1 (MQ=255)
                                                             cTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                                                                                   >  2:4034578/1‑150 (MQ=255)
                                                             cTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                                                                                   >  2:1764711/1‑150 (MQ=255)
                                                             cTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                                                                                   >  2:1700068/1‑150 (MQ=255)
                                                              tACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCa                                                                                  >  1:3063167/1‑150 (MQ=255)
                                                              tACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCACCGACCTTAACCa                                                                                  >  1:1530172/1‑150 (MQ=255)
                                                                 gTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTc                                                                               >  1:2765815/1‑150 (MQ=255)
                                                                  tCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCt                                                                              <  1:3428425/150‑1 (MQ=255)
                                                                     cATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGaa                                                                           <  2:94200/150‑1 (MQ=255)
                                                                     cATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGaa                                                                           <  2:3316535/150‑1 (MQ=255)
                                                                      aTTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAt                                                                          <  2:4136087/150‑1 (MQ=255)
                                                                       ttGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtg                                                                         <  2:2498950/150‑1 (MQ=255)
                                                                        tGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgt                                                                        >  2:1922852/1‑150 (MQ=255)
                                                                        tGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgt                                                                        >  1:2945239/1‑150 (MQ=255)
                                                                         ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgtg                                                                       <  1:2203795/150‑1 (MQ=255)
                                                                          gTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGa                                                                      <  2:2877429/150‑1 (MQ=255)
                                                                          gTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGa                                                                      <  2:1501555/150‑1 (MQ=255)
                                                                            gCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGAtt                                                                    <  1:4488334/150‑1 (MQ=255)
                                                                                 gATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGa                                                               >  2:626749/1‑150 (MQ=255)
                                                                                   tCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACg                                                             >  2:3206829/1‑150 (MQ=255)
                                                                                    ccTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACgg                                                            <  2:1249238/150‑1 (MQ=255)
                                                                                      tCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGTGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTg                                                          <  2:840265/150‑1 (MQ=255)
                                                                                          tGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCttt                                                      >  1:674387/1‑150 (MQ=255)
                                                                                            cgcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttct                                                    >  2:588362/1‑150 (MQ=255)
                                                                                             gcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttctt                                                   >  2:1898520/1‑150 (MQ=255)
                                                                                             gcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttctt                                                   <  1:2619878/150‑1 (MQ=255)
                                                                                             gcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttctt                                                   >  1:1081045/1‑150 (MQ=255)
                                                                                                gATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttcttctt                                                >  1:2787198/1‑150 (MQ=255)
                                                                                                gATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttcttctt                                                >  2:1608942/1‑150 (MQ=255)
                                                                                                gATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTttcttctt                                                <  1:4105013/150‑1 (MQ=255)
                                                                                                   cGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCCTTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCt                                             >  1:4355866/1‑150 (MQ=255)
                                                                                                   cGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCt                                             <  2:2846624/150‑1 (MQ=255)
                                                                                                    ggCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTg                                            <  2:4476577/150‑1 (MQ=255)
                                                                                                      cACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGcc                                          <  2:595352/150‑1 (MQ=255)
                                                                                                      cACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGcc                                          >  1:2663580/1‑150 (MQ=255)
                                                                                                        ccATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGc                                        >  2:1186141/1‑150 (MQ=255)
                                                                                                        ccATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGc                                        >  1:4536116/1‑150 (MQ=255)
                                                                                                         cATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCt                                       <  1:3288191/150‑1 (MQ=255)
                                                                                                            gTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGAt                                    <  2:4317439/150‑1 (MQ=255)
                                                                                                                ccGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTc                                <  1:992874/150‑1 (MQ=255)
                                                                                                                ccGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTc                                >  1:189171/1‑150 (MQ=255)
                                                                                                                ccGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTc                                <  2:141076/150‑1 (MQ=255)
                                                                                                                 cGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg                               <  1:84488/150‑1 (MQ=255)
                                                                                                                 cGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcg                               <  1:3756152/150‑1 (MQ=255)
                                                                                                                   cTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTcgcg                             <  1:139812/150‑1 (MQ=255)
                                                                                                                     gATCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGc                           <  1:687584/150‑1 (MQ=255)
                                                                                                                      aTCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCa                          <  1:2248990/150‑1 (MQ=255)
                                                                                                                       tCCCGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCAt                         >  2:3624156/1‑150 (MQ=255)
                                                                                                                        cccGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATc                        >  1:1673977/1‑150 (MQ=255)
                                                                                                                         ccGGCGATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCt                       >  1:4265502/1‑150 (MQ=255)
                                                                                                                              gATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTgccgc                  >  2:3147052/1‑150 (MQ=255)
                                                                                                                              gATTATCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTgccgc                  >  2:2658416/1‑150 (MQ=255)
                                                                                                                                 tATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCAt               >  1:3099391/1‑150 (MQ=255)
                                                                                                                                   tCGGCGGATCGATGGTCAAACTGCTGGCAANGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa             >  1:4284830/1‑150 (MQ=255)
                                                                                                                                   tCGGCGGATCGATGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCa             >  2:2215875/1‑150 (MQ=255)
                                                                                                                                      gcggATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAt          <  2:274091/150‑1 (MQ=255)
                                                                                                                                       cggATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAAtt         >  1:3561699/1‑150 (MQ=255)
                                                                                                                                         gATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTa       >  1:4008955/1‑150 (MQ=255)
                                                                                                                                              aTGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc  <  1:3255101/150‑1 (MQ=255)
                                                                                                                                              aTGGTAAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAAcc  >  2:4404741/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGTCAGCCTGCTGCCGGGAGATATGCAGCCCGCACAGCCCGTGGGTAAACCCAAACTCACGCTACGTCGCATTGGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGACGGTGCTTTCTTCTTCCTGCCGCTGATGGTCGCGGCATCTGCCGCCATCAAATTTAAAACC  >  CP009273/2833116‑2833407

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: