Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 274,735 A→T 100% pseudogene (138/348 nt) afuB ← pseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273274,7350AT95.7% 680.6 / 13.4 206pseudogene (138/348 nt)afuBpseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (6/3);  new base T (106/90);  total (113/93)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

GAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGA  >  CP009273/274588‑274881
                                                                                                                                                   |                                                                                                                                                  
  aaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATgg                                                                                                                                                >  2:4319136/1‑150 (MQ=255)
   aataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgc                                                                                                                                               <  2:4452812/150‑1 (MQ=255)
     taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgc                                                                                                                                             >  2:500612/1‑150 (MQ=255)
     taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgc                                                                                                                                             >  2:500666/1‑150 (MQ=255)
      aatCGCCAGCATCACCACGATCAGAATGGAACCGGAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgcg                                                                                                                                            <  2:593629/150‑1 (MQ=255)
        tCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGca                                                                                                                                          >  2:2845994/1‑150 (MQ=255)
        tCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGca                                                                                                                                          <  2:1790345/150‑1 (MQ=255)
         cGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcac                                                                                                                                         >  2:3008128/1‑150 (MQ=255)
           ccAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACaa                                                                                                                                       >  2:2874445/1‑150 (MQ=255)
            cAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACaaa                                                                                                                                      >  2:1458819/1‑150 (MQ=255)
             aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAg                                                                                                                                     <  2:325939/150‑1 (MQ=255)
             aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAg                                                                                                                                     >  1:4501034/1‑150 (MQ=255)
              gCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGc                                                                                                                                    >  1:3255707/1‑150 (MQ=255)
              gCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGc                                                                                                                                    >  1:2680297/1‑150 (MQ=255)
               cATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGAATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCt                                                                                                                                   >  1:4462024/1‑150 (MQ=255)
                aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                  >  1:3551992/1‑150 (MQ=255)
                aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                  >  2:3815695/1‑150 (MQ=255)
                aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                  >  2:2210139/1‑150 (MQ=255)
                  caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTa                                                                                                                                >  1:3211386/1‑150 (MQ=255)
                  caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTa                                                                                                                                >  1:3211390/1‑150 (MQ=255)
                    ccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTaaa                                                                                                                              >  2:97137/1‑150 (MQ=255)
                    ccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTGATGGCGCGCACAAAGCTGTaaa                                                                                                                              >  2:3901148/1‑150 (MQ=255)
                     cacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAAt                                                                                                                             >  2:3070647/1‑150 (MQ=255)
                       cGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCa                                                                                                                           >  1:2655384/1‑150 (MQ=255)
                       cGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCa                                                                                                                           >  1:1537682/1‑150 (MQ=255)
                         aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                         >  1:2881246/1‑150 (MQ=255)
                           cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                       >  1:2770096/1‑150 (MQ=255)
                             gAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCg                                                                                                                     <  1:4250859/150‑1 (MQ=255)
                             gAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCg                                                                                                                     <  1:4248583/150‑1 (MQ=255)
                              aaTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCga                                                                                                                    <  1:1071611/150‑1 (MQ=255)
                               aTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCgag                                                                                                                   >  1:1914099/1‑150 (MQ=255)
                               aTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCgag                                                                                                                   >  1:771604/1‑150 (MQ=255)
                                 ggAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGagag                                                                                                                 <  1:4238760/150‑1 (MQ=255)
                                 ggAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGagag                                                                                                                 <  1:3749545/150‑1 (MQ=255)
                                  gAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAgg                                                                                                                <  2:3079337/150‑1 (MQ=255)
                                  gAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAgg                                                                                                                <  2:1752843/150‑1 (MQ=255)
                                  gAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAgg                                                                                                                >  2:1500660/1‑150 (MQ=255)
                                     ccGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATc                                                                                                             <  1:1420715/150‑1 (MQ=255)
                                      cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                            >  2:4111351/1‑150 (MQ=255)
                                      cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                            >  2:1858898/1‑150 (MQ=255)
                                        tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGATCACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                          <  1:4416354/150‑1 (MQ=255)
                                          ggCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCgg                                                                                                        <  1:1626941/150‑1 (MQ=255)
                                           gCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCggg                                                                                                       <  1:280010/150‑1 (MQ=255)
                                            cAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgc                                                                                                      >  2:3386860/1‑150 (MQ=255)
                                             aaTCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgcg                                                                                                     >  1:998831/1‑150 (MQ=255)
                                                 gCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGACCACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagc                                                                                                 >  2:1329541/1‑150 (MQ=255)
                                                 gCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagc                                                                                                 >  2:403082/1‑150 (MQ=255)
                                                  cTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagca                                                                                                <  1:2998160/150‑1 (MQ=255)
                                                     cacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCTAGAGGATCGCCGGGCGCAGCAGCg                                                                                             >  1:2393301/1‑150 (MQ=255)
                                                     cacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCg                                                                                             <  1:1827057/150‑1 (MQ=255)
                                                       ccGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgc                                                                                           >  2:656229/1‑150 (MQ=255)
                                                         gTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcag                                                                                         >  1:1226138/1‑150 (MQ=255)
                                                           aTTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                       >  2:2460265/1‑150 (MQ=255)
                                                           aTTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                       >  1:2634604/1‑150 (MQ=255)
                                                            ttCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGcagcag                                                                                      >  2:3820367/1‑150 (MQ=255)
                                                               gCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAt                                                                                   <  1:2144467/150‑1 (MQ=255)
                                                                 cATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgt                                                                                 >  2:1475416/1‑150 (MQ=255)
                                                                        cACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAAt                                                                          <  2:3214129/150‑1 (MQ=255)
                                                                           gcgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGt                                                                       <  2:3504832/150‑1 (MQ=255)
                                                                             ggTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTAc                                                                     >  1:4035748/1‑150 (MQ=255)
                                                                                tGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCa                                                                  >  2:2852001/1‑150 (MQ=255)
                                                                                 gAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAg                                                                 <  1:2439142/150‑1 (MQ=255)
                                                                                 gAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAg                                                                 <  1:476727/150‑1 (MQ=255)
                                                                                   ggATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAggg                                                               >  2:1868126/1‑150 (MQ=255)
                                                                                       gTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAcc                                                           >  1:4354598/1‑150 (MQ=255)
                                                                                         aaGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCgg                                                         <  1:2262526/150‑1 (MQ=255)
                                                                                           gCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcg                                                       >  1:3100709/1‑150 (MQ=255)
                                                                                           gCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcg                                                       <  2:3689208/150‑1 (MQ=255)
                                                                                                cGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAg                                                  >  2:1089128/1‑150 (MQ=255)
                                                                                                 gCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAgg                                                 >  2:1432820/1‑150 (MQ=255)
                                                                                                   tACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCt                                               <  2:509796/150‑1 (MQ=255)
                                                                                                                  cGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                >  1:1769850/1‑150 (MQ=255)
                                                                                                                  cGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                <  2:793703/150‑1 (MQ=255)
                                                                                                                  cGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                >  2:2530751/1‑150 (MQ=255)
                                                                                                                   gTGACGAGGAACACACTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCg                               <  2:2849368/150‑1 (MQ=255)
                                                                                                                   gTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCg                               <  1:2193726/150‑1 (MQ=255)
                                                                                                                        gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                          <  2:1206124/150‑1 (MQ=255)
                                                                                                                          ggAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAtt                        <  1:1609715/150‑1 (MQ=255)
                                                                                                                           gAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAttt                       <  1:1128383/150‑1 (MQ=255)
                                                                                                                              cacaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                    <  1:2211645/150‑1 (MQ=255)
                                                                                                                               acaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCg                   >  2:1949203/1‑150 (MQ=255)
                                                                                                                               acaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCg                   <  1:3200298/150‑1 (MQ=255)
                                                                                                                                 aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                 <  1:500666/150‑1 (MQ=255)
                                                                                                                                 aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                 <  1:500612/150‑1 (MQ=255)
                                                                                                                                 aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                 >  1:3523743/1‑150 (MQ=255)
                                                                                                                                   tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt               <  2:2062520/150‑1 (MQ=255)
                                                                                                                                   tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt               <  1:328302/150‑1 (MQ=255)
                                                                                                                                    ggCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTg              <  1:739446/150‑1 (MQ=255)
                                                                                                                                       gcTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGcc           >  1:4321196/1‑150 (MQ=255)
                                                                                                                                        cTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGccc          >  1:3729141/1‑150 (MQ=255)
                                                                                                                                         tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         <  2:2402927/150‑1 (MQ=255)
                                                                                                                                         tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         <  1:1459478/150‑1 (MQ=255)
                                                                                                                                         tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         <  1:2845994/150‑1 (MQ=255)
                                                                                                                                         tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         >  2:2580206/1‑150 (MQ=255)
                                                                                                                                           aCGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGt       <  2:3936356/150‑1 (MQ=255)
                                                                                                                                             ggtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTaa     <  1:3757649/150‑1 (MQ=255)
                                                                                                                                             ggtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTaa     <  2:2227847/150‑1 (MQ=255)
                                                                                                                                              gtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt    >  2:72581/1‑150 (MQ=255)
                                                                                                                                               tggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATg   <  1:3422948/150‑1 (MQ=255)
                                                                                                                                                ggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGa  >  1:3766744/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGA  >  CP009273/274588‑274881

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: