Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,327,193 G→T 7.6% C630F (TGC→TTC)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,327,1930GT 7.6% 1194.0 / 68.2 315C630F (TGC→TTC) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (137/154);  new base T (17/7);  total (154/161)
Fisher's exact test for biased strand distribution p-value = 3.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 318 total aligned reads displayed.

AGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAG  >  CP009273/1327044‑1327341
                                                                                                                                                     |                                                                                                                                                    
 gATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGc                                                                                                                                                     >  1:4425359/1‑150 (MQ=255)
 gATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGc                                                                                                                                                     >  1:1838017/1‑150 (MQ=255)
     cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                                                                                                                                 >  1:5264554/1‑150 (MQ=255)
      ggAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                                                                                                                                                <  2:3995227/150‑1 (MQ=255)
       gAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacc                                                                                                                                               <  2:5321075/150‑1 (MQ=255)
        aaGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                                                                                                                                              >  2:3371343/1‑150 (MQ=255)
        aaGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                                                                                                                                              >  2:1345389/1‑150 (MQ=255)
         agagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccac                                                                                                                                             <  2:4671440/150‑1 (MQ=255)
         agagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccac                                                                                                                                             >  1:4727991/1‑150 (MQ=255)
         agagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAaccac                                                                                                                                             >  2:5117698/1‑150 (MQ=255)
           agGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAAccacca                                                                                                                                           >  1:1389139/1‑150 (MQ=255)
            gggTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCAt                                                                                                                                          <  2:2439327/150‑1 (MQ=255)
            gggTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCAt                                                                                                                                          <  1:2141983/150‑1 (MQ=255)
               tggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTaa                                                                                                                                       <  1:2835300/150‑1 (MQ=255)
                 gtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                                                                                                                                     <  2:5156921/150‑1 (MQ=255)
                  tATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCt                                                                                                                                    >  2:579399/1‑150 (MQ=255)
                        cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGcc                                                                                                                              <  1:3291604/150‑1 (MQ=255)
                         cctaaCCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCg                                                                                                                             >  2:774421/4‑150 (MQ=255)
                         ccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCg                                                                                                                             <  1:3389894/150‑1 (MQ=255)
                          cGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCgg                                                                                                                            <  1:5284965/150‑1 (MQ=255)
                           gAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGa                                                                                                                           <  2:4880092/150‑1 (MQ=255)
                           gAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGa                                                                                                                           <  1:1492501/150‑1 (MQ=255)
                              ccAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaaa                                                                                                                        >  2:3117715/1‑150 (MQ=255)
                               cAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAc                                                                                                                       >  2:2067634/1‑150 (MQ=255)
                                   tGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAg                                                                                                                   >  2:2250173/1‑150 (MQ=255)
                                    ggTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGt                                                                                                                  <  2:3295560/150‑1 (MQ=255)
                                    ggTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGt                                                                                                                  >  1:3572268/1‑150 (MQ=255)
                                        cTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTg                                                                                                              <  1:2680610/150‑1 (MQ=255)
                                         tGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                                             >  1:4189867/1‑150 (MQ=255)
                                            ccAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACg                                                                                                          <  1:4832079/150‑1 (MQ=255)
                                                cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                                                                      <  2:2011420/150‑1 (MQ=255)
                                                cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                                                                      <  1:1217836/150‑1 (MQ=255)
                                                 aTTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTg                                                                                                     <  2:2563007/150‑1 (MQ=255)
                                                     aCTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAg                                                                                                 >  2:5455388/1‑150 (MQ=255)
                                                      cTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAgg                                                                                                <  2:2845207/150‑1 (MQ=255)
                                                       tGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGc                                                                                               <  1:629544/150‑1 (MQ=255)
                                                        gCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCg                                                                                              >  2:5018020/1‑150 (MQ=255)
                                                           cGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAg                                                                                           >  1:1973632/1‑150 (MQ=255)
                                                            gACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGa                                                                                          >  2:2104136/1‑150 (MQ=255)
                                                            gACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGATATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGa                                                                                          >  2:3259655/1‑150 (MQ=255)
                                                              cTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATg                                                                                        <  1:2995265/150‑1 (MQ=255)
                                                                 gtgGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTg                                                                                     >  1:1754110/1‑150 (MQ=255)
                                                                      cGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAAcc                                                                                >  2:826817/1‑150 (MQ=255)
                                                                       gCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCa                                                                               <  2:1843915/150‑1 (MQ=255)
                                                                         aaaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAc                                                                             <  2:4534093/150‑1 (MQ=255)
                                                                          aaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAcg                                                                            <  1:5003484/150‑1 (MQ=255)
                                                                             tGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCGTTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACgcgc                                                                         >  1:3658376/1‑150 (MQ=255)
                                                                             tGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCGTTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACgcgc                                                                         >  1:3658609/1‑150 (MQ=255)
                                                                              ggggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGTCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCt                                                                        <  2:4113706/150‑1 (MQ=255)
                                                                               gggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTg                                                                       >  1:3462733/1‑150 (MQ=255)
                                                                                 gATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTgcg                                                                     >  2:1640705/1‑150 (MQ=255)
                                                                                  aTTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTgcgc                                                                    >  1:331634/1‑150 (MQ=255)
                                                                                      gCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCaa                                                                >  1:755647/1‑150 (MQ=255)
                                                                                          aGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAAcg                                                            >  1:2085399/1‑150 (MQ=255)
                                                                                               gCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGtt                                                       <  2:2526039/150‑1 (MQ=255)
                                                                                                  ccGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGcc                                                    <  2:4931999/150‑1 (MQ=255)
                                                                                                   cGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGccc                                                   <  1:976672/150‑1 (MQ=255)
                                                                                                    ggggTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCg                                                  >  1:244182/1‑150 (MQ=255)
                                                                                                      ggTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGaa                                                >  2:13132/1‑150 (MQ=255)
                                                                                                        tATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAACACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAAt                                              <  2:1262070/150‑1 (MQ=255)
                                                                                                         aTTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATg                                             >  2:5028976/1‑150 (MQ=255)
                                                                                                           tCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCg                                           <  2:499043/150‑1 (MQ=255)
                                                                                                              ttGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCa                                        <  1:3863838/150‑1 (MQ=255)
                                                                                                               tGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCAc                                       <  1:3618802/150‑1 (MQ=255)
                                                                                                               tGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCAc                                       <  1:3040860/150‑1 (MQ=255)
                                                                                                                ggCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACg                                      <  1:234828/150‑1 (MQ=255)
                                                                                                                   tGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCg                                   >  1:2514588/1‑150 (MQ=255)
                                                                                                                   tGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGGTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCg                                   <  2:4331447/150‑1 (MQ=255)
                                                                                                                    gTTCTGGCTATGCGCTGCCGCCGTAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGa                                  >  2:359918/1‑150 (MQ=255)
                                                                                                                     ttCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGAt                                 >  1:5171650/1‑150 (MQ=255)
                                                                                                                       cTGGCTATGCGCTGCCGCCGAAAGAGCGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATgg                               >  2:883000/1‑150 (MQ=255)
                                                                                                                       cTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATgg                               >  1:4477335/1‑150 (MQ=255)
                                                                                                                         ggCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGAc                             <  1:4069459/150‑1 (MQ=255)
                                                                                                                         ggCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGAc                             <  1:428297/150‑1 (MQ=255)
                                                                                                                                cgcTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTAtc                      >  1:5007933/1‑150 (MQ=255)
                                                                                                                                cgcTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTAtc                      <  2:4832130/150‑1 (MQ=255)
                                                                                                                                 gcTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTAtct                     <  1:961135/150‑1 (MQ=255)
                                                                                                                                   tGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCa                   <  1:5438094/150‑1 (MQ=255)
                                                                                                                                     ccgccgAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCGACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATc                 <  2:1056540/150‑1 (MQ=255)
                                                                                                                                        ccgAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGAt              >  1:4931404/1‑150 (MQ=255)
                                                                                                                                         cgAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATc             <  2:3636543/150‑1 (MQ=255)
                                                                                                                                          gAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATcc            <  2:3631373/150‑1 (MQ=255)
                                                                                                                                           aaaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCg           <  2:3533598/150‑1 (MQ=255)
                                                                                                                                            aaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGa          >  2:4446845/1‑150 (MQ=255)
                                                                                                                                            aaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGa          >  2:3470872/1‑150 (MQ=255)
                                                                                                                                            aaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGa          >  1:4652855/1‑150 (MQ=255)
                                                                                                                                            aaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGa          >  2:3806798/1‑150 (MQ=255)
                                                                                                                                            aaGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGAGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGa          <  1:4715529/150‑1 (MQ=255)
                                                                                                                                             agagCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGaa         <  2:5267482/150‑1 (MQ=255)
                                                                                                                                               agCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAAc       <  2:661535/150‑1 (MQ=255)
                                                                                                                                                 cGTTTCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGt     >  2:2291017/1‑150 (MQ=255)
                                                                                                                                                 cGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGt     <  2:1733364/150‑1 (MQ=255)
                                                                                                                                                  gTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTa    <  2:4374688/150‑1 (MQ=255)
                                                                                                                                                  gTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTa    >  2:2736568/1‑150 (MQ=255)
                                                                                                                                                  gTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTa    <  1:1188884/150‑1 (MQ=255)
                                                                                                                                                  gTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTa    <  2:1609386/150‑1 (MQ=255)
                                                                                                                                                   ttGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTaa   <  1:2083150/150‑1 (MQ=255)
                                                                                                                                                    tGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAg  <  1:3999301/150‑1 (MQ=255)
                                                                                                                                                    tGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAg  <  1:1236057/150‑1 (MQ=255)
                                                                                                                                                    tGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAg  <  1:1236012/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAG  >  CP009273/1327044‑1327341

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: