Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,401,341 C→A 7.1% D238Y (GAT→TAT)  ydaM ← diguanylate cyclase, csgD regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,401,3410CA 7.1% 932.5 / 6.7 311D238Y (GAT→TAT) ydaMdiguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base C (113/159);  new base A (37/2);  total (150/161)
Fisher's exact test for biased strand distribution p-value = 6.42e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

AAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAG  >  CP009273/1401193‑1401490
                                                                                                                                                    |                                                                                                                                                     
aaGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCa                                                                                                                                                      <  1:4599128/150‑1 (MQ=255)
   cTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGc                                                                                                                                                   >  2:1367653/1‑150 (MQ=255)
    tGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGca                                                                                                                                                  >  2:5278334/1‑150 (MQ=255)
      tAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGcaca                                                                                                                                                <  2:1455137/150‑1 (MQ=255)
       aaTCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAATATCATGCACaa                                                                                                                                               >  1:933883/1‑150 (MQ=255)
            tGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACAt                                                                                                                                          <  1:2004132/150‑1 (MQ=255)
             gAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATa                                                                                                                                         <  2:288868/150‑1 (MQ=255)
              aGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAAAACATaa                                                                                                                                        >  2:2246411/1‑150 (MQ=255)
               gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAAc                                                                                                                                       <  1:1477848/150‑1 (MQ=255)
                 gATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACAt                                                                                                                                     <  2:545628/150‑1 (MQ=255)
                  aTGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATg                                                                                                                                    <  2:2102543/150‑1 (MQ=255)
                  aTGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCTCAGTAATATCATGCACAATACATAACATg                                                                                                                                    >  1:2184170/1‑150 (MQ=255)
                  aTGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCGCAGTAATATAATGCACAATACATAACAtt                                                                                                                                    >  1:1696231/1‑149 (MQ=255)
                    ggCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAg                                                                                                                                  <  2:1216412/150‑1 (MQ=255)
                    ggCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAAAAACATGCACAAAACATaaaaagat                                                                                                                                  >  1:3765201/1‑144 (MQ=255)
                    ggCGAGATGCTGAATTTTGACTGGTTCCGTAATGTTATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATTCACAATACATAACATGAg                                                                                                                                  <  1:4257894/150‑1 (MQ=255)
                       gagaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCtt                                                                                                                               <  2:160705/150‑1 (MQ=255)
                       gagaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCtt                                                                                                                               <  2:1003074/150‑1 (MQ=255)
                         gaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATTCACAATACATAACATGAGCTTGt                                                                                                                             >  1:1189349/1‑150 (MQ=255)
                         gaTGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGt                                                                                                                             >  2:3341374/1‑150 (MQ=255)
                           tgctgcATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCg                                                                                                                           <  2:1075989/150‑1 (MQ=255)
                                cATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATa                                                                                                                      >  1:3058851/1‑150 (MQ=255)
                                  tttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATaaa                                                                                                                    <  1:410772/150‑1 (MQ=255)
                                   ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACAAAACATGAGCTTGTCGCCATAAAt                                                                                                                   >  1:3399015/1‑150 (MQ=255)
                                   ttGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCTCAGTAATATCATGCACCATACATAACATGTGCTTGTCGCCATAAAt                                                                                                                   >  1:2542995/1‑150 (MQ=255)
                                     ggCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCTCAGTAATATCATGCACAATACATAACAAGAGCTTGTCGCCATAAAttt                                                                                                                 >  1:334102/1‑150 (MQ=255)
                                        cTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCaa                                                                                                              >  2:156352/1‑150 (MQ=255)
                                          ggTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATc                                                                                                            >  2:4716219/1‑150 (MQ=255)
                                             tCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTTGTGAGTGATATAAAACACAAAACAAAACCAGAGCTGGTCGCCCAAAAATTAAAtatgg                                                                                                         >  1:2267267/1‑146 (MQ=255)
                                               cGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTcc                                                                                                       <  2:3501116/150‑1 (MQ=255)
                                                gTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCg                                                                                                      <  2:3394972/150‑1 (MQ=255)
                                                     gtgATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATa                                                                                                 <  2:4586736/150‑1 (MQ=255)
                                                      tgATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATTCACAATACATAACATGAGCTTGTCGCCATAAAATTCAATCGGTCCGGCATAg                                                                                                >  1:2674901/1‑150 (MQ=255)
                                                       gATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTGTCAGTAAAATCAATCACAAAAAAAAAAAAGAGCTTGTCGCCCAAAAATTAAAACGGGCCGGccaagg                                                                                               >  1:4962214/1‑145 (MQ=255)
                                                       gATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAATATAATTCACAATACATAAAATGAGCTTGTCGCCAAAAATTTCAATCGGTCCGGCATAgg                                                                                               >  1:2068354/1‑150 (MQ=255)
                                                         tAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTc                                                                                             >  2:2790190/1‑150 (MQ=255)
                                                          aaaaCTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCt                                                                                            <  1:5324913/150‑1 (MQ=255)
                                                           aaaCTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTg                                                                                           <  1:2610555/150‑1 (MQ=255)
                                                            caCTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGc                                                                                          <  2:835278/149‑1 (MQ=255)
                                                              cTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATAATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGcac                                                                                        <  2:5087061/150‑1 (MQ=255)
                                                                gTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACAt                                                                                      <  2:4501763/150‑1 (MQ=255)
                                                                   gCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGGTCAGTAATATAATGCACAATAAATAACAAGAGCTTGTAGCCAAAAATTTAAATAGGTCCGGCATAGGGCTGCACAAGAc                                                                                   >  1:3635711/1‑150 (MQ=255)
                                                                           aGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATTCACAATACAAAACAAGAGCTTGTCGCCCAAAATTTCAATCGGTCCGGCCTATGTCTGCCCAAGACGAGTCGaa                                                                           >  1:1127286/1‑150 (MQ=255)
                                                                            gTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAAc                                                                          <  1:1420527/150‑1 (MQ=255)
                                                                            gTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCACTAATATAAATCACAATACAAAACAAGAGCGTGTCGCCCAAAATTTAAAACGGTCCGGCATAGGTCTGCACAAGACGAGTAGAAc                                                                          >  1:2720109/1‑150 (MQ=255)
                                                                             tAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAAcc                                                                         <  1:4510996/150‑1 (MQ=255)
                                                                               aTCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCAt                                                                       <  1:5154339/150‑1 (MQ=255)
                                                                               aTCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGGTGAGGAATATAAAGCACAAAACAAAACAAGAGCTTGTCGCCCAAAATTTCAAAAGGTGCGGCCGAGGTCTGCACAAGACGAGTCGAACCaa                                                                       >  1:1365164/1‑149 (MQ=255)
                                                                                 ccGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATcc                                                                     <  1:692054/150‑1 (MQ=255)
                                                                                 ccGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATcc                                                                     <  2:316822/150‑1 (MQ=255)
                                                                                 ccGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCGTGTAGCCATAAATTTAAATCGGTCCGGCCAAGGGCTGCACATGTCGAGTCGaaccaacc                                                                     >  1:4401932/1‑150 (MQ=255)
                                                                                 ccGGTCATCGCGTAATGGTGTGAAGCATGTTCCAGCTTCTCCTACAGCCTTTTTTGCTCAGTAATATAATGAACAATACATAACCTGAGCTTGTCGCCATTAATTTCAATCGGTCCGGCATAGGTCTGCACATGAAGAGTCGAACCATcc                                                                     <  2:3990870/150‑1 (MQ=255)
                                                                                   ggTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAATATCATTCACAATACATAACATGAGCTTTTCGCCATAAATTTAAAACGGTCCGGCAAAGGTCTGCACAAGACGAGTCGAACCATCCGc                                                                   >  1:3377568/1‑150 (MQ=255)
                                                                                    gTCATCGCGTCATGGTGAGCAGCTTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGcc                                                                  <  2:1768946/150‑1 (MQ=255)
                                                                                    gTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGGTGAGTAATATAATGCACAATACATAACAAGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATTTCGTGTCGGACCATCCGcc                                                                  >  1:52554/1‑150 (MQ=255)
                                                                                    gTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGCGCAGTAAAATAAAGCACAAAAAAAAACAAGAGCTTGTAGCCAAAAAATTCAAACGGGACGGCATAGGTGTGCACATGACGAGTCGAACCaaacccc                                                                  >  2:1812803/1‑144 (MQ=255)
                                                                                       aTCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGt                                                               <  2:776569/150‑1 (MQ=255)
                                                                                       aTCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGt                                                               <  2:4343973/150‑1 (MQ=255)
                                                                                        tCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                                                              <  1:919229/150‑1 (MQ=255)
                                                                                        tCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGtt                                                              >  2:2999733/1‑150 (MQ=255)
                                                                                         cgcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                                                             >  2:91818/1‑150 (MQ=255)
                                                                                         cgcgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCAGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGttt                                                             <  1:3908386/150‑1 (MQ=255)
                                                                                           cgTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTAt                                                           <  2:4287777/150‑1 (MQ=255)
                                                                                            gTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGGTCAGTAATATAATGCACAAAACAAAACATGAGCTTGTCGCCAGAAATTTAAATCGGTCCGGCCTAGGTCTGCACAAGTCGCGTCGAACCATCCGCCAGGTTATg                                                          >  1:4666769/1‑150 (MQ=255)
                                                                                              cATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGGTAAGTAAAATAAATCACAATACATAAAATGAGCGTTTCGCCAGAAATTTCAATCGGTCCGGCAAAGGTCTGCACAAGAAGAGTCGAACCATACGCCAGTTTATGaa                                                        >  1:5189727/1‑150 (MQ=255)
                                                                                               aTGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAAc                                                       <  2:401997/150‑1 (MQ=255)
                                                                                               aTGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTGTTTGGTCAGTAATATAATTCACAAAAAATAAAAAGAGCTTGTGGCCAAAAAATTAAATAGGTCCGGCAGAGGTCTGCACATGACGAGTCGAACCAACCGCCAGTTTATTAAc                                                       >  1:2421594/1‑150 (MQ=255)
                                                                                                 ggTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACaa                                                     <  2:4044596/150‑1 (MQ=255)
                                                                                                  gTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACaaa                                                    <  1:3005625/150‑1 (MQ=255)
                                                                                                     agcagcATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAAtt                                                 <  2:3284882/150‑1 (MQ=255)
                                                                                                        agcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAATATCATGCACAAAACATAACATGAGCTTGTCGCCAAAAAATTAAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaa                                              >  1:2627871/1‑150 (MQ=255)
                                                                                                        agcaTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATAGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCaa                                              >  2:4738027/1‑150 (MQ=255)
                                                                                                           aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAgg                                           <  1:1880593/150‑1 (MQ=255)
                                                                                                           aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAgg                                           <  2:3885134/150‑1 (MQ=255)
                                                                                                           aTGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATAATGCACAAAACATAACAAGAGCTTGTCGCCAAAAATTTAAATCGGTCCGGCATAGGTCTGCACAAGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAgg                                           >  1:2108781/1‑150 (MQ=255)
                                                                                                            tGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCCTAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCCGTTTATGAACAAAATTCAAAGGt                                          >  2:200114/1‑150 (MQ=255)
                                                                                                              ttCCAGCTGCTCCTCCAGCCGTTTTTTCTCAGTAAAATAAATCACAATACATAAAAAGAGCATGTGGCCATAAATTTAAAAAGGTCCGGCAAATGTCTGCACATGACGAGTAGAAACATACGCCAGTTTATGAAAAAAATTAAAAAGttt                                        >  1:48628/1‑150 (MQ=255)
                                                                                                              ttCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGGAATATCATGCACAATACATAACATTAGCTTGTAGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGATTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGgtt                                        >  2:4128496/1‑150 (MQ=255)
                                                                                                                ccAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTAt                                      <  2:3058851/150‑1 (MQ=255)
                                                                                                                ccAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTAt                                      <  1:420244/150‑1 (MQ=255)
                                                                                                                   gctgctCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGac                                   <  1:4446685/150‑1 (MQ=255)
                                                                                                                     tgctCCTCCAGCCGTTTTTGCTCAGTAAAATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTAAAAGGTTTATGacca                                 >  2:2423112/1‑150 (MQ=255)
                                                                                                                            ccAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTa                          <  2:5118694/150‑1 (MQ=255)
                                                                                                                             cAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTaa                         <  2:1151098/150‑1 (MQ=255)
                                                                                                                             cAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTaa                         <  1:4754676/150‑1 (MQ=255)
                                                                                                                             cAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTaa                         <  1:5087199/150‑1 (MQ=255)
                                                                                                                             cAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTaa                         <  1:2045606/150‑1 (MQ=255)
                                                                                                                               gCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAAt                       <  2:3765201/150‑1 (MQ=255)
                                                                                                                                 cGTTTTTGCTCATTAATAACATTCAAAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGAAGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAATGGTTTATGACCACCGGGTAAATGc                     <  1:1615700/150‑1 (MQ=255)
                                                                                                                                            aGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatg          <  1:5286923/150‑1 (MQ=255)
                                                                                                                                             gTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgc         >  1:1224479/1‑150 (MQ=255)
                                                                                                                                               aaTATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgcat       <  1:3120344/150‑1 (MQ=255)
                                                                                                                                               aaTATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATACGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgcat       <  1:1040019/150‑1 (MQ=255)
                                                                                                                                                atatCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgcatg      <  1:3223578/150‑1 (MQ=255)
                                                                                                                                                atatCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTcatgcatg      <  1:1837608/150‑1 (MQ=255)
                                                                                                                                                 tatCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGa     >  1:2931781/1‑150 (MQ=255)
                                                                                                                                                  atCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGAt    <  2:3635711/150‑1 (MQ=255)
                                                                                                                                                  atCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGAt    >  2:2552154/1‑150 (MQ=255)
                                                                                                                                                  atCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGAt    >  1:4806520/1‑150 (MQ=255)
                                                                                                                                                   tCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGAGTTCATGCATGATa   <  2:748143/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGCTGTAATCCTGAGCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAATCCGGTCATCGCGTCATGGTGAGCAGCATGTTCCAGCTGCTCCTCCAGCCGTTTTTGCTCAGTAATATCATGCACAATACATAACATGAGCTTGTCGCCATAAATTTCAATCGGTCCGGCATAGGTCTGCACATGACGAGTCGAACCATCCGCCAGTTTATGAACAAAATTCAAAGGTTTATGACCACCGGGTAAATGCGAGATTTCATGCATGATAG  >  CP009273/1401193‑1401490

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