Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,968,711 A→T 10.5% Q69L (CAG→CTG)  wecE → TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,968,7110AT10.5% 883.1 / 72.1 276Q69L (CAG→CTG) wecETDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase
Reads supporting (aligned to +/- strand):  ref base A (128/119);  new base T (10/19);  total (138/138)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 277 total aligned reads displayed.

GGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGAC  >  CP009273/3968562‑3968859
                                                                                                                                                     |                                                                                                                                                    
ggCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCa                                                                                                                                                      >  1:4211238/1‑150 (MQ=255)
    aTGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCt                                                                                                                                                  >  2:4206140/1‑150 (MQ=255)
      gggTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTgg                                                                                                                                                <  2:1793534/150‑1 (MQ=255)
      gggTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTgg                                                                                                                                                >  2:1396667/1‑150 (MQ=255)
       ggTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGc                                                                                                                                               <  2:2413019/150‑1 (MQ=255)
        gTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCg                                                                                                                                              >  1:1564387/1‑150 (MQ=255)
        gTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCg                                                                                                                                              <  1:880297/150‑1 (MQ=255)
           gCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATg                                                                                                                                           <  2:1089778/150‑1 (MQ=255)
            cGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGa                                                                                                                                          >  2:4615952/1‑150 (MQ=255)
                aaaCTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTg                                                                                                                                      <  2:1128317/150‑1 (MQ=255)
                     gtgtgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCAt                                                                                                                                 <  2:2999302/150‑1 (MQ=255)
                     gtgtgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCAt                                                                                                                                 <  1:2674256/150‑1 (MQ=255)
                     gtgtgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCAt                                                                                                                                 <  1:987235/150‑1 (MQ=255)
                      tgtgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGATGGAGCAACGTTTAGGCAGCGCCAAAGTGTAACTGACGCCGTCCTGCACAGCTTCGATGTAGATGGCGGCGCAGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATg                                                                                                                                >  2:3333252/1‑150 (MQ=255)
                        tgGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGcc                                                                                                                              <  2:3811679/150‑1 (MQ=255)
                         ggCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGAGCATGCCg                                                                                                                             >  1:4445281/1‑150 (MQ=255)
                            gATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGc                                                                                                                          <  1:578111/150‑1 (MQ=255)
                            gATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGc                                                                                                                          >  2:2725252/1‑150 (MQ=255)
                             aTGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCt                                                                                                                         <  2:2588587/150‑1 (MQ=255)
                              tGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTa                                                                                                                        <  2:1750326/150‑1 (MQ=255)
                               ggCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTac                                                                                                                       <  1:3153972/150‑1 (MQ=255)
                                gCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTaca                                                                                                                      >  2:1625377/1‑150 (MQ=255)
                                  ggTTTTACCCGTCGCTGCCAGAAGTGGCTGGAGCAACGTTTTGGCAGCGCAAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCTTGCCGAGCTACAcc                                                                                                                    <  2:750745/150‑1 (MQ=255)
                                     tttACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCttt                                                                                                                 <  1:3002590/150‑1 (MQ=255)
                                     tttACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCttt                                                                                                                 <  2:2898043/150‑1 (MQ=255)
                                      ttACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTg                                                                                                                <  1:3062398/150‑1 (MQ=255)
                                        aCCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtc                                                                                                              <  2:3749926/150‑1 (MQ=255)
                                         cccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtct                                                                                                             <  2:2308714/150‑1 (MQ=255)
                                         cccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtct                                                                                                             <  2:2768499/150‑1 (MQ=255)
                                          ccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtctc                                                                                                            >  1:1266484/1‑150 (MQ=255)
                                          ccGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGtctc                                                                                                            >  1:2954650/1‑150 (MQ=255)
                                           cGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTcc                                                                                                           >  1:176864/1‑150 (MQ=255)
                                           cGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTcc                                                                                                           >  1:2772442/1‑150 (MQ=255)
                                             tCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCAc                                                                                                         >  1:3691736/1‑150 (MQ=255)
                                             tCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCAc                                                                                                         >  1:1981822/1‑150 (MQ=255)
                                                 tGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGcc                                                                                                     >  1:2868626/1‑150 (MQ=255)
                                                 tGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCACTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGcc                                                                                                     >  2:910941/1‑150 (MQ=255)
                                                  gCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCa                                                                                                    >  2:3253297/1‑150 (MQ=255)
                                                   ccAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCaa                                                                                                   >  1:3314876/1‑150 (MQ=255)
                                                     agcagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATg                                                                                                 <  2:2993615/150‑1 (MQ=255)
                                                       cagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGcc                                                                                               <  2:2279011/150‑1 (MQ=255)
                                                       cagTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGcc                                                                                               >  2:2966544/1‑150 (MQ=255)
                                                        agTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCt                                                                                              <  1:450111/150‑1 (MQ=255)
                                                           ggCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtg                                                                                           <  1:171202/150‑1 (MQ=255)
                                                             cTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGACGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTtgtg                                                                                         >  2:1867516/1‑150 (MQ=255)
                                                                gAGCAACGTTGTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACAGCGTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGtgctg                                                                                      <  2:4472315/150‑1 (MQ=255)
                                                                 aGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGtgctgc                                                                                     >  1:1718436/1‑150 (MQ=255)
                                                                   cAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGt                                                                                   <  1:836508/150‑1 (MQ=255)
                                                                    aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTg                                                                                  >  1:1037418/1‑150 (MQ=255)
                                                                    aaCGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTAGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTg                                                                                  >  2:4516466/1‑150 (MQ=255)
                                                                      cGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGgc                                                                                >  1:4598617/1‑150 (MQ=255)
                                                                       gTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGgcg                                                                               >  2:93901/1‑150 (MQ=255)
                                                                         tttGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCa                                                                             >  2:3591476/1‑150 (MQ=255)
                                                                             gCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCaaaaa                                                                         <  1:3917082/150‑1 (MQ=255)
                                                                               aGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATc                                                                       <  2:1233439/150‑1 (MQ=255)
                                                                                   ccAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGttt                                                                   <  1:2661004/150‑1 (MQ=255)
                                                                                     aaaGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGttttt                                                                 >  2:831862/1‑150 (MQ=255)
                                                                                     aaaGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGttttt                                                                 <  1:4529427/150‑1 (MQ=255)
                                                                                     aaaGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGttttt                                                                 >  2:1280212/1‑150 (MQ=255)
                                                                                     aaaGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGttttt                                                                 <  1:972412/150‑1 (MQ=255)
                                                                                       aGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTtgt                                                               <  2:510523/150‑1 (MQ=255)
                                                                                           ttACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGAt                                                           <  1:830462/150‑1 (MQ=255)
                                                                                           ttACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGAt                                                           >  2:590786/1‑150 (MQ=255)
                                                                                              cTGACGCCGTCCTGCACAGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGtt                                                        <  1:163208/150‑1 (MQ=255)
                                                                                               tGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTc                                                       <  2:2925603/150‑1 (MQ=255)
                                                                                               tGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTc                                                       >  2:2658542/1‑150 (MQ=255)
                                                                                               tGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTc                                                       <  2:3684088/150‑1 (MQ=255)
                                                                                                   gCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGccc                                                   >  2:1486012/1‑150 (MQ=255)
                                                                                                    ccGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCg                                                  <  2:2060704/150‑1 (MQ=255)
                                                                                                      gTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGa                                                <  2:3102528/150‑1 (MQ=255)
                                                                                                      gTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGa                                                <  2:785152/150‑1 (MQ=255)
                                                                                                           gCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCa                                           >  1:1373800/1‑150 (MQ=255)
                                                                                                           gCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCa                                           >  1:2736976/1‑150 (MQ=255)
                                                                                                           gCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCa                                           >  1:3418418/1‑150 (MQ=255)
                                                                                                              ccGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGa                                        <  1:3547009/150‑1 (MQ=255)
                                                                                                                  ttCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACAt                                    >  1:2350543/1‑150 (MQ=255)
                                                                                                                   tCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATc                                   <  1:1460685/150‑1 (MQ=255)
                                                                                                                      cTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATcgac                                <  1:3873577/150‑1 (MQ=255)
                                                                                                                      cTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATcgac                                <  1:3366854/150‑1 (MQ=255)
                                                                                                                        ggAGATGGCGGCGCTGCTGCTCGATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATcgacga                              >  2:3932545/1‑150 (MQ=255)
                                                                                                                            aTGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACg                          <  1:4263087/150‑1 (MQ=255)
                                                                                                                              ggcggcGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCt                        <  1:681414/150‑1 (MQ=255)
                                                                                                                               gcggcgCTGCTGCTCGATATCCAGCCTGGCGATGTAGTGATCATGCCGATCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGTCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTg                       <  1:3201564/150‑1 (MQ=255)
                                                                                                                               gcggcgCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTg                       <  1:88449/150‑1 (MQ=255)
                                                                                                                               gcggcgCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTg                       <  1:1396667/150‑1 (MQ=255)
                                                                                                                                 ggcgCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGAt                     >  1:1830916/1‑150 (MQ=255)
                                                                                                                                  gcgcTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACTCCATGAACATCGACGAAACGCTGAtt                    <  1:206089/150‑1 (MQ=255)
                                                                                                                                   cgcTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTg                   <  2:1547452/150‑1 (MQ=255)
                                                                                                                                      tgctgctCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAg                >  2:3136773/1‑150 (MQ=255)
                                                                                                                                      tgctgctCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAg                >  2:1957739/1‑150 (MQ=255)
                                                                                                                                       gctgctCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAgc               <  1:2497225/150‑1 (MQ=255)
                                                                                                                                         tgctCGATATCCAGCCTGGCGATGAAGTGATCACGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAgcgg             >  1:2706846/1‑150 (MQ=255)
                                                                                                                                           ctCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAgcggcg           >  1:3601792/1‑150 (MQ=255)
                                                                                                                                            tCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGa          <  2:978239/150‑1 (MQ=255)
                                                                                                                                            tCGATATCCAGCCTGGCGATGAAGTGATCACGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGa          >  2:4138717/1‑150 (MQ=255)
                                                                                                                                              gATATCCTGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATc        <  2:1453078/150‑1 (MQ=255)
                                                                                                                                              gATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATc        <  1:4377317/150‑1 (MQ=255)
                                                                                                                                              gATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATc        >  2:2495996/1‑150 (MQ=255)
                                                                                                                                               atatCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCa       >  2:3471244/1‑150 (MQ=255)
                                                                                                                                                    cAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGAc  >  2:4167487/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCAATGGGTAGCGGCAAACTGTGTGGCGATGGCGGTTTTACCCGTCGCTGCCAGCAGTGGCTGGAGCAACGTTTTGGCAGCGCCAAAGTGTTACTGACGCCGTCCTGCACCGCTTCGCTGGAGATGGCGGCGCTGCTGCTCGATATCCAGCCTGGCGATGAAGTGATCATGCCGAGCTACACCTTTGTCTCCACCGCCAATGCCTTTGTGCTGCGTGGCGCAAAAATCGTTTTTGTGGATGTTCGCCCGGACACCATGAACATCGACGAAACGCTGATTGAAGCGGCGATCACCGAC  >  CP009273/3968562‑3968859

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: