Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,495,918 C→A 32.6% 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,9180CA32.6% 582.6 / 305.5 328V683F (GTT→TTT) fryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (120/101);  new base A (62/45);  total (182/146)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

TCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATAC  >  CP009273/2495770‑2496067
                                                                                                                                                    |                                                                                                                                                     
tCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACg                                                                                                                                                      <  1:1474288/150‑1 (MQ=255)
tCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACg                                                                                                                                                      <  1:2606532/150‑1 (MQ=255)
 cccGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGt                                                                                                                                                     >  1:2585876/1‑150 (MQ=255)
   cGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGtct                                                                                                                                                   <  2:767011/150‑1 (MQ=255)
    gCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtctt                                                                                                                                                  <  2:4215088/150‑1 (MQ=255)
      tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                                                                                >  2:2785817/1‑150 (MQ=255)
      tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                                                                                >  2:2888128/1‑150 (MQ=255)
         cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTcc                                                                                                                                             <  1:4409777/150‑1 (MQ=255)
         cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTcc                                                                                                                                             >  2:2244703/1‑150 (MQ=255)
         cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTcc                                                                                                                                             >  2:2926577/1‑150 (MQ=255)
          aGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCg                                                                                                                                            >  2:3960286/1‑150 (MQ=255)
          aGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCg                                                                                                                                            >  1:1896188/1‑150 (MQ=255)
          aGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAAATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCg                                                                                                                                            <  2:3243880/150‑1 (MQ=255)
            aCATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGc                                                                                                                                          >  1:4098370/1‑150 (MQ=255)
               tcttctTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTa                                                                                                                                       >  1:4752629/1‑150 (MQ=255)
                cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCTCAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                                                                                      <  2:3500001/150‑1 (MQ=255)
                 ttcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGACCGTCTTCTTCCGGCGTaaa                                                                                                                                     >  1:522500/1‑150 (MQ=255)
                  tcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAg                                                                                                                                    <  1:1006461/150‑1 (MQ=255)
                        aGCTCAAACGGTTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              <  1:3249977/150‑1 (MQ=255)
                        aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              >  2:1106327/1‑150 (MQ=255)
                        aGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                              >  2:501508/1‑150 (MQ=255)
                         gCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGa                                                                                                                             <  1:774623/150‑1 (MQ=255)
                           tCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGt                                                                                                                           >  1:3867123/1‑150 (MQ=255)
                           tCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGt                                                                                                                           <  1:3930272/150‑1 (MQ=255)
                            cAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTa                                                                                                                          <  1:3495560/150‑1 (MQ=255)
                            cAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTa                                                                                                                          <  1:3495619/150‑1 (MQ=255)
                             aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                                                                                                                         >  2:4422844/1‑150 (MQ=255)
                             aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                                                                                                                         <  1:3680313/150‑1 (MQ=255)
                             aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                                                                                                                         >  2:4620338/1‑150 (MQ=255)
                                cGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGACAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGc                                                                                                                      >  2:4276414/1‑150 (MQ=255)
                                cGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGc                                                                                                                      >  2:1209788/1‑150 (MQ=255)
                                   aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTc                                                                                                                   <  2:3884972/150‑1 (MQ=255)
                                    tGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAACAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCa                                                                                                                  <  2:2333152/150‑1 (MQ=255)
                                      ttCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAt                                                                                                                <  2:177521/150‑1 (MQ=255)
                                      ttCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAt                                                                                                                >  2:21129/1‑150 (MQ=255)
                                      ttCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAAt                                                                                                                >  2:2947439/1‑150 (MQ=255)
                                           tGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCt                                                                                                           >  2:2616281/1‑150 (MQ=255)
                                           tGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCt                                                                                                           >  1:1128923/1‑150 (MQ=255)
                                                   ttAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTg                                                                                                   >  1:2023004/1‑150 (MQ=255)
                                                   ttAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTg                                                                                                   >  2:3580989/1‑150 (MQ=255)
                                                    tAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGc                                                                                                  >  1:4335922/1‑150 (MQ=255)
                                                    tAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATGTCCTGCGCACTGc                                                                                                  >  2:2846458/1‑150 (MQ=255)
                                                     aaCGCCGAGGTTGCCGCACTGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCTGTGAGTAACGCTTCAATTTCCTGCGCACTGct                                                                                                 <  2:2182624/150‑2 (MQ=255)
                                                        gCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCa                                                                                              <  2:3836337/150‑1 (MQ=255)
                                                        gCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCa                                                                                              <  1:3179013/150‑1 (MQ=255)
                                                           gAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTc                                                                                           >  1:2982860/1‑150 (MQ=255)
                                                             ggTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcac                                                                                         <  1:1566280/150‑1 (MQ=255)
                                                             ggTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTATTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcac                                                                                         <  2:657307/150‑1 (MQ=255)
                                                              gTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcaca                                                                                        >  2:1228321/1‑150 (MQ=255)
                                                               ttGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACAt                                                                                       <  1:4143476/150‑1 (MQ=255)
                                                                     cacaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                                 <  1:2703440/150‑1 (MQ=255)
                                                                      acaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAc                                                                                >  2:4758821/1‑150 (MQ=255)
                                                                      acaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGTGGCATTCACATGCCTGAc                                                                                <  1:36336/150‑1 (MQ=255)
                                                                        aGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcgc                                                                              >  1:59855/1‑150 (MQ=255)
                                                                           aaCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGcc                                                                           >  1:680801/1‑150 (MQ=255)
                                                                                 aTCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTcc                                                                     >  1:2825696/1‑150 (MQ=255)
                                                                                 aTCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATGTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTcc                                                                     <  2:4296005/150‑1 (MQ=255)
                                                                                  tCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTccc                                                                    >  1:2308403/1‑150 (MQ=255)
                                                                                   cGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCg                                                                   <  1:4516112/150‑1 (MQ=255)
                                                                                   cGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCg                                                                   <  2:1438899/150‑1 (MQ=255)
                                                                                     ccTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGac                                                                 >  1:4283652/1‑150 (MQ=255)
                                                                                           ctTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                           <  1:3388222/150‑1 (MQ=255)
                                                                                           ctTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                           >  1:1335883/1‑150 (MQ=255)
                                                                                             ttATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCg                                                         <  1:944600/150‑1 (MQ=255)
                                                                                                  gCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATc                                                    <  1:1960364/150‑1 (MQ=255)
                                                                                                    tAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCa                                                  >  1:902359/1‑150 (MQ=255)
                                                                                                    tAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCa                                                  >  1:4234603/1‑150 (MQ=255)
                                                                                                     aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAg                                                 >  2:4030155/1‑150 (MQ=255)
                                                                                                     aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGGAAAGGCGGTGAGTAACGCTTCAATGTCCTGCGCACTGCGGCATGCACATGCCTGAGGCGCCAGTTCCCGACACGCCTCGCTCTcccg                                                 >  2:1083601/1‑148 (MQ=255)
                                                                                                      aaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGt                                                >  2:2452946/1‑150 (MQ=255)
                                                                                                       aaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGtt                                               <  2:680801/150‑1 (MQ=255)
                                                                                                        aTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTg                                              >  2:1661366/1‑150 (MQ=255)
                                                                                                        aTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTg                                              >  1:3471062/1‑150 (MQ=255)
                                                                                                          ccTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGAc                                            <  2:1959434/150‑1 (MQ=255)
                                                                                                           cTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACg                                           <  1:2926577/150‑1 (MQ=255)
                                                                                                                  acaaAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTgg                                    >  2:2547099/1‑150 (MQ=255)
                                                                                                                   caaAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTCTCCAGTTGACGAAGCTGGc                                   >  2:4684995/1‑150 (MQ=255)
                                                                                                                   caaAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGc                                   >  2:815887/1‑150 (MQ=255)
                                                                                                                    aaaGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCt                                  >  1:1623251/1‑150 (MQ=255)
                                                                                                                     aaGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACTCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCtt                                 <  1:858984/150‑1 (MQ=255)
                                                                                                                     aaGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCtt                                 >  1:3095414/1‑150 (MQ=255)
                                                                                                                       gATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCtttt                               <  2:3867123/150‑1 (MQ=255)
                                                                                                                        atatTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTc                              >  1:563715/1‑150 (MQ=255)
                                                                                                                          atTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAc                            >  1:3702476/1‑150 (MQ=255)
                                                                                                                           ttCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAcc                           >  2:1861966/1‑150 (MQ=255)
                                                                                                                             ctcCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAccgc                         >  2:3159947/1‑150 (MQ=255)
                                                                                                                               ccAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAccgccg                       >  1:3897547/1‑150 (MQ=255)
                                                                                                                                   cgcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAAt                   <  1:354036/150‑1 (MQ=255)
                                                                                                                                    gcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATa                  <  1:3739305/150‑1 (MQ=255)
                                                                                                                                    gcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATa                  >  2:613069/1‑150 (MQ=255)
                                                                                                                                    gcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATa                  <  1:1537135/150‑1 (MQ=255)
                                                                                                                                    gcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATa                  >  2:880755/1‑150 (MQ=255)
                                                                                                                                     ccAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAc                 <  2:4558277/150‑1 (MQ=255)
                                                                                                                                       agcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc               >  2:1290933/1‑150 (MQ=255)
                                                                                                                                       agcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc               >  2:1642394/1‑150 (MQ=255)
                                                                                                                                        gcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgcg              <  1:2785817/150‑1 (MQ=255)
                                                                                                                                               gCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtact       <  2:3741444/150‑1 (MQ=255)
                                                                                                                                                cGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactc      >  1:2426128/1‑150 (MQ=255)
                                                                                                                                                 gAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactca     >  2:4353137/1‑150 (MQ=255)
                                                                                                                                                    cGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatac  >  2:201547/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATAC  >  CP009273/2495770‑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: