Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,906,578 C→A 7.9% A86S (GCT→TCT)  htpX ← putative endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,906,5780CA 7.9% 1072.1 / 49.9 297A86S (GCT→TCT) htpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (119/154);  new base A (11/13);  total (130/167)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 304 total aligned reads displayed.

CGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCA  >  CP009273/1906429‑1906727
                                                                                                                                                     |                                                                                                                                                     
cGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGc                                                                                                                                                       <  1:1515677/150‑1 (MQ=255)
  ccTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCt                                                                                                                                                     <  1:2164212/150‑1 (MQ=255)
   cTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTg                                                                                                                                                    >  2:2404188/1‑150 (MQ=255)
    tCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTgg                                                                                                                                                   <  1:4809516/150‑1 (MQ=255)
     cGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTggg                                                                                                                                                  >  1:442321/1‑150 (MQ=255)
     cGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTggg                                                                                                                                                  >  2:2456203/1‑150 (MQ=255)
      ggCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGt                                                                                                                                                 >  2:2912371/1‑150 (MQ=255)
       gCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGtt                                                                                                                                                >  1:1162537/1‑150 (MQ=255)
        cTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTg                                                                                                                                               >  2:5223392/1‑150 (MQ=255)
         ttCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGc                                                                                                                                              <  1:2261449/150‑1 (MQ=255)
           cATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTa                                                                                                                                            <  1:1852770/150‑1 (MQ=255)
              ccGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAg                                                                                                                                         <  1:1805364/150‑1 (MQ=255)
               cGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGACTGGGTTGCTACAGt                                                                                                                                        >  1:2984379/1‑150 (MQ=255)
                  gCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAtt                                                                                                                                     <  2:361317/150‑1 (MQ=255)
                   cTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTg                                                                                                                                    >  1:842169/1‑150 (MQ=255)
                   cTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTg                                                                                                                                    <  2:265581/150‑1 (MQ=255)
                    tCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                                                                                                   >  1:1714209/1‑150 (MQ=255)
                       tGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCa                                                                                                                                <  1:2591031/150‑1 (MQ=255)
                        gTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAg                                                                                                                               <  2:4273440/150‑1 (MQ=255)
                            tGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCaa                                                                                                                           <  1:4787028/150‑1 (MQ=255)
                            tGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGACTGGGTTGCTACAGTATTGACCAGCCaa                                                                                                                           <  2:354876/150‑1 (MQ=255)
                             gcagcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAAc                                                                                                                          >  1:1290320/1‑150 (MQ=255)
                                gcaAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGtt                                                                                                                       >  1:3625003/1‑150 (MQ=255)
                                  aaaCCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTcc                                                                                                                     >  2:2485208/1‑150 (MQ=255)
                                    aCCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCt                                                                                                                   <  1:158/150‑1 (MQ=255)
                                     ccGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                                                                                                  <  2:3473993/150‑1 (MQ=255)
                                       ggTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTc                                                                                                                <  1:3210740/150‑1 (MQ=255)
                                         tGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGt                                                                                                              >  2:803794/1‑150 (MQ=255)
                                           cTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTa                                                                                                            >  2:3103809/1‑150 (MQ=255)
                                             gACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTAcg                                                                                                          <  2:490464/150‑1 (MQ=255)
                                               cAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTAcgcg                                                                                                        <  1:5120980/150‑1 (MQ=255)
                                                 gCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGt                                                                                                      <  1:1752567/150‑1 (MQ=255)
                                                     ccAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCt                                                                                                  <  2:1104724/150‑1 (MQ=255)
                                                      cAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                                                                                                 <  2:734836/150‑1 (MQ=255)
                                                      cAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                                                                                                 <  2:3895578/150‑1 (MQ=255)
                                                        gagagGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGa                                                                                               <  2:2067937/150‑1 (MQ=255)
                                                           agGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCa                                                                                            <  2:4120561/150‑1 (MQ=255)
                                                           agGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCa                                                                                            <  1:3331374/150‑1 (MQ=255)
                                                           agGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTAGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTCTTGACCAGCCAACGTTCCCTATCGTTACGCGGTTGCTCGATCa                                                                                            >  2:4055434/1‑150 (MQ=255)
                                                            ggCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCAc                                                                                           <  1:3949239/150‑1 (MQ=255)
                                                               aTCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTc                                                                                        >  2:4740883/1‑150 (MQ=255)
                                                                 cACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTccc                                                                                      <  2:4382731/150‑1 (MQ=255)
                                                                    ggCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGc                                                                                   >  1:4763735/1‑150 (MQ=255)
                                                                      cgcgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCa                                                                                 >  2:156465/1‑150 (MQ=255)
                                                                      cgcgcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCa                                                                                 >  1:4237580/1‑150 (MQ=255)
                                                                         gcACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGACTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACa                                                                              <  1:261792/150‑1 (MQ=255)
                                                                           aCCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGa                                                                            >  2:2515114/1‑150 (MQ=255)
                                                                           aCCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGa                                                                            >  2:5022738/1‑150 (MQ=255)
                                                                            ccGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGAt                                                                           <  2:4538354/150‑1 (MQ=255)
                                                                             cGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATc                                                                          <  1:3175131/150‑1 (MQ=255)
                                                                             cGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATc                                                                          <  2:2349043/150‑1 (MQ=255)
                                                                              ggTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCg                                                                         >  2:128451/1‑150 (MQ=255)
                                                                                 tGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTaa                                                                      >  1:3683086/1‑150 (MQ=255)
                                                                                 tGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTaa                                                                      >  1:1392189/1‑150 (MQ=255)
                                                                                   cAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATg                                                                    <  1:3877799/150‑1 (MQ=255)
                                                                                      aaGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGcca                                                                 <  2:3997002/150‑1 (MQ=255)
                                                                                      aaGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGcca                                                                 >  1:2427003/1‑150 (MQ=255)
                                                                                         cGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGccatcc                                                              <  1:2146287/150‑1 (MQ=255)
                                                                                          gTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGccatcca                                                             >  1:764840/1‑150 (MQ=255)
                                                                                          gTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGccatcca                                                             >  1:3790232/1‑150 (MQ=255)
                                                                                            tGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCAtt                                                           >  1:4529950/1‑150 (MQ=255)
                                                                                            tGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCAtt                                                           >  1:2365396/1‑150 (MQ=255)
                                                                                            tGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCAtt                                                           <  1:943553/150‑1 (MQ=255)
                                                                                             gATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCAttt                                                          <  1:87889/150‑1 (MQ=255)
                                                                                                 tCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGAc                                                      <  2:1800118/150‑1 (MQ=255)
                                                                                                 tCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGAc                                                      <  2:3451624/150‑1 (MQ=255)
                                                                                                  ccGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACa                                                     >  2:1379369/1‑150 (MQ=255)
                                                                                                       gcATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGa                                                <  1:3549290/150‑1 (MQ=255)
                                                                                                        cATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGaa                                               >  2:3311038/1‑150 (MQ=255)
                                                                                                         aTGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAg                                              >  1:622061/1‑150 (MQ=255)
                                                                                                             tAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGaa                                          >  1:3040381/1‑150 (MQ=255)
                                                                                                              aGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCCTATGGACATCAGAAGCGaaa                                         >  2:4890638/1‑150 (MQ=255)
                                                                                                               gATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAAc                                        <  1:67878/150‑1 (MQ=255)
                                                                                                                 tAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGa                                      >  2:3329307/1‑150 (MQ=255)
                                                                                                                  aGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGaa                                     >  2:4430686/1‑150 (MQ=255)
                                                                                                                  aGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGaa                                     <  2:912495/150‑1 (MQ=255)
                                                                                                                    ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGAAATGCCATCCATTTGGACAACAGAAGCGAAACGAAgg                                   >  1:1003723/1‑150 (MQ=255)
                                                                                                                      aCTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACCTCAGAAGCGAAACGAAGGaa                                 >  2:3191373/1‑150 (MQ=255)
                                                                                                                         tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTATCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAacca                              >  2:178487/1‑150 (MQ=255)
                                                                                                                          gCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAaccac                             <  1:4228433/150‑1 (MQ=255)
                                                                                                                          gCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAaccac                             <  1:2801961/150‑1 (MQ=255)
                                                                                                                            ggCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAAccacca                           <  1:4056571/150‑1 (MQ=255)
                                                                                                                              cATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCaaa                         >  2:1845272/1‑150 (MQ=255)
                                                                                                                              cATAGCGATCCCCGCCTGACGAGACTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCaaa                         <  2:2294279/150‑1 (MQ=255)
                                                                                                                               aTAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGACGGAACCACCaaaa                        >  1:2732644/1‑150 (MQ=255)
                                                                                                                                 aGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAAcc                      >  1:2654464/1‑150 (MQ=255)
                                                                                                                                    gATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGaa                   >  2:670036/1‑150 (MQ=255)
                                                                                                                                    gATCCCCGCCTGACGAGCCTGGGGTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTGCGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAAAGCCATCCATTTGGACATAATAAGCGAAACAAAGGAACCACAAAAACCGaa                   >  2:3137661/1‑150 (MQ=255)
                                                                                                                                     aTCCCCGCCTGACTAGCCTGGGTTGCTACAGTATTGACAAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCAGCCAACAGATCGTAATGACATCCATTAGTACATCAGAATCGAAACGAAGGAACCACCAAAAACGAAc                  >  1:1967100/1‑150 (MQ=255)
                                                                                                                                        cccGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGc               >  1:5037566/1‑150 (MQ=255)
                                                                                                                                         ccGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCa              <  1:4482681/150‑1 (MQ=255)
                                                                                                                                           gCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAg            >  1:4808367/1‑150 (MQ=255)
                                                                                                                                           gCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAg            <  2:5185117/150‑1 (MQ=255)
                                                                                                                                             cTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGc          <  2:379772/150‑1 (MQ=255)
                                                                                                                                                aCGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCAt       <  1:987949/150‑1 (MQ=255)
                                                                                                                                                   aGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGat    <  1:5206915/150‑1 (MQ=255)
                                                                                                                                                    gCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGatc   <  1:4996747/150‑1 (MQ=255)
                                                                                                                                                    gCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGatc   <  2:1714209/150‑1 (MQ=255)
                                                                                                                                                    gCCTGGGTTGCTACAGTAGTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACGTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGatc   <  1:1192851/150‑1 (MQ=255)
                                                                                                                                                     ccTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGatca  >  1:4944747/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCA  >  CP009273/1906429‑1906727

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: