Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,481,732 G→A 15.9% A81T (GCC→ACC)  ydcF → DUF218 superfamily protein, SAM‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,481,7320GA15.9% 711.6 / 109.3 257A81T (GCC→ACC) ydcFDUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base G (112/104);  new base A (18/23);  total (130/127)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 261 total aligned reads displayed.

GTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAACTGCGGTGAAAACGC  >  CP009273/1481583‑1481880
                                                                                                                                                     |                                                                                                                                                    
    ggTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGAGTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCAt                                                                                                                                                  >  2:2966589/1‑150 (MQ=255)
     gTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATc                                                                                                                                                 <  1:253959/150‑1 (MQ=255)
       gCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGc                                                                                                                                               >  1:2749209/1‑150 (MQ=255)
        ccGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGca                                                                                                                                              >  1:240588/1‑150 (MQ=255)
          gTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGcaca                                                                                                                                            >  2:675355/1‑150 (MQ=255)
           tATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAg                                                                                                                                           <  1:1988628/150‑1 (MQ=255)
             tCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCa                                                                                                                                         <  2:3315511/150‑1 (MQ=255)
               aGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATc                                                                                                                                       >  2:3025577/1‑150 (MQ=255)
                ggCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATcc                                                                                                                                      >  2:1823054/1‑150 (MQ=255)
                   cGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCa                                                                                                                                   <  2:3188352/150‑1 (MQ=255)
                    gATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCAc                                                                                                                                  >  1:4303538/1‑150 (MQ=255)
                      ttGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTa                                                                                                                                >  1:3837318/1‑150 (MQ=255)
                          gTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaac                                                                                                                            >  1:2859663/1‑150 (MQ=255)
                           tGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTacaaca                                                                                                                           <  2:2594425/150‑1 (MQ=255)
                            gATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAacac                                                                                                                          <  1:1397134/150‑1 (MQ=255)
                            gATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAacac                                                                                                                          <  2:2219992/150‑1 (MQ=255)
                            gATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAacac                                                                                                                          <  2:2220015/150‑1 (MQ=255)
                              tCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTa                                                                                                                        <  2:3700562/150‑1 (MQ=255)
                                cTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATc                                                                                                                      >  2:2770009/1‑150 (MQ=255)
                                cTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATc                                                                                                                      >  2:2565518/1‑150 (MQ=255)
                                  tGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCg                                                                                                                    >  1:3442360/1‑150 (MQ=255)
                                              ggTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcag                                                                                                        >  1:2198315/1‑150 (MQ=255)
                                               gTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcaga                                                                                                       <  1:1229774/150‑1 (MQ=255)
                                                 tATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagc                                                                                                     <  1:3742994/150‑1 (MQ=255)
                                                  aTGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGgcagagca                                                                                                    <  2:347248/150‑1 (MQ=255)
                                                       gACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGc                                                                                               <  1:4183079/150‑1 (MQ=255)
                                                         cTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGa                                                                                             <  2:4146804/150‑1 (MQ=255)
                                                            tCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCa                                                                                          >  2:3243617/1‑150 (MQ=255)
                                                            tCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCa                                                                                          >  2:2110656/1‑150 (MQ=255)
                                                            tCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCa                                                                                          >  1:2824698/1‑150 (MQ=255)
                                                              gATGCGTCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATc                                                                                        <  2:3584279/150‑1 (MQ=255)
                                                               aTGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATcc                                                                                       <  2:708179/150‑1 (MQ=255)
                                                               aTGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATcc                                                                                       >  1:3848308/1‑150 (MQ=255)
                                                               aTGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATcc                                                                                       >  2:3810381/1‑150 (MQ=255)
                                                                  cGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTgg                                                                                    >  1:366171/1‑150 (MQ=255)
                                                                   ggCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGc                                                                                   <  2:1027339/150‑1 (MQ=255)
                                                                    gCATGTAAGATGGCCCGCGATCAGCAATTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCg                                                                                  <  1:3276160/150‑1 (MQ=255)
                                                                       tGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGat                                                                               >  2:4350901/1‑150 (MQ=255)
                                                                       tGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGat                                                                               >  1:1750705/1‑150 (MQ=255)
                                                                       tGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGat                                                                               >  2:3661946/1‑150 (MQ=255)
                                                                       tGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGat                                                                               <  2:2832197/150‑1 (MQ=255)
                                                                         tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGatat                                                                             >  2:2575128/1‑150 (MQ=255)
                                                                         tAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGatat                                                                             <  2:1000618/150‑1 (MQ=255)
                                                                          aaGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATc                                                                            >  2:2324296/1‑150 (MQ=255)
                                                                          aaGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATc                                                                            >  2:2324284/1‑150 (MQ=255)
                                                                          aaGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATc                                                                            >  1:2080655/1‑150 (MQ=255)
                                                                           aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCg                                                                           <  2:3022192/150‑1 (MQ=255)
                                                                           aGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCg                                                                           <  1:750934/150‑1 (MQ=255)
                                                                             aTTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCt                                                                         <  1:3249516/150‑1 (MQ=255)
                                                                               tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtca                                                                       >  1:622852/1‑150 (MQ=255)
                                                                               tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtca                                                                       >  1:4285157/1‑150 (MQ=255)
                                                                               tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtca                                                                       >  1:3058582/1‑150 (MQ=255)
                                                                               tGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtca                                                                       >  1:3512131/1‑150 (MQ=255)
                                                                                  ccGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtcatca                                                                    <  1:2886705/150‑1 (MQ=255)
                                                                                   cgcgATCAGCAAATTCCTTTCCTGATTTGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATCGCGCCTTCGCAAAGCATCCGCACTACACCACTATCCGCACCACTGGCAGAGAAGAAGCGACCATCCTGGCGGATATCGCTCATCAg                                                                   <  1:3047464/150‑1 (MQ=255)
                                                                                   cgcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAg                                                                   <  1:3629733/150‑1 (MQ=255)
                                                                                   cgcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAg                                                                   >  2:116031/1‑150 (MQ=255)
                                                                                    gcgATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGt                                                                  <  2:890759/150‑1 (MQ=255)
                                                                                       aTCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCt                                                               <  2:405300/150‑1 (MQ=255)
                                                                                       aTCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCt                                                               >  1:2261033/1‑150 (MQ=255)
                                                                                       aTCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCt                                                               <  2:2669881/150‑1 (MQ=255)
                                                                                         cAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTgg                                                             <  1:675355/150‑1 (MQ=255)
                                                                                          aGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGc                                                            >  2:1749729/1‑150 (MQ=255)
                                                                                           gCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGca                                                           >  2:495096/1‑150 (MQ=255)
                                                                                           gCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGca                                                           <  2:553454/150‑1 (MQ=255)
                                                                                              aaTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACAt                                                        >  2:1392483/1‑150 (MQ=255)
                                                                                               aTTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACAtt                                                       <  2:3479889/150‑1 (MQ=255)
                                                                                                 tCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTcc                                                     <  1:1289790/150‑1 (MQ=255)
                                                                                                  ccTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACCCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCg                                                    <  1:3174409/150‑1 (MQ=255)
                                                                                                    tttACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa                                                  <  2:3837062/150‑1 (MQ=255)
                                                                                                      tACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATg                                                >  2:2533613/1‑150 (MQ=255)
                                                                                                       aCTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGa                                               <  1:3991058/150‑1 (MQ=255)
                                                                                                           aTTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGaaaaa                                           >  1:536425/1‑150 (MQ=255)
                                                                                                             tAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAAt                                         <  1:3649456/150‑1 (MQ=255)
                                                                                                                tggtggTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTg                                      <  2:204518/150‑1 (MQ=255)
                                                                                                                 ggtggtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTgg                                     <  2:3783613/150‑1 (MQ=255)
                                                                                                                  gtggtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGa                                    >  1:3645241/1‑150 (MQ=255)
                                                                                                                   tggtATCGGTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGAt                                   <  1:1517111/150‑1 (MQ=255)
                                                                                                                     gtATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg                                 >  1:828145/1‑150 (MQ=255)
                                                                                                                        tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCTCATTCCGCATGAAAAAATCTGGATTGAAg                              <  2:1366902/150‑1 (MQ=255)
                                                                                                                        tCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAg                              <  2:2617225/150‑1 (MQ=255)
                                                                                                                         cGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGa                             >  1:187165/1‑150 (MQ=255)
                                                                                                                          ggTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAc                            >  2:3241849/1‑150 (MQ=255)
                                                                                                                           gTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAcc                           <  1:442400/150‑1 (MQ=255)
                                                                                                                           gTCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAcc                           <  1:2163252/150‑1 (MQ=255)
                                                                                                                            tCACTCGACAACTTTTTTTTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa                          <  1:4331272/150‑1 (MQ=255)
                                                                                                                            tCACTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa                          <  2:2749187/150‑1 (MQ=255)
                                                                                                                              aCTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGt                        >  1:1206230/1‑150 (MQ=255)
                                                                                                                              aCTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGt                        >  2:109316/1‑150 (MQ=255)
                                                                                                                               cTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTc                       >  2:2607400/1‑150 (MQ=255)
                                                                                                                               cTCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTc                       >  1:1925825/1‑150 (MQ=255)
                                                                                                                                tCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTca                      >  1:1011997/1‑150 (MQ=255)
                                                                                                                                tCGACAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTca                      <  2:4303523/150‑1 (MQ=255)
                                                                                                                                   aCAACTTTTTTGTATAGCACCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTcaaca                   >  2:3847332/1‑150 (MQ=255)
                                                                                                                                    cAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTcaacaa                  >  1:4510367/1‑150 (MQ=255)
                                                                                                                                    cAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTcaacaa                  >  2:2636722/1‑150 (MQ=255)
                                                                                                                                       cTTTTTTGTATATCGCCATCGCACTGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGAACATTCCGCTTTAAAAAATCTGGATTGAAGACCAGTCAACAAACt               <  2:1752271/150‑1 (MQ=255)
                                                                                                                                       cTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAACt               <  1:3816585/150‑1 (MQ=255)
                                                                                                                                        ttttttGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAACTg              >  2:238822/1‑150 (MQ=255)
                                                                                                                                                    cgcCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAACTGCGGTGAAAacgc  >  2:2319642/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGAGGTGCCGTATCAGGCCGATTGCGTGATCCTTGCAGGCAATGCGGTTATGCCGACTATCGATGCGGCATGTAAGATTGCCCGCGATCAGCAAATTCCTTTACTGATTAGTGGTGGTATCGGTCACTCGACAACTTTTTTGTATAGCGCCATCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAACTGCGGTGAAAACGC  >  CP009273/1481583‑1481880

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: