Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,533,023 A→G K446R (AAG→AGG)  RPPX_06795 → GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,533,0230AG99.6% 865.2 / ‑4.7 248K446R (AAG→AGG) RPPX_06795GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (108/138);  total (109/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 248 total aligned reads displayed.

GCGACGAATTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAAGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACGGAATA  >  CP009974/1532876‑1533170
                                                                                                                                                   |                                                                                                                                                   
  gacgaATTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTTATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGCCCAGGCTgg                                                                                                                                                >  2:3248790/1‑151 (MQ=255)
     gaATTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCaa                                                                                                                                             >  1:4688259/1‑151 (MQ=255)
      aaTTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAg                                                                                                                                            >  1:4120583/1‑151 (MQ=255)
      aaTTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAg                                                                                                                                            >  1:3053154/1‑151 (MQ=255)
           gCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCctgc                                                                                                                                       >  1:3634982/1‑151 (MQ=255)
             cAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCctgctg                                                                                                                                     <  1:3134606/151‑1 (MQ=255)
               aCCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGa                                                                                                                                   <  2:4995780/151‑1 (MQ=255)
                 cTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGAtt                                                                                                                                 >  1:5001323/1‑151 (MQ=255)
                 cTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGAtt                                                                                                                                 >  2:254640/1‑151 (MQ=255)
                   ggCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCg                                                                                                                               <  1:4773892/151‑1 (MQ=255)
                    gCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGc                                                                                                                              <  2:3419773/151‑1 (MQ=255)
                       cGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCgg                                                                                                                           <  2:3485702/151‑1 (MQ=255)
                           aCAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAAcc                                                                                                                       <  2:1447492/151‑1 (MQ=255)
                            cagcGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACcg                                                                                                                      <  1:1787486/148‑1 (MQ=255)
                            cAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACcg                                                                                                                      >  1:3550838/1‑151 (MQ=255)
                             aaCGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACcgc                                                                                                                     <  2:415478/151‑1 (MQ=255)
                             aaCGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACcgc                                                                                                                     <  2:5213020/151‑1 (MQ=255)
                                 gCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGgtt                                                                                                                 <  2:1691876/151‑1 (MQ=255)
                                   gTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTTATGGCACCCAGCCGCCAGGCACTCGGCTCTGCCCCGACCCGGCTGTCAAGCCTCCTGGATTCCCCGGAACCCCggccgc                                                                                                               >  1:1050307/1‑147 (MQ=255)
                                   gTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGgttgt                                                                                                               <  1:2427514/151‑1 (MQ=255)
                                    tAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGgttgtt                                                                                                              >  1:5143442/1‑151 (MQ=255)
                                     aaCAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCACGCCTGCTGGATTCGCCGGAACCGCGgttgtcc                                                                                                             >  2:2688640/1‑149 (MQ=255)
                                      aCAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCt                                                                                                            <  1:4969052/151‑1 (MQ=255)
                                        aGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTAc                                                                                                          <  1:2030991/151‑1 (MQ=255)
                                         gTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGCTTCGCCGGAACCGCGGTTGTTCTACt                                                                                                         >  2:3710615/1‑151 (MQ=255)
                                             tCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGCACCGCGGTTGTTCTACTGcaa                                                                                                     >  2:3039374/1‑151 (MQ=255)
                                             tCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGaaa                                                                                                     <  1:4853030/151‑1 (MQ=255)
                                               ggCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGc                                                                                                   <  1:1149519/150‑1 (MQ=255)
                                               aGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGc                                                                                                   >  1:893154/1‑151 (MQ=255)
                                               aGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGc                                                                                                   <  1:3230135/151‑1 (MQ=255)
                                                    cAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGc                                                                                              >  1:935378/1‑151 (MQ=255)
                                                    cAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGc                                                                                              <  1:1498857/151‑1 (MQ=255)
                                                    cAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATCCGCCGGAACCGCGGTTGTTCTACTGAAAGCGACGc                                                                                              >  2:4204448/1‑151 (MQ=255)
                                                     aGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGcc                                                                                             >  1:1657773/1‑151 (MQ=255)
                                                     aGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGcc                                                                                             <  1:1327381/151‑1 (MQ=255)
                                                        cAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCAAATCTCTCAGATGGCACCCAGCCGACAGGCCCTGCGCTTTGCAATGACCAGGCTGGCAAGCCCGCTGGATTCGCCGGAACCGCGGTTCTTCTACAGAAGGCGCTGcacga                                                                                          >  2:3020299/1‑147 (MQ=255)
                                                         aGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGAt                                                                                         >  2:1663463/1‑151 (MQ=255)
                                                           ttCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATAc                                                                                       >  2:4278759/1‑151 (MQ=255)
                                                              cTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATAcgcg                                                                                    <  2:3806281/151‑1 (MQ=255)
                                                               tGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCgg                                                                                   <  2:1557711/151‑1 (MQ=255)
                                                                gCCAGGACGCAGCGACGCGAGGAGCGGGGGGAGCAGCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCGGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTGGTTCTGCTGAAAGCGATGCCAGATACGCGGt                                                                                  <  2:4019405/151‑1 (MQ=255)
                                                                 cgaGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGtg                                                                                 <  2:574825/149‑1 (MQ=255)
                                                                 ccATAACGCATCCGCGCGCCCGGCGGCGTGCGCAGCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGtg                                                                                 <  2:4378798/151‑1 (MQ=255)
                                                                  cAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAAGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGtgc                                                                                <  2:809519/151‑1 (MQ=255)
                                                                    gAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGtgctg                                                                              <  2:425210/151‑1 (MQ=255)
                                                                     aaCGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGtgctgc                                                                             >  1:2122871/1‑151 (MQ=255)
                                                                      aCGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGcc                                                                            <  2:2705614/151‑1 (MQ=255)
                                                                        gCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCa                                                                          <  1:962883/151‑1 (MQ=255)
                                                                         cggccacgccacCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCAt                                                                         <  2:454401/149‑1 (MQ=255)
                                                                          aGCCACGCCACCCGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTCGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATa                                                                        <  1:2966905/151‑1 (MQ=255)
                                                                          aGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATa                                                                        >  2:3595703/1‑151 (MQ=255)
                                                                          aGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATa                                                                        >  1:2468107/1‑151 (MQ=255)
                                                                           gccacgccacCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATaa                                                                       <  1:759349/151‑1 (MQ=255)
                                                                           gccacgccacCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATaa                                                                       <  2:2295470/151‑1 (MQ=255)
                                                                            ccacgccacCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATaaa                                                                      >  2:946576/1‑151 (MQ=255)
                                                                                   accaGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAAccc                                                               <  2:4585898/151‑1 (MQ=255)
                                                                                   accaGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAAccc                                                               >  1:2680331/1‑151 (MQ=255)
                                                                                      aGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAgg                                                            <  1:4221787/151‑1 (MQ=255)
                                                                                       gCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAggg                                                           <  2:1936075/151‑1 (MQ=255)
                                                                                         ggTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGCGCa                                                         >  2:3526676/1‑151 (MQ=255)
                                                                                         ggTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCACACAAAACCCCAgcgcc                                                         >  2:1803353/1‑150 (MQ=255)
                                                                                         ggTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCACGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCa                                                         <  1:130548/151‑1 (MQ=255)
                                                                                          gtgtGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTTCTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCaa                                                        >  1:4318491/1‑151 (MQ=255)
                                                                                            gtgCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACt                                                      <  1:2888483/151‑1 (MQ=255)
                                                                                             tgCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCAGAAAAAAACCCAGGGCACCTg                                                     >  2:5096070/1‑151 (MQ=255)
                                                                                             tgCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGAAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTg                                                     >  1:706003/1‑151 (MQ=255)
                                                                                                gcATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGccc                                                  <  1:2233872/151‑1 (MQ=255)
                                                                                                  aTCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTg                                                >  2:1970089/1‑151 (MQ=255)
                                                                                                  aTCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTg                                                <  1:1448289/151‑1 (MQ=255)
                                                                                                       tctGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCCAGCCTGCTGGATTCGCCGGCACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATCAAAAAACCCAGCGCAACTGCCCTggcctt                                           >  2:4416354/1‑147 (MQ=255)
                                                                                                        ctGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTc                                          <  1:3390088/151‑1 (MQ=255)
                                                                                                         tGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCg                                         <  1:4373019/151‑1 (MQ=255)
                                                                                                             ggCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCtt                                     <  1:3916442/151‑1 (MQ=255)
                                                                                                                 cccAGCCGCCAGGCACTGGGCTTTGCCCTGACCGGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGCCAACTGCCCTGGATTTCGTCTTGCtc                                 >  2:3966986/1‑150 (MQ=255)
                                                                                                                 cccAGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTg                                 <  1:2242666/151‑1 (MQ=255)
                                                                                                                    aGCCGCCAGGCACTGGGCTTTGCCCTGGCGAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAGACACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAAc                              <  2:2406492/151‑1 (MQ=255)
                                                                                                                    aGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAAc                              <  2:5153286/151‑1 (MQ=255)
                                                                                                                    aGCCGCCAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAAc                              <  1:2708396/151‑1 (MQ=255)
                                                                                                                     gccgccAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACg                             >  1:4726599/1‑151 (MQ=255)
                                                                                                                         ccAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCt                         <  1:3302543/151‑1 (MQ=255)
                                                                                                                          cAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAGAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGACCGCCCTa                        >  2:3732126/1‑151 (MQ=255)
                                                                                                                          cAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTa                        <  2:290163/151‑1 (MQ=255)
                                                                                                                          cAGGCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTa                        >  1:4380253/1‑151 (MQ=255)
                                                                                                                            ggCACTGGGCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGGTTCGTCTTGCTGACCGCCCTAcc                      >  2:4839624/1‑151 (MQ=255)
                                                                                                                                ttGGGCTTTGCGCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCAAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCtt                  <  2:1092382/150‑1 (MQ=255)
                                                                                                                                 ggggCTGTGCCGGGAGCAGGCGGGGGAGCGTGCTGGGTGCGCCGGAATGGCGGTTGTTCTACTGAAAGAGATGCAAGATACGCGGTGCTGCCCATAAAAAGACGCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTCCCGCttt                 <  2:1050307/150‑1 (MQ=255)
                                                                                                                                  gggCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTg                <  1:4007507/151‑1 (MQ=255)
                                                                                                                                   ggCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGa               >  2:2818804/1‑151 (MQ=255)
                                                                                                                                    gCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGAt              >  1:3289579/1‑151 (MQ=255)
                                                                                                                                    gCTTTGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGAt              <  2:2579739/151‑1 (MQ=255)
                                                                                                                                        tGCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGa          >  1:895247/1‑151 (MQ=255)
                                                                                                                                         gCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGccc                          >  1:304025/1‑134 (MQ=255)
                                                                                                                                         gCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGaa         >  1:5221221/1‑151 (MQ=255)
                                                                                                                                         gCCCTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGaa         <  1:4188843/151‑1 (MQ=255)
                                                                                                                                          ccatGACCAGGCTGGAAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGCTCAGAAc        <  1:1174817/148‑1 (MQ=255)
                                                                                                                                            cTGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACgg      >  1:3492093/1‑151 (MQ=255)
                                                                                                                                             tGACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACGGa     <  1:3707365/151‑1 (MQ=255)
                                                                                                                                              gACCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAGCCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACGGaa    <  2:479276/151‑1 (MQ=255)
                                                                                                                                               aCCAGGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGGTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACGGAat   <  2:624161/151‑1 (MQ=255)
                                                                                                                                                ccAGGCCGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTCCTACTGAAAGCGATGCCAGATACGCGGTGCTCCCCATAAAACAACCCAGGGCAACTGCCCTGGGTTCCGTCCTGCTGATCGCCCTACCGCGTAGACCAGCACggggta  >  2:3339172/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGACGAATTTGCCAACCTGGCACGTGACAACGGCGTAACAGTGGTCAGCGGCAGCCAGTTCCTGCCAGAACGCAGCCACGCCACCAGCGGTGTGCGCATCGCTCTGATGGCACCCAGCCGCCAGGCACTGGGCTTTGCCCTGACCAAGCTGGCAAGCCTGCTGGATTCGCCGGAACCGCGGTTGTTCTACTGAAAGCGATGCCAGATACGCGGTGCTGCCCATAAAAAAACCCAGGGCAACTGCCCTGGGTTTCGTCTTGCTGAACGCCCTACCGCTTTGATCAGAACGGAATA  >  CP009974/1532876‑1533170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: