Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,495,918 C→A 10.2% V683F (GTT→TTT)  fryA ← putative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,495,9180CA10.2% 921.0 / 52.4 261V683F (GTT→TTT) fryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (120/114);  new base A (16/10);  total (137/124)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 264 total aligned reads displayed.

TTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATAC  >  CP009273/2495769‑2496067
                                                                                                                                                     |                                                                                                                                                     
ttCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGTATGTTCAGTGCGCCCGTTAACGACGAGGTTGCCGCACAGGAACTGGATAGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCAATCGCCAGCAGAGGGCGAAc                                                                                                                                                       >  1:2440145/1‑150 (MQ=255)
 tCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACg                                                                                                                                                      >  2:784599/1‑150 (MQ=255)
    cGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtct                                                                                                                                                   <  1:2085928/150‑1 (MQ=255)
       tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                                                                                >  2:2785250/1‑150 (MQ=255)
       tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                                                                                >  1:4077126/1‑150 (MQ=255)
            gACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCgg                                                                                                                                           <  2:4401064/150‑1 (MQ=255)
              cATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCg                                                                                                                                         <  2:1674300/150‑1 (MQ=255)
              cATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATAGTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCg                                                                                                                                         >  1:3672632/1‑150 (MQ=255)
                 tttcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                                                                                      >  1:3111475/2‑150 (MQ=255)
                 cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                                                                                      <  2:2462490/150‑1 (MQ=255)
                      tcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCg                                                                                                                                 <  1:1871790/150‑1 (MQ=255)
                        cAGCTCAATCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGt                                                                                                                               >  2:4332204/1‑150 (MQ=255)
                        cAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGt                                                                                                                               >  2:4332255/1‑150 (MQ=255)
                         aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              <  1:3072895/150‑1 (MQ=255)
                         aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              <  2:3246180/150‑1 (MQ=255)
                         aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              <  2:3246165/150‑1 (MQ=255)
                                aCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACg                                                                                                                       >  1:932813/1‑150 (MQ=255)
                                    aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCATTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTc                                                                                                                   <  2:3108291/150‑1 (MQ=255)
                                    aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTc                                                                                                                   <  2:3107954/150‑1 (MQ=255)
                                      gTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCaa                                                                                                                 >  2:2243631/1‑150 (MQ=255)
                                        tCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAtt                                                                                                               >  2:3392250/1‑150 (MQ=255)
                                         cAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAttt                                                                                                              >  2:2571781/1‑150 (MQ=255)
                                             gcgcCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTg                                                                                                          <  2:3888160/150‑1 (MQ=255)
                                               gcCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTGCAATTTCCTgcg                                                                                                        >  2:3729405/1‑150 (MQ=255)
                                                cccGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTgcgc                                                                                                       >  2:4163844/1‑150 (MQ=255)
                                                 ccGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                                                                                      >  2:968602/1‑150 (MQ=255)
                                                   gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACt                                                                                                    <  2:111600/150‑1 (MQ=255)
                                                   gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACt                                                                                                    >  1:4302479/1‑150 (MQ=255)
                                                   gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACt                                                                                                    <  2:2564391/150‑1 (MQ=255)
                                                    ttAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTg                                                                                                   <  1:3880321/150‑1 (MQ=255)
                                                     tAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGc                                                                                                  >  1:2311967/1‑150 (MQ=255)
                                                      aaCGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCg                                                                                                 <  2:576526/150‑1 (MQ=255)
                                                          ccGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAt                                                                                             <  1:1668367/150‑1 (MQ=255)
                                                           cGAGGTTGCCGGACAGGAACTGGATCTCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAtt                                                                                            >  2:3667267/1‑150 (MQ=255)
                                                           cGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAtt                                                                                            >  1:3749204/1‑150 (MQ=255)
                                                                ttGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACAt                                                                                       <  2:1190361/150‑1 (MQ=255)
                                                                  gCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGc                                                                                     >  2:70240/1‑150 (MQ=255)
                                                                    cGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCt                                                                                   >  1:2525889/1‑150 (MQ=255)
                                                                      cacaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                                 <  2:2063639/150‑1 (MQ=255)
                                                                      cacaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                                 <  2:2882142/150‑1 (MQ=255)
                                                                      cacaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                                 >  1:633313/1‑150 (MQ=255)
                                                                      cacaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                                 <  2:2065012/150‑1 (MQ=255)
                                                                       acaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAc                                                                                <  2:4493305/150‑1 (MQ=255)
                                                                        caGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                                                                               <  1:2275452/150‑1 (MQ=255)
                                                                        caGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                                                                               <  1:2275604/150‑1 (MQ=255)
                                                                            aaCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGcc                                                                           >  2:26123/1‑150 (MQ=255)
                                                                              cTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAg                                                                         <  1:3012796/150‑1 (MQ=255)
                                                                               tGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGt                                                                        >  2:1600230/1‑150 (MQ=255)
                                                                               tGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGt                                                                        >  2:1599506/1‑150 (MQ=255)
                                                                                 gATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTc                                                                      >  1:2132291/1‑150 (MQ=255)
                                                                                  aTCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTcc                                                                     >  1:3825718/1‑150 (MQ=255)
                                                                                   tCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTccc                                                                    >  1:4363575/1‑150 (MQ=255)
                                                                                    cGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCg                                                                   >  1:1059883/1‑150 (MQ=255)
                                                                                      ccTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGac                                                                 >  2:1291267/1‑150 (MQ=255)
                                                                                           tctTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGcc                                                            >  1:328433/1‑150 (MQ=255)
                                                                                            ctTTATTGCTAAAATCTTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                           >  2:2950765/1‑150 (MQ=255)
                                                                                            ctTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                           <  2:2625897/150‑1 (MQ=255)
                                                                                            ctTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                           <  1:4030935/150‑1 (MQ=255)
                                                                                             tttATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTc                                                          <  1:1225772/150‑1 (MQ=255)
                                                                                               tATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGc                                                        >  2:1419455/1‑150 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       <  1:2819273/150‑1 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       <  1:1132379/150‑1 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       <  1:2922026/150‑1 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       <  2:1511237/150‑1 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       >  1:1425365/1‑150 (MQ=255)
                                                                                                aTTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCt                                                       >  2:2984799/1‑150 (MQ=255)
                                                                                                  tGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTAt                                                     >  1:4153998/1‑150 (MQ=255)
                                                                                                   gCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATc                                                    <  1:2684786/150‑1 (MQ=255)
                                                                                                   gCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATc                                                    <  1:2684842/150‑1 (MQ=255)
                                                                                                   gCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATc                                                    <  1:2684764/150‑1 (MQ=255)
                                                                                                   gCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATc                                                    <  1:665666/150‑1 (MQ=255)
                                                                                                      aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTGAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAg                                                 >  2:812660/1‑150 (MQ=255)
                                                                                                      aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAg                                                 >  2:1852683/1‑150 (MQ=255)
                                                                                                      aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAg                                                 >  2:811038/1‑150 (MQ=255)
                                                                                                      aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGAGTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAg                                                 >  2:4308688/1‑150 (MQ=255)
                                                                                                        aaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGtt                                               >  1:3165819/1‑150 (MQ=255)
                                                                                                           ccTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGAc                                            <  1:817993/150‑1 (MQ=255)
                                                                                                             tGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGa                                          >  1:127943/1‑150 (MQ=255)
                                                                                                                tCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGc                                       >  2:1699638/1‑150 (MQ=255)
                                                                                                                 caacaaAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCt                                      <  1:3667267/150‑1 (MQ=255)
                                                                                                                 caacaaAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCt                                      <  1:1402265/150‑1 (MQ=255)
                                                                                                                  aacaaAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTg                                     >  2:3564840/1‑150 (MQ=255)
                                                                                                                  aacaaAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTg                                     >  2:587407/1‑150 (MQ=255)
                                                                                                                  aacaaAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTg                                     <  2:2311967/150‑1 (MQ=255)
                                                                                                                       acatagTCTCCAGCGCCAGCAGTGGGCGACCGTCGTCTTCCGGCGTAACGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCttt                                <  1:4443421/148‑1 (MQ=255)
                                                                                                                         atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTc                              <  1:3025/150‑1 (MQ=255)
                                                                                                                         atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTc                              <  1:3031/150‑1 (MQ=255)
                                                                                                                         atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTc                              >  2:3007566/1‑150 (MQ=255)
                                                                                                                              ctcCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAccgc                         >  2:1747934/1‑150 (MQ=255)
                                                                                                                                    cgcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAAt                   >  2:2205170/1‑150 (MQ=255)
                                                                                                                                    cgcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAAt                   >  2:2427165/1‑150 (MQ=255)
                                                                                                                                      ccAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAc                 <  1:3965168/150‑1 (MQ=255)
                                                                                                                                      ccAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAc                 >  2:1816896/1‑150 (MQ=255)
                                                                                                                                         gcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgcg              <  2:4077368/150‑1 (MQ=255)
                                                                                                                                         gcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgcg              <  1:3564647/150‑1 (MQ=255)
                                                                                                                                              gggCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCta         <  1:1011793/150‑1 (MQ=255)
                                                                                                                                              gggCGAACGTCTTCGTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCta         <  1:4122262/150‑1 (MQ=255)
                                                                                                                                                gCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtact       <  2:3749204/150‑1 (MQ=255)
                                                                                                                                                   aaCGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcat    <  1:3439357/150‑1 (MQ=255)
                                                                                                                                                     cGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatac  >  2:3595622/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATAC  >  CP009273/2495769‑2496067

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: