Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,327,148 T→G 15.8% V615G (GTA→GGA)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,327,1480TG15.8% 642.3 / 126.8 259V615G (GTA→GGA) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base T (130/88);  new base G (22/19);  total (152/107)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 261 total aligned reads displayed.

TCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAAT  >  CP009273/1326999‑1327297
                                                                                                                                                     |                                                                                                                                                     
 cACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTa                                                                                                                                                      >  2:4298431/1‑150 (MQ=255)
   cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTAtt                                                                                                                                                    <  1:1795535/150‑1 (MQ=255)
    ttcttcTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGGGGTCGCAAAATGGGGATTCGCACAGCGAGCCCCGGGGGATTc                                                                                                                                                   >  2:3092976/1‑150 (MQ=255)
          tcGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGc                                                                                                                                             <  2:557211/150‑1 (MQ=255)
          tcGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGc                                                                                                                                             <  1:1625768/150‑1 (MQ=255)
           cGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCt                                                                                                                                            >  2:1386963/1‑150 (MQ=255)
           cGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCt                                                                                                                                            >  1:2753315/1‑150 (MQ=255)
            ggATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTg                                                                                                                                           >  1:1462283/1‑150 (MQ=255)
            ggATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTg                                                                                                                                           >  2:1272437/1‑150 (MQ=255)
             gATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGt                                                                                                                                          >  1:4394168/1‑150 (MQ=255)
             gATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGt                                                                                                                                          >  1:2185417/1‑150 (MQ=255)
                ttCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCt                                                                                                                                       <  2:3977049/150‑1 (MQ=255)
                  cACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTgg                                                                                                                                     <  1:3474946/150‑1 (MQ=255)
                      cagcagTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTAt                                                                                                                                 >  1:969334/1‑150 (MQ=255)
                       agcagTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATg                                                                                                                                >  1:1955530/1‑150 (MQ=255)
                         cagTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcg                                                                                                                              >  1:2281758/1‑150 (MQ=255)
                         cagTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcg                                                                                                                              <  2:2638578/150‑1 (MQ=255)
                          agTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcgc                                                                                                                             >  2:1638583/1‑150 (MQ=255)
                           gTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCt                                                                                                                            <  1:2977207/150‑1 (MQ=255)
                            ttAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTg                                                                                                                           >  2:3582610/1‑150 (MQ=255)
                             tAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgc                                                                                                                          >  1:2643369/1‑150 (MQ=255)
                               gATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTgccg                                                                                                                        <  1:1994561/150‑1 (MQ=255)
                                aTAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgc                                                                                                                       >  2:1353201/1‑150 (MQ=255)
                                 tAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgcc                                                                                                                      >  2:2718157/1‑150 (MQ=255)
                                 tAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgcc                                                                                                                      >  1:2370993/1‑150 (MQ=255)
                                  aaaGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGccgccg                                                                                                                     >  1:2507891/1‑150 (MQ=255)
                                   aaGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGa                                                                                                                    >  2:2246284/1‑150 (MQ=255)
                                    aGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGaa                                                                                                                   >  2:4045241/1‑150 (MQ=255)
                                       tGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAaga                                                                                                                >  2:723840/1‑150 (MQ=255)
                                        gAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAagag                                                                                                               <  1:4297640/150‑1 (MQ=255)
                                         aaaaaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGc                                                                                                              <  1:259900/150‑1 (MQ=255)
                                          aaaaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAAGGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCg                                                                                                             >  2:4483270/1‑150 (MQ=255)
                                           aaaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGt                                                                                                            >  1:4237218/1‑150 (MQ=255)
                                            aaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                                                                                           <  2:2579805/150‑1 (MQ=255)
                                            aaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                                                                                           <  2:2019354/150‑1 (MQ=255)
                                            aaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                                                                                           <  2:2153242/150‑1 (MQ=255)
                                            aaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGtt                                                                                                           <  2:2226516/150‑1 (MQ=255)
                                                tCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaa                                                                                                       <  2:1700551/150‑1 (MQ=255)
                                                tCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaa                                                                                                       >  2:668893/1‑150 (MQ=255)
                                                tCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaa                                                                                                       <  2:1726508/150‑1 (MQ=255)
                                                  cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                                                                                     <  1:1717282/150‑1 (MQ=255)
                                                  cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                                                                                     <  1:1424743/150‑1 (MQ=255)
                                                  cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                                                                                     <  1:2022417/150‑1 (MQ=255)
                                                  cGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCACACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaaaa                                                                                                     >  1:4071839/1‑150 (MQ=255)
                                                    gAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacc                                                                                                   >  2:195677/1‑150 (MQ=255)
                                                     aaGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAacca                                                                                                  >  1:3173316/1‑150 (MQ=255)
                                                             ggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAc                                                                                          <  2:853789/150‑1 (MQ=255)
                                                             ggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAc                                                                                          >  2:2524919/1‑150 (MQ=255)
                                                              gtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                                                                                         <  2:950049/150‑1 (MQ=255)
                                                                 tGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTgg                                                                                      <  1:2687225/150‑1 (MQ=255)
                                                                      ccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCg                                                                                 >  1:888462/1‑150 (MQ=255)
                                                                          aCCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaa                                                                             <  1:1053458/150‑1 (MQ=255)
                                                                           ccAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaaa                                                                            >  1:3229174/1‑150 (MQ=255)
                                                                            cAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAAc                                                                           <  2:158406/150‑1 (MQ=255)
                                                                             aGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg                                                                          >  2:4452196/1‑150 (MQ=255)
                                                                             aGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg                                                                          >  1:4122927/1‑150 (MQ=255)
                                                                             aGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg                                                                          >  2:2802771/1‑150 (MQ=255)
                                                                             aGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg                                                                          >  1:2956303/1‑150 (MQ=255)
                                                                              gATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGa                                                                         >  1:125284/1‑150 (MQ=255)
                                                                                  gTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTg                                                                     <  1:1291480/150‑1 (MQ=255)
                                                                                  gTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTg                                                                     >  2:4391575/1‑150 (MQ=255)
                                                                                   ttCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGc                                                                    >  2:3000998/1‑150 (MQ=255)
                                                                                     cTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTg                                                                  >  1:2189012/1‑150 (MQ=255)
                                                                                      tGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                 <  1:708352/150‑1 (MQ=255)
                                                                                       gACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGaa                                                                >  1:3773022/1‑150 (MQ=255)
                                                                                        aCCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGGTGCAAAACCACCATTAACCTGGTGCCGGAAAACGACGTGCTGAAc                                                               >  2:3357909/1‑150 (MQ=255)
                                                                                          cAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGt                                                             <  2:375307/150‑1 (MQ=255)
                                                                                          cAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGt                                                             >  1:2406419/1‑150 (MQ=255)
                                                                                             cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                          >  2:2900373/1‑150 (MQ=255)
                                                                                             cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                          >  2:697503/1‑150 (MQ=255)
                                                                                             cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATACCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                          >  1:288093/1‑150 (MQ=255)
                                                                                               ttGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTgg                                                        >  1:3822273/1‑150 (MQ=255)
                                                                                                 gACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGaa                                                      >  1:212917/1‑150 (MQ=255)
                                                                                                  aCTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAg                                                     <  1:238472/150‑1 (MQ=255)
                                                                                                     gCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCg                                                  <  2:3804923/150‑1 (MQ=255)
                                                                                                           cTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTAGTCCTGGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATg                                            <  2:1823227/150‑1 (MQ=255)
                                                                                                            ttGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGc                                           <  1:2802771/150‑1 (MQ=255)
                                                                                                            ttGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGc                                           >  1:1204043/1‑150 (MQ=255)
                                                                                                               tgGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGa                                        >  2:502844/1‑150 (MQ=255)
                                                                                                                   cGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAAcc                                    >  2:3573494/1‑150 (MQ=255)
                                                                                                                   cGCAAAATGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAAcc                                    >  2:2374052/1‑150 (MQ=255)
                                                                                                                      aaaaTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAc                                 >  2:3544755/1‑150 (MQ=255)
                                                                                                                         aTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAcgcg                              >  1:3086590/1‑150 (MQ=255)
                                                                                                                         aTGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAcgcg                              <  1:2316888/150‑1 (MQ=255)
                                                                                                                         aTGGGGATTCGCACAGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAAcgcg                              <  2:1068116/150‑1 (MQ=255)
                                                                                                                          tGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACgcgc                             >  2:3341237/1‑150 (MQ=255)
                                                                                                                           ggggATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCt                            >  2:233536/1‑150 (MQ=255)
                                                                                                                               aTTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGACACCAACGCGCTgcgc                        >  2:2964461/1‑150 (MQ=255)
                                                                                                                                ttCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTgcgcg                       >  1:648890/1‑150 (MQ=255)
                                                                                                                                ttCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTgcgcg                       <  2:3133/150‑1 (MQ=255)
                                                                                                                                       aGCGAGCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAAcg                >  1:3555773/1‑150 (MQ=255)
                                                                                                                                          gAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAAcgtcg             >  1:734311/1‑150 (MQ=255)
                                                                                                                                            gCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGtt           <  2:323999/150‑1 (MQ=255)
                                                                                                                                            gCACCGGGGGATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGtt           >  2:1556564/1‑150 (MQ=255)
                                                                                                                                             cACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTg          >  1:4043475/1‑150 (MQ=255)
                                                                                                                                              aCCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGc         <  2:1482508/150‑1 (MQ=255)
                                                                                                                                              aCCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGc         <  1:4506939/150‑1 (MQ=255)
                                                                                                                                                cGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGccc       >  2:969271/1‑150 (MQ=255)
                                                                                                                                                cGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGccc       >  1:1094792/1‑150 (MQ=255)
                                                                                                                                                   ggTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGaa    >  1:195721/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAAT  >  CP009273/1326999‑1327297

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: ‑