Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,756,668 C→G intergenic (+550/+90) RPPX_16775 → / ← RPPX_16780 hypothetical protein/phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,756,6680CG76.2% 783.2 / 266.1 399intergenic (+550/+90)RPPX_16775/RPPX_16780hypothetical protein/phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (39/56);  new base G (158/146);  total (197/202)
Fisher's exact test for biased strand distribution p-value = 7.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 407 total aligned reads displayed.

CAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGAATGCCAACTCCAGCGACTG  >  CP009974/3756518‑3756814
                                                                                                                                                      |                                                                                                                                                  
  aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGTCCGTCATCGCCCGCCGCAAGgcc                                                                                                                                                  >  1:1003216/1‑151 (MQ=17)
  aaaaCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGgcc                                                                                                                                                  >  1:4903876/1‑151 (MQ=255)
      ccGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTgcg                                                                                                                                              <  2:3829693/151‑1 (MQ=255)
      ccGTCGGCTCCCATCAGCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTgcg                                                                                                                                              <  2:5238397/151‑1 (MQ=255)
        gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCt                                                                                                                                            <  1:1809755/151‑1 (MQ=255)
        gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCt                                                                                                                                            <  2:3339967/151‑1 (MQ=255)
        gTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCt                                                                                                                                            <  1:5764925/151‑1 (MQ=255)
         tCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTc                                                                                                                                           >  1:3636861/1‑151 (MQ=255)
          cGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCg                                                                                                                                          >  2:2013690/1‑151 (MQ=255)
           ggCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCgg                                                                                                                                         >  2:236621/1‑151 (MQ=255)
              tCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCCCTGCACGCCCCATCCACGGGGCCCCACCGCTTGACGGGCATCCATACCCGTCATCTCCCTCCGCAAGGCCGCCGCTAGTcct                                                                                                                                      >  1:3521357/1‑151 (MQ=14)
                 catcatCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGcctcct                                                                                                                                   >  1:1913838/1‑151 (MQ=255)
                   tcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGa                                                                                                                                 >  2:4916919/1‑151 (MQ=255)
                   tcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCGCCCTTCGCAAGCCCTGCGCTCGGcctgctcc                                                                                                                                 >  1:343525/1‑149 (MQ=255)
                     atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCCCGCGCTGCAGGCCCCATCCATGGGGCCTCCGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCCAGCCCTGCGCCCGGCCTCCTGaaa                                                                                                                               >  2:2532926/1‑150 (MQ=14)
                     atcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGacg                                                                                                                               >  2:5870377/1‑149 (MQ=255)
                      tcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGt                                                                                                                              <  2:758127/151‑1 (MQ=255)
                         ggCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGc                                                                                                                           <  2:448900/151‑1 (MQ=255)
                          gCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCa                                                                                                                          <  1:1411298/151‑1 (MQ=255)
                          gCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCa                                                                                                                          >  2:5287287/1‑151 (MQ=255)
                              gtACGCCGCGCTCCGGGGTCCTCTCGCTCCGGCCTCGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTAGGCCTCCTGACGTCGCAATCt                                                                                                                      <  2:3130393/150‑1 (MQ=14)
                               tACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTg                                                                                                                     <  1:3867699/151‑1 (MQ=255)
                               tACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCCCCTGAAGTCGGaagcac                                                                                                                     >  1:2913869/1‑147 (MQ=255)
                                aCGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATGCATGCCCGTCATCTCCCTCCGCACGGCCTGCGCTCGCCCTCCTGAAGTCGCGATCCGc                                                                                                                    >  2:1626684/1‑151 (MQ=14)
                                     gcgcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTGCTGAAGTCGCAATCTGCGTCGc                                                                                                               >  2:3964222/1‑151 (MQ=255)
                                       gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                                                                             <  2:5689378/151‑1 (MQ=255)
                                       gcTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                                                                             >  1:6444856/1‑151 (MQ=255)
                                        cTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCCTGGGGCCTCAGCGCTTGACGGGCGTCACTGCCCGCCCTCTCCCTCCGCAGGGCCTGCGCGCGGCCTCCCGAAGTCGCACTCTGCGTCGccgg                                                                                                            >  2:2742883/1‑149 (MQ=14)
                                        cTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTg                                                                                                            <  1:2553882/151‑1 (MQ=255)
                                           cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAAc                                                                                                         >  1:809486/1‑151 (MQ=255)
                                           cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAAc                                                                                                         <  1:55137/151‑1 (MQ=255)
                                              ggTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGCAGTCGCAAGCTGCGCCGCCTGAACCAc                                                                                                      >  1:1536318/1‑151 (MQ=11)
                                               gTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGACCAAcc                                                                                                     >  1:52339/1‑151 (MQ=255)
                                                 ccccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                                                                   <  1:973320/151‑1 (MQ=255)
                                                 ccccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                                                                   <  1:759808/151‑1 (MQ=255)
                                                  cccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcg                                                                                                  <  1:4744330/151‑1 (MQ=255)
                                                   ccTCGCTCCGGCCTTGCCCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                                                                                 <  1:3073101/151‑1 (MQ=255)
                                                    cTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgcg                                                                                                <  1:4245688/151‑1 (MQ=255)
                                                     tCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCACGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCgcgcgc                                                                                               >  2:841200/1‑151 (MQ=255)
                                                       gCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCtt                                                                                             >  1:2030132/1‑151 (MQ=255)
                                                         tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCGATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAg                                                                                           <  2:2283996/151‑1 (MQ=255)
                                                         tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAg                                                                                           <  1:2207371/151‑1 (MQ=255)
                                                           cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGaa                                                                                         <  1:6507799/151‑1 (MQ=255)
                                                           cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGaa                                                                                         >  2:4170416/1‑151 (MQ=255)
                                                           cGGCCTTGCTCCCGGGAGCACCCGCTGCAGGCCCCATCCATGGAGCCTCAGCGCTTGACGGGCATACATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGaa                                                                                         <  2:1918110/151‑1 (MQ=11)
                                                              ccTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCa                                                                                      >  1:510291/1‑151 (MQ=255)
                                                               cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCaa                                                                                     >  2:3195409/1‑151 (MQ=255)
                                                                ttGCTCCCGGGAGCACGCGCTGCGGGCCCCGTCCAGGGGGCCTCAGCGCTGGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAg                                                                                    <  2:473439/151‑1 (MQ=14)
                                                                ttGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAg                                                                                    >  1:1470868/1‑151 (MQ=255)
                                                                ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTCGCAGCAAg                                                                                    >  2:6196538/1‑151 (MQ=11)
                                                                ttGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAg                                                                                           <  1:5343575/144‑1 (MQ=255)
                                                                 tGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGa                                                                                   <  2:3494896/151‑1 (MQ=255)
                                                                 tGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGa                                                                                   <  1:6167982/151‑1 (MQ=255)
                                                                   cgcccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATc                                                                                 <  1:5625630/149‑1 (MQ=255)
                                                                      ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaa                                                                              <  1:1179924/151‑1 (MQ=255)
                                                                        gggAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaag                                                                            <  1:3700852/151‑1 (MQ=255)
                                                                        gggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaag                                                                            >  1:4476523/1‑151 (MQ=255)
                                                                         ggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagc                                                                           <  1:4658019/151‑1 (MQ=255)
                                                                         ggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagc                                                                           <  1:1736397/151‑1 (MQ=255)
                                                                           aGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcga                                                                         >  1:594348/1‑151 (MQ=255)
                                                                           aGCACGCGATACAGGCCCCATCCACGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagcga                                                                         <  1:4975436/151‑1 (MQ=17)
                                                                             cACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCCTAGAAGCAAGATCaaaagcgaga                                                                       >  2:6525385/1‑149 (MQ=255)
                                                                             cACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCCCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagcgaga                                                                       <  1:441035/151‑3 (MQ=255)
                                                                               cgcgCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaag                                                                     >  1:2462018/1‑151 (MQ=255)
                                                                                gcgcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaagc                                                                    <  1:2823426/151‑1 (MQ=255)
                                                                                  gcTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaagcga                                                                  >  2:5923131/1‑151 (MQ=255)
                                                                                     gCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaagcgaatt                                                               >  2:6391878/1‑149 (MQ=255)
                                                                                     gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCCTCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGCAGTCGCAATCTGCCTCGCCTGAACTACCGCGCGCTTAGCAGCACGATCaaaagcgaaaagcgaaaa                                                               >  2:4599368/1‑151 (MQ=14)
                                                                                     gCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgaaaagcgaaaa                                                               >  1:3977054/1‑151 (MQ=255)
                                                                                      cAGGCCCCATCCGTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgaaaagcgaaaag                                                              <  2:5423607/151‑1 (MQ=255)
                                                                                      cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgaaaagcgaaaag                                                              >  1:5699640/1‑151 (MQ=255)
                                                                                      cAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCaaaagcgaaa                                                                       >  2:371531/1‑142 (MQ=255)
                                                                                       aGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaagcgaattcg                                                             <  1:5317854/151‑5 (MQ=255)
                                                                                         gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaagcgagaagcgaattcggt                                                           <  2:3628162/151‑7 (MQ=17)
                                                                                         gCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAgcgc                                                           >  1:5834126/1‑151 (MQ=255)
                                                                                          ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCt                                                          >  1:3987214/1‑151 (MQ=255)
                                                                                           cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTc                                                         <  2:4405166/151‑1 (MQ=255)
                                                                                            ccatCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCa                                                        >  1:1086992/4‑151 (MQ=255)
                                                                                            ccatCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCa                                                        >  2:4351490/4‑151 (MQ=255)
                                                                                              atCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAgc                                                      >  1:3752514/2‑151 (MQ=255)
                                                                                                ccATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGAGCAAAAGCGAAAAGCGAAAAGCGCTCAgcgc                                                    >  2:2823464/1‑151 (MQ=255)
                                                                                                 cATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGcgcc                                                   <  2:1086992/151‑1 (MQ=255)
                                                                                                   tGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGcgccgc                                                 >  1:604925/1‑151 (MQ=255)
                                                                                                    ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCa                                                >  1:4139988/1‑151 (MQ=255)
                                                                                                        ccTCGGCGCTTGACGGGCACCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAACGAAAAGCGAAAAGCGCTCAGCGCCGCACCTg                                            <  1:3206903/151‑1 (MQ=255)
                                                                                                        ccTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTg                                            >  1:833179/1‑151 (MQ=255)
                                                                                                           cAGCGCTGGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCtt                                         <  2:4029521/151‑1 (MQ=255)
                                                                                                           agggaggtGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCtt                                         <  2:4559773/144‑1 (MQ=255)
                                                                                                              cgcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAg                                      >  1:1241533/1‑151 (MQ=255)
                                                                                                                   gACGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGttt                                 <  2:2926961/151‑1 (MQ=255)
                                                                                                                    aCGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTc                                <  2:3057336/151‑1 (MQ=255)
                                                                                                                     cGGGCATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCg                               >  2:2067183/1‑151 (MQ=255)
                                                                                                                         cATCCATGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCaa                           <  2:5097846/151‑1 (MQ=255)
                                                                                                                           tCCATGCCCGTCATCTCCCGCCGCGAGGCCGGCGCTCGGCCTCCTGAAGTCGCAATCTGCGGCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATc                         <  2:6073676/151‑1 (MQ=25)
                                                                                                                               tGCCCGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGa                     <  2:6348948/151‑1 (MQ=255)
                                                                                                                                  ccGTCATCTCCCTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGAATg                  >  2:1504149/1‑151 (MQ=255)
                                                                                                                                       agctcACTCCGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAAAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGACGCACCTGCTTGAGGGTATCGGCAATAAGGAATGCCAAc             <  2:1208956/149‑1 (MQ=25)
                                                                                                                                         ctccctccGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGAATGCCAACTc           <  2:2995202/151‑1 (MQ=255)
                                                                                                                                             ctccGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGAATGCCAACTCCAGc       <  2:1697797/148‑1 (MQ=255)
                                                                                                                                               ccGCAAGGCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCACTCAGGAATGCCAACTCCAGCGa     >  1:1836216/2‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCGTCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTTCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAAAAGCGAAAAGCGCTCAGCGCCGCACCTGCTTGAGGGTTTCGGCAATCAGGAATGCCAACTCCAGCGACTG  >  CP009974/3756518‑3756814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: