Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,495,516 C→A 10.4% S297* (TCG→TAG)  frlB → fructoselysine‑6‑P‑deglycase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,495,5160CA10.4% 1262.2 / 76.0 354S297* (TCG→TAG) frlBfructoselysine‑6‑P‑deglycase
Reads supporting (aligned to +/- strand):  ref base C (165/151);  new base A (23/15);  total (188/166)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 358 total aligned reads displayed.

CCGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTCCCGGCCCG  >  CP009273/3495367‑3495665
                                                                                                                                                     |                                                                                                                                                     
 cGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCg                                                                                                                                                      >  1:1493500/1‑150 (MQ=255)
 cGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCg                                                                                                                                                      >  2:3993781/1‑150 (MQ=255)
  gCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGc                                                                                                                                                     >  1:4474980/1‑150 (MQ=255)
    cATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCaa                                                                                                                                                   >  2:4521317/1‑150 (MQ=255)
     aTGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAg                                                                                                                                                  >  2:5102338/1‑150 (MQ=255)
        gCCCGCGGGAGATTGTCGAACCGGGCGGTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGc                                                                                                                                               <  2:427710/150‑1 (MQ=255)
         cccGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCt                                                                                                                                              <  2:1241800/150‑1 (MQ=255)
           cGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGc                                                                                                                                            >  1:2922922/1‑150 (MQ=255)
            gCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCa                                                                                                                                           <  2:5098898/150‑1 (MQ=255)
            gCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCa                                                                                                                                           <  1:5136330/150‑1 (MQ=255)
              tGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCAcc                                                                                                                                         >  1:1839669/1‑150 (MQ=255)
              tGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCAcc                                                                                                                                         >  2:4934389/1‑150 (MQ=255)
                 agaTTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGt                                                                                                                                      <  2:2247981/150‑1 (MQ=255)
                  gaTTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACCGCGTACTGACAACGTGAGCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGtg                                                                                                                                     >  2:704662/1‑150 (MQ=255)
                   aTTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGtgg                                                                                                                                    <  1:1910300/150‑1 (MQ=255)
                     tGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCCAAGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGAACTGACAACGTGATCGTCATCGATTAAGCCGAAATTTCGCAAGGGCTGCACCCGtggct                                                                                                                                  >  1:4998163/1‑150 (MQ=255)
                     tGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGtggct                                                                                                                                  >  2:237735/1‑150 (MQ=255)
                     tGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGtggct                                                                                                                                  >  1:1195811/1‑150 (MQ=255)
                        cGAACCGGGCGTTCCGTTCCTGGTCCTGCTCGGCACTGATGAAAGTCGCCACACCACCGAACGCGCCATTAACGGTGGTAAACCGCGTACTGACACCGTGATCGTCATCGCTTCCGCCGAACGTTCGCAAGGGCGGCACCCGtggctggc                                                                                                                               >  2:2652196/1‑150 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              <  1:1265373/150‑1 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              >  2:16695/1‑150 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              >  2:604502/1‑150 (MQ=255)
                         gAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCa                                                                                                                              >  2:3258777/1‑150 (MQ=255)
                            ccGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCg                                                                                                                           >  2:3652627/1‑150 (MQ=255)
                              gggCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGtt                                                                                                                         <  2:2530979/150‑1 (MQ=255)
                               ggCGTTCCGTTCCTGTTCCAGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTc                                                                                                                        >  2:4130698/1‑150 (MQ=255)
                                gCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTCAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTcc                                                                                                                       <  1:3916192/150‑1 (MQ=255)
                                    tccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGAt                                                                                                                   <  2:1008916/150‑1 (MQ=255)
                                    tccgttccTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGAt                                                                                                                   >  2:2052148/1‑150 (MQ=255)
                                       gttcctgttcctGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGtt                                                                                                                <  2:1554986/150‑1 (MQ=255)
                                         tcctgttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCg                                                                                                              <  1:4958935/150‑1 (MQ=255)
                                           ctgttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGTCACCGTTCCTGATGTTCGTg                                                                                                            <  2:4103179/150‑1 (MQ=255)
                                              ttcctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCa                                                                                                         <  2:2669472/150‑1 (MQ=255)
                                                cctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAAt                                                                                                       <  2:3761682/150‑1 (MQ=255)
                                                 ctgCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATg                                                                                                      >  2:2084713/1‑150 (MQ=255)
                                                    cTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAg                                                                                                   <  1:3339154/150‑1 (MQ=255)
                                                       ggCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACACCGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTgg                                                                                                <  2:2769899/150‑1 (MQ=255)
                                                       ggCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTgg                                                                                                >  1:4742585/1‑150 (MQ=255)
                                                       ggCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTgg                                                                                                >  2:4311718/1‑150 (MQ=255)
                                                       ggCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTgg                                                                                                >  2:3214386/1‑150 (MQ=255)
                                                         cAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGct                                                                                              >  1:1724614/1‑150 (MQ=255)
                                                           atgatgAAAGTCGACTCACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGAGTACGCCGAAATTTCGCAAGGTCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCTATGGAGTGGctct                                                                                            <  2:423632/150‑1 (MQ=255)
                                                            tgatgaAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTg                                                                                           <  2:2862898/150‑1 (MQ=255)
                                                             gatgaAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGc                                                                                          >  2:3852686/1‑150 (MQ=255)
                                                             gatgaAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGc                                                                                          >  1:4897275/1‑150 (MQ=255)
                                                                  aaGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGctacta                                                                                     <  2:4711662/150‑1 (MQ=255)
                                                                  aaGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGctacta                                                                                     >  2:1343488/1‑150 (MQ=255)
                                                                      cGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTg                                                                                 <  1:3689763/150‑1 (MQ=255)
                                                                         cacaCCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCt                                                                              <  2:4036103/150‑1 (MQ=255)
                                                                          acacCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTa                                                                             <  2:1141172/150‑1 (MQ=255)
                                                                          acacCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTa                                                                             >  2:2002582/1‑150 (MQ=255)
                                                                              caccGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTa                                                                         >  2:3419144/1‑150 (MQ=255)
                                                                                  gAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCACGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACaaa                                                                     <  2:1778711/150‑1 (MQ=255)
                                                                                  gAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACaaa                                                                     >  1:572258/1‑150 (MQ=255)
                                                                                       cgcCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATca                                                                <  2:3496417/150‑1 (MQ=255)
                                                                                        gcCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATcac                                                               <  1:2268825/150‑1 (MQ=255)
                                                                                         ccATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATcaca                                                              >  1:4484911/1‑150 (MQ=255)
                                                                                          cATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACaa                                                             >  2:3517413/1‑150 (MQ=255)
                                                                                          cATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACaa                                                             <  2:1478839/150‑1 (MQ=255)
                                                                                          cATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACaa                                                             >  1:2759563/1‑150 (MQ=255)
                                                                                            ttAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAcc                                                           <  2:4077148/150‑1 (MQ=255)
                                                                                            ttAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAcc                                                           >  1:3280217/1‑150 (MQ=255)
                                                                                            ttAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAcc                                                           >  1:3043805/1‑150 (MQ=255)
                                                                                             tAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAAccc                                                          >  1:1424620/1‑150 (MQ=255)
                                                                                              aaCTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCg                                                         <  2:2599380/150‑1 (MQ=255)
                                                                                              aaCTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCg                                                         >  1:4723831/1‑150 (MQ=255)
                                                                                              aaCTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCg                                                         <  1:923779/150‑1 (MQ=255)
                                                                                                 tttGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGAt                                                      <  1:4920135/150‑1 (MQ=255)
                                                                                                  ttGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATg                                                     >  1:4023988/1‑150 (MQ=255)
                                                                                                    gTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGaa                                                   <  1:721482/150‑1 (MQ=255)
                                                                                                      tAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcg                                                 >  1:674610/1‑150 (MQ=255)
                                                                                                       aaaCAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgc                                                <  1:237735/150‑1 (MQ=255)
                                                                                                        aaCAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgcc                                               >  2:295324/1‑150 (MQ=255)
                                                                                                         aCAGTGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAAcgccg                                              >  1:2567101/1‑150 (MQ=255)
                                                                                                           aGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCt                                            <  2:4243230/150‑1 (MQ=255)
                                                                                                             cGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTAt                                          <  1:1004511/150‑1 (MQ=255)
                                                                                                              gTACTGACACCGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTTGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTAtt                                         <  2:2291669/150‑1 (MQ=255)
                                                                                                                aCTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTAc                                       >  2:3840840/1‑150 (MQ=255)
                                                                                                                aCTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTAc                                       >  1:4004027/1‑150 (MQ=255)
                                                                                                                    aCAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACggtg                                   <  2:1061457/150‑1 (MQ=255)
                                                                                                                      aaCGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACggtggt                                 <  1:4067237/150‑1 (MQ=255)
                                                                                                                         gTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtg                              >  1:2441129/1‑150 (MQ=255)
                                                                                                                          tGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtgg                             >  1:2051349/1‑150 (MQ=255)
                                                                                                                          tGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtgg                             >  1:2206812/1‑150 (MQ=255)
                                                                                                                             tCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtggtgg                          >  2:905003/1‑150 (MQ=255)
                                                                                                                             tCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCtggtgg                          <  1:3090053/150‑1 (MQ=255)
                                                                                                                                 cATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAata                      >  2:1937079/1‑150 (MQ=255)
                                                                                                                                   tCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGACTAtt                    >  2:1196349/1‑150 (MQ=255)
                                                                                                                                   tCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATAtt                    >  1:4331702/1‑150 (MQ=255)
                                                                                                                                    cGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTa                   >  2:2980381/1‑150 (MQ=255)
                                                                                                                                     gATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTaa                  >  2:1488529/1‑150 (MQ=255)
                                                                                                                                      aTTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAAt                 >  2:4404010/1‑150 (MQ=255)
                                                                                                                                      aTTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAAt                 >  1:213127/1‑150 (MQ=255)
                                                                                                                                       ttACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATc                >  1:3526286/1‑150 (MQ=255)
                                                                                                                                       ttACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATc                >  1:4160326/1‑150 (MQ=255)
                                                                                                                                           gCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCctc            <  1:4668099/150‑1 (MQ=255)
                                                                                                                                              gAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcc         <  2:4675412/150‑1 (MQ=255)
                                                                                                                                                aaTTTAGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccg       <  2:3058072/150‑1 (MQ=255)
                                                                                                                                                 aTTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccgg      >  1:1336884/1‑150 (MQ=255)
                                                                                                                                                  tttCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTcccggc     <  2:480378/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCCATGGCCCGCTGGAGATTGTCGAACCGGGCGTTCCGTTCCTGTTCCTGCTCGGCAATGATGAAAGTCGCCACACCACCGAACGCGCCATTAACTTTGTTAAACAGCGTACTGACAACGTGATCGTCATCGATTACGCCGAAATTTCGCAAGGGCTGCACCCGTGGCTGGCACCGTTCCTGATGTTCGTGCCAATGGAGTGGCTCTGCTACTACCTGTCTATTTACAAAGATCACAACCCGGATGAACGCCGCTATTACGGTGGTCTGGTGGAATATTAATCCCTCTCCCGGCCCG  >  CP009273/3495367‑3495665

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: