Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,593,189 T→G 8.0% A426A (GCA→GCC lsrK ← autoinducer‑2 (AI‑2) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,593,1890TG 8.0% 920.6 / 68.3 308A426A (GCA→GCClsrKautoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base T (158/125);  new base G (9/16);  total (167/141)
Fisher's exact test for biased strand distribution p-value = 6.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 311 total aligned reads displayed.

CTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTTAATAAAGGA  >  CP009273/1593041‑1593338
                                                                                                                                                    |                                                                                                                                                     
cTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTg                                                                                                                                                      >  1:3649790/1‑150 (MQ=255)
 tGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGc                                                                                                                                                     >  2:4051293/1‑150 (MQ=255)
  gCCTTTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACATGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTTCa                                                                                                                                                    <  2:3525483/150‑1 (MQ=255)
  gCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCa                                                                                                                                                    <  2:1579662/150‑1 (MQ=255)
   ccATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCCGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATACTTTCCCTTGTGAACCTCCGCCTGCaa                                                                                                                                                   >  2:2870380/1‑150 (MQ=255)
    cATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCaaa                                                                                                                                                  >  2:725354/1‑150 (MQ=255)
       tgatgaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGAc                                                                                                                                               <  1:39628/150‑1 (MQ=255)
       tgatgaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGAc                                                                                                                                               <  2:1860326/150‑1 (MQ=255)
        gatgaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCGGCAAAGACt                                                                                                                                              >  2:1619394/1‑150 (MQ=255)
         atgaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTa                                                                                                                                             <  2:2629157/150‑1 (MQ=255)
           gaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAAc                                                                                                                                           >  1:4082437/1‑150 (MQ=255)
           gaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAAc                                                                                                                                           <  2:1024604/150‑1 (MQ=255)
           gaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAAc                                                                                                                                           <  2:4309980/150‑1 (MQ=255)
             aaaaaTTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGa                                                                                                                                         <  1:2638906/150‑1 (MQ=255)
              aaaaTTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGAt                                                                                                                                        <  2:110968/150‑1 (MQ=255)
              aaaaTTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGAt                                                                                                                                        <  1:38326/150‑1 (MQ=255)
                aaTTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCGGCAAAGACTAACGATGa                                                                                                                                      <  1:2208330/150‑1 (MQ=255)
                  ttCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAg                                                                                                                                    <  1:14705/150‑1 (MQ=255)
                     cGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGAt                                                                                                                                 <  2:592652/150‑1 (MQ=255)
                      ggCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATg                                                                                                                                >  2:2531853/1‑150 (MQ=255)
                      ggCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATg                                                                                                                                >  1:1523103/1‑150 (MQ=255)
                      ggCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCGGCAAAGACTAACGATGAAGGATg                                                                                                                                <  2:1176867/150‑1 (MQ=255)
                        cACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGaa                                                                                                                              >  2:2053903/1‑150 (MQ=255)
                         aCCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCCAAGACTAACGATGAAGGATGAat                                                                                                                             >  1:3878434/1‑150 (MQ=255)
                         aCCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAat                                                                                                                             >  2:4811186/1‑150 (MQ=255)
                          ccGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAata                                                                                                                            >  2:3149936/1‑150 (MQ=255)
                          ccGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAata                                                                                                                            <  1:429303/150‑1 (MQ=255)
                            gACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATAtt                                                                                                                          <  1:2762775/150‑1 (MQ=255)
                              cGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCg                                                                                                                        <  2:2285964/150‑1 (MQ=255)
                                ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCgag                                                                                                                      >  2:851139/1‑150 (MQ=255)
                                  aGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCGGCAAAGACTAACGATGAAGGATGAATATTCGAGaa                                                                                                                    <  2:924326/150‑1 (MQ=255)
                                      gCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCa                                                                                                                >  1:3936036/1‑150 (MQ=255)
                                        aaTGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGc                                                                                                              >  2:470918/1‑150 (MQ=255)
                                        aaTGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGc                                                                                                              >  1:4789018/1‑150 (MQ=255)
                                           ggCACATCCTTATGCTGTGGCTTCTTTGAACACCGGAATATTGACGGGTAATCCCGAGACATCAGCTAGAATTTGACTCCATAATTTCCCTTTTGTACATACGACTGCAAAGACTATCGTTGAAGGTTGAATATTCGAGAAATCAGCAAt                                                                                                           <  1:919051/150‑1 (MQ=255)
                                           ggCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCCAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAAt                                                                                                           >  2:3803354/1‑150 (MQ=255)
                                              acaTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTtg                                                                                                        >  1:101083/1‑150 (MQ=255)
                                              acaTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTtg                                                                                                        >  1:3137583/1‑150 (MQ=255)
                                              acaTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTtg                                                                                                        >  1:550856/1‑150 (MQ=255)
                                               caTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTtgc                                                                                                       <  1:2275332/150‑1 (MQ=255)
                                                 tCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTtgctg                                                                                                     <  1:2957737/150‑1 (MQ=255)
                                                   cTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCACAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCa                                                                                                   >  2:3664951/1‑150 (MQ=255)
                                                       tGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGtt                                                                                               <  1:5088337/150‑1 (MQ=255)
                                                       tGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGtt                                                                                               >  2:2705876/1‑150 (MQ=255)
                                                         cAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTac                                                                                             >  2:1320993/1‑150 (MQ=255)
                                                              gCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCt                                                                                        >  2:30401/1‑150 (MQ=255)
                                                                ttcttTGACCACCGGTATATTGTCGGGTAATCACGTGACTTCAGCGAGAATTTGACTCCATAATTTACATTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGAAATTTGCTGCAAGTTACACGCTGa                                                                                      <  2:1139595/150‑1 (MQ=255)
                                                                ttcttTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGa                                                                                      >  1:2746031/1‑150 (MQ=255)
                                                                 tcttTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGAt                                                                                     >  2:4441017/1‑150 (MQ=255)
                                                                    ttGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACa                                                                                  <  1:725354/150‑1 (MQ=255)
                                                                     tGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACaa                                                                                 >  2:4773205/1‑150 (MQ=255)
                                                                     tGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACaa                                                                                 >  2:2392966/1‑150 (MQ=255)
                                                                       accaccGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATc                                                                               >  2:341998/1‑150 (MQ=255)
                                                                        ccaccGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATcg                                                                              <  1:854069/150‑1 (MQ=255)
                                                                        ccaccGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATcg                                                                              >  1:2047815/1‑150 (MQ=255)
                                                                           ccGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATcgccg                                                                           <  2:1963083/150‑1 (MQ=255)
                                                                              gAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAt                                                                        <  1:1854719/150‑1 (MQ=255)
                                                                              gAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAt                                                                        >  1:3439775/1‑150 (MQ=255)
                                                                                atatTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAttt                                                                      >  2:1704179/1‑150 (MQ=255)
                                                                                atatTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAttt                                                                      >  1:2359701/1‑150 (MQ=255)
                                                                                 tatTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAtttt                                                                     <  2:1918239/150‑1 (MQ=255)
                                                                                  atTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttc                                                                    >  1:449124/1‑150 (MQ=255)
                                                                                    tGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCACTTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttctt                                                                  >  2:2666496/1‑150 (MQ=255)
                                                                                    tGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttctt                                                                  <  1:3050822/150‑1 (MQ=255)
                                                                                          gTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAgcg                                                            >  2:155323/1‑150 (MQ=255)
                                                                                            aaTCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCa                                                          >  2:3849550/1‑150 (MQ=255)
                                                                                               cccGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACgg                                                       >  2:4993882/1‑150 (MQ=255)
                                                                                                ccGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGa                                                      >  1:2443505/1‑150 (MQ=255)
                                                                                                    gaCATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACaa                                                  >  2:3367276/1‑150 (MQ=255)
                                                                                                     aCATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAAt                                                 <  1:582170/150‑1 (MQ=255)
                                                                                                      cATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACTATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTGTCTTCCAGCGCACGGAACAATg                                                >  2:3219389/1‑150 (MQ=255)
                                                                                                        tCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTc                                              >  2:4164994/1‑150 (MQ=255)
                                                                                                        tCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTc                                              >  2:5000401/1‑150 (MQ=255)
                                                                                                            cGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCtt                                          >  1:2901871/1‑150 (MQ=255)
                                                                                                             gagaATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCttt                                         <  2:4639354/150‑1 (MQ=255)
                                                                                                             gagaATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCttt                                         <  2:1376976/150‑1 (MQ=255)
                                                                                                             gagaATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCGGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCttt                                         <  1:4596541/150‑1 (MQ=255)
                                                                                                               gaATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGt                                       >  2:852589/1‑150 (MQ=255)
                                                                                                                 aTTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa                                     <  1:470918/150‑1 (MQ=255)
                                                                                                                     gACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACAtt                                 <  1:1167653/150‑1 (MQ=255)
                                                                                                                      aCTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACAttt                                <  1:4849522/150‑1 (MQ=255)
                                                                                                                        tCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTAt                              <  1:1282741/150‑1 (MQ=255)
                                                                                                                        tCCATAATTTCCCTTTTGAACCTACGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTAt                              <  1:1985906/150‑1 (MQ=255)
                                                                                                                          cATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATcc                            >  2:3592655/1‑150 (MQ=255)
                                                                                                                           aTAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCg                           <  1:4351557/150‑1 (MQ=255)
                                                                                                                            tAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCgg                          <  2:1495133/150‑1 (MQ=255)
                                                                                                                              aTTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGt                        >  2:4933519/1‑150 (MQ=255)
                                                                                                                               tttCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTc                       >  1:2559091/1‑150 (MQ=255)
                                                                                                                               tttCCCTTTTGAACCTCCGCCGGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTc                       >  1:1665297/1‑150 (MQ=255)
                                                                                                                                    cTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATgg                  >  1:3295403/1‑150 (MQ=255)
                                                                                                                                       ttGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACa               <  2:344210/150‑1 (MQ=255)
                                                                                                                                       ttGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACa               >  1:3382/1‑150 (MQ=255)
                                                                                                                                          aaCCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGt            <  1:1849651/150‑1 (MQ=255)
                                                                                                                                          aaCCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGt            >  2:523029/1‑150 (MQ=255)
                                                                                                                                           aCCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTGTGTTACATTTATCCGGGTCAATGGACAAGtg           >  1:3527902/1‑149 (MQ=255)
                                                                                                                                               ccGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGGCAATGGACAAGTtaata       >  2:1358601/1‑150 (MQ=255)
                                                                                                                                                cGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTtaataa      <  1:3203947/150‑1 (MQ=255)
                                                                                                                                                cGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTtaataa      >  2:2078487/1‑150 (MQ=255)
                                                                                                                                                   cTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTTAATAaagg   >  1:969097/1‑150 (MQ=255)
                                                                                                                                                   cTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTTAATAaagg   >  2:2772204/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTCAATGGACAAGTTAATAAAGGA  >  CP009273/1593041‑1593338

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: