Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,413,119 G→T 26.1% pseudogene (654/960 nt) yhdW → pseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,413,1190GT26.1% 478.8 / 134.9 218pseudogene (654/960 nt)yhdWpseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (78/83);  new base T (36/21);  total (114/104)
Fisher's exact test for biased strand distribution p-value = 6.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 219 total aligned reads displayed.

TGCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCTGAAG  >  CP009273/3412973‑3413266
                                                                                                                                                  |                                                                                                                                                   
tGCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGa                                                                                                                                                    <  2:832877/148‑1 (MQ=255)
 gCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGAt                                                                                                                                                   <  1:1072418/148‑1 (MQ=255)
  cGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGAtt                                                                                                                                                  >  1:400411/1‑148 (MQ=255)
   gATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTg                                                                                                                                                 >  1:3563726/1‑148 (MQ=255)
    aTACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGt                                                                                                                                                >  2:2143433/1‑148 (MQ=255)
        gCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGc                                                                                                                                            >  2:2195032/1‑148 (MQ=255)
            gCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGa                                                                                                                                        >  2:1502878/1‑148 (MQ=255)
            gCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGa                                                                                                                                        >  2:3142661/1‑148 (MQ=255)
                cGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCt                                                                                                                                    >  1:2605451/1‑148 (MQ=255)
                    tcaatcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTc                                                                                                                                >  2:1044782/1‑148 (MQ=255)
                      aatcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGc                                                                                                                              >  2:3093101/1‑148 (MQ=255)
                       atcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGcc                                                                                                                             >  1:3092318/1‑148 (MQ=255)
                        tcaCAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCa                                                                                                                            >  1:3107740/1‑148 (MQ=255)
                          acaACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatg                                                                                                                          <  2:139871/148‑1 (MQ=255)
                          acaACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatg                                                                                                                          <  2:1198331/148‑1 (MQ=255)
                           caACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgc                                                                                                                         <  1:1346806/148‑1 (MQ=255)
                               tGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCatgctgaa                                                                                                                     >  1:644418/1‑148 (MQ=255)
                               tGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaa                                                                                                                     <  2:2600938/148‑1 (MQ=255)
                                gTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaat                                                                                                                    >  2:85761/1‑148 (MQ=255)
                                gTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaat                                                                                                                    >  1:841942/1‑148 (MQ=255)
                                 tATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatg                                                                                                                   <  1:898461/148‑1 (MQ=255)
                                   tGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatgct                                                                                                                 >  2:509693/1‑148 (MQ=255)
                                     cccTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCatgctgaatgctga                                                                                                               >  1:14227/1‑148 (MQ=255)
                                     cccTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTAGCCatgctgaattctga                                                                                                               >  1:14380/1‑148 (MQ=255)
                                        tGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAaga                                                                                                            >  2:547059/1‑148 (MQ=255)
                                          cgcATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTNGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAgaga                                                                                                          <  1:1808766/148‑1 (MQ=255)
                                          cgcATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAgaga                                                                                                          <  1:1186995/148‑1 (MQ=255)
                                              tCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATggg                                                                                                      >  2:1869510/1‑148 (MQ=255)
                                              tCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATggg                                                                                                      >  2:2005242/1‑148 (MQ=255)
                                              tCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATggg                                                                                                      >  1:1090861/1‑148 (MQ=255)
                                              tCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATggg                                                                                                      >  2:2919069/1‑148 (MQ=255)
                                              tCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATggg                                                                                                      >  2:450914/1‑148 (MQ=255)
                                                 aaTTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCAt                                                                                                   <  2:1000731/148‑1 (MQ=255)
                                                   ttAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCa                                                                                                 >  2:2208410/1‑148 (MQ=255)
                                                    tAAGCAACCCAGCTGAATGGATTGTCTTACCTGAAGTTATCTCAAAATAACCGCTTTGTCCGGTAGTTCGTCGTGTCGATGATGAATTGTTCTCTATTGTACGCTGGACGCTTTTCTCCATGCTTAATGCTGAAGAGATGGGCATCaa                                                                                                >  2:189138/1‑148 (MQ=255)
                                                    tAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCaa                                                                                                >  2:585867/1‑148 (MQ=255)
                                                     aaGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAAt                                                                                               <  2:1886152/148‑1 (MQ=255)
                                                      aGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAAtt                                                                                              >  1:217279/1‑148 (MQ=255)
                                                        cAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTcc                                                                                            <  1:2304026/148‑1 (MQ=255)
                                                        cAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTcc                                                                                            <  2:761861/148‑1 (MQ=255)
                                                           cccAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAg                                                                                         >  2:3045166/1‑148 (MQ=255)
                                                            ccAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGa                                                                                        >  1:878994/1‑148 (MQ=255)
                                                               gCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACg                                                                                     >  1:2480909/1‑148 (MQ=255)
                                                                  gAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCg                                                                                  <  1:250220/148‑1 (MQ=255)
                                                                  gAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCg                                                                                  <  2:1547805/148‑1 (MQ=255)
                                                                   aaTGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGa                                                                                 >  2:1943551/1‑148 (MQ=255)
                                                                   aaTGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGa                                                                                 >  2:3022657/1‑148 (MQ=255)
                                                                     tGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATg                                                                               >  2:2888217/1‑148 (MQ=255)
                                                                      ggATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGa                                                                              >  1:3297341/1‑148 (MQ=255)
                                                                      ggATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGa                                                                              <  2:1243443/148‑1 (MQ=255)
                                                                       gATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGaa                                                                             >  2:3418920/1‑148 (MQ=255)
                                                                       gATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGaa                                                                             >  1:2966689/1‑148 (MQ=255)
                                                                          tGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGaaaaa                                                                          <  2:2457516/148‑1 (MQ=255)
                                                                           gTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCCGAACGTCGATGAAAAAg                                                                         >  1:17338/1‑148 (MQ=255)
                                                                            tCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGc                                                                        <  1:1854711/148‑1 (MQ=255)
                                                                             cTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCg                                                                       >  1:158402/1‑148 (MQ=255)
                                                                                 ccGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTa                                                                   <  1:2283342/148‑1 (MQ=255)
                                                                                 ccGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTa                                                                   <  1:2195032/148‑1 (MQ=255)
                                                                                  cGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAATCGGCTaa                                                                  >  1:2827568/1‑148 (MQ=255)
                                                                                    gAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATc                                                                <  1:1464788/148‑1 (MQ=255)
                                                                                    gAAGTTATATCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGAAGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCAAGTTGAATGCTGAAGAGATGGCCCTCAAATCCCACAACGTCGATGAAAAAGCGGCTAAta                                                                >  2:3379533/1‑147 (MQ=255)
                                                                                      aGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCa                                                              >  2:2287355/1‑148 (MQ=255)
                                                                                      aGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCa                                                              <  1:2922675/148‑1 (MQ=255)
                                                                                        ttATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCCGAACGTCGTTGAAATAGCGGCTAATCCAGc                                                            >  1:3398302/1‑148 (MQ=255)
                                                                                         tATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCa                                                           <  2:1555689/148‑1 (MQ=255)
                                                                                         tATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCa                                                           <  2:56427/148‑1 (MQ=255)
                                                                                          aTCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCaa                                                          >  2:3482177/1‑148 (MQ=255)
                                                                                            ctcAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACg                                                        >  2:2640551/1‑148 (MQ=255)
                                                                                             tcAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGc                                                       <  2:721057/148‑1 (MQ=255)
                                                                                               aaaaGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCt                                                     <  1:2207922/148‑1 (MQ=255)
                                                                                                aaaGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTg                                                    >  2:1469282/1‑148 (MQ=255)
                                                                                                  aGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGat                                                  >  1:1385581/1‑148 (MQ=255)
                                                                                                  aGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGat                                                  >  1:2392345/1‑148 (MQ=255)
                                                                                                    aaCCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGatat                                                >  2:2842162/1‑148 (MQ=255)
                                                                                                    aaCCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGatat                                                <  2:2488657/148‑1 (MQ=255)
                                                                                                     aCCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATg                                               >  2:3473102/1‑148 (MQ=255)
                                                                                                          ttGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGcaca                                          >  1:2249653/1‑148 (MQ=255)
                                                                                                           tGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGGACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACAt                                         <  2:725534/148‑1 (MQ=255)
                                                                                                              tCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATctg                                      <  2:216753/148‑1 (MQ=255)
                                                                                                                  gTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTgg                                  <  2:2685925/148‑1 (MQ=255)
                                                                                                                  gTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTgg                                  <  2:18494/148‑1 (MQ=255)
                                                                                                                   tAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTggg                                 >  2:2583251/1‑148 (MQ=255)
                                                                                                                      ttCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTaa                              <  1:2005242/148‑1 (MQ=255)
                                                                                                                        cgtcgtGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaag                            <  2:3203324/148‑1 (MQ=255)
                                                                                                                        cgtcgtGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaag                            <  2:3300502/148‑1 (MQ=255)
                                                                                                                        cgtcgtGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaag                            <  1:2160476/148‑1 (MQ=255)
                                                                                                                          tcgtGGCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaagaa                          <  1:547059/148‑1 (MQ=255)
                                                                                                                          tcgtGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAaagaa                          <  2:3064383/148‑1 (MQ=255)
                                                                                                                             tGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGc                       >  2:2185469/1‑148 (MQ=255)
                                                                                                                              ggCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAAAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCg                      >  2:1642111/1‑148 (MQ=255)
                                                                                                                               gCGATGATGAATGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGa                     <  1:585867/148‑1 (MQ=255)
                                                                                                                                cGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGAt                    <  2:1435434/148‑1 (MQ=255)
                                                                                                                                 gatgatGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGAtt                   <  1:2951707/148‑1 (MQ=255)
                                                                                                                                  atgatgAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTa                  >  1:1920534/1‑148 (MQ=255)
                                                                                                                                       gaATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCa             <  1:696985/148‑1 (MQ=255)
                                                                                                                                         aTGGTTCTCTATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAg           >  1:1106954/1‑148 (MQ=255)
                                                                                                                                         aTGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAg           >  1:818375/1‑148 (MQ=255)
                                                                                                                                         aTGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAg           <  1:1215686/148‑1 (MQ=255)
                                                                                                                                             ttCTCGATTGTACGCTTGACTCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATc       <  1:3344687/148‑1 (MQ=255)
                                                                                                                                              tctcGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCt      >  1:2862218/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCGATACGCTGGCCTCGGATCAATCACAACTGTATGCCCTGCGCATCAAATTAAGCAACCCAGCTGAATGGATTGTCTTACCGGAAGTTATCTCAAAAGAACCGCTTGGTCCGGTAGTTCGTCGTGGCGATGATGAATGGTTCTCGATTGTACGCTGGACGCTTTTCGCCATGCTGAATGCTGAAGAGATGGGCATCAATTCCCAGAACGTCGATGAAAAAGCGGCTAATCCAGCAACGCCTGATATGGCACATCTGCTGGGTAAAGAAGGCGATTACGGCAAGGATCTGAAG  >  CP009273/3412973‑3413266

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: