Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,856,944 C→A 13.0% pseudogene (20/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,856,9440CA13.0% 976.5 / 89.9 297pseudogene (20/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (115/143);  new base A (21/17);  total (136/161)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 299 total aligned reads displayed.

CAATAATGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCTTTCATG  >  CP009273/3856795‑3857092
                                                                                                                                                     |                                                                                                                                                    
cAATAATGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGc                                                                                                                                                      <  1:3220777/150‑1 (MQ=255)
  ataatGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGt                                                                                                                                                    >  1:757146/1‑150 (MQ=255)
  ataatGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGt                                                                                                                                                    <  1:3323984/150‑1 (MQ=255)
  ataatGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGt                                                                                                                                                    >  1:4428400/1‑150 (MQ=255)
     atGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGa                                                                                                                                                 <  2:2424437/150‑1 (MQ=255)
       gtgtACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAAt                                                                                                                                               <  1:1365076/150‑1 (MQ=255)
       gtgtACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAAt                                                                                                                                               >  2:1344310/1‑150 (MQ=255)
       gtgtACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAAt                                                                                                                                               <  2:1461582/150‑1 (MQ=255)
       gtgtACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAAt                                                                                                                                               >  2:1169840/1‑150 (MQ=255)
        tgtACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAAtt                                                                                                                                              <  2:1481870/150‑1 (MQ=255)
          tACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTg                                                                                                                                            >  1:4398244/1‑150 (MQ=255)
           aCGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGa                                                                                                                                           <  2:3357036/150‑1 (MQ=255)
              aTTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTg                                                                                                                                        >  1:360960/1‑150 (MQ=255)
                ttGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAg                                                                                                                                      <  2:2937345/150‑1 (MQ=255)
                  gcgcGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCa                                                                                                                                    <  2:3114714/150‑1 (MQ=255)
                   cgcgAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCAt                                                                                                                                   <  1:11449/150‑1 (MQ=255)
                          aaaCTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCcaca                                                                                                                            >  2:3068954/1‑150 (MQ=255)
                              tGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCct                                                                                                                        >  2:3898964/1‑150 (MQ=255)
                               gTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCctc                                                                                                                       <  1:3437781/150‑1 (MQ=255)
                                  cGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCttt                                                                                                                    <  1:4522574/150‑1 (MQ=255)
                                  cGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCttt                                                                                                                    >  2:4424239/1‑150 (MQ=255)
                                   ggATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCtttt                                                                                                                   <  2:2472233/150‑1 (MQ=255)
                                     aTCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTtct                                                                                                                 <  1:2848751/150‑1 (MQ=255)
                                     aTCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTtct                                                                                                                 >  2:4519914/1‑150 (MQ=255)
                                        ggTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCaa                                                                                                              >  2:2496729/1‑150 (MQ=255)
                                          tCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAAtt                                                                                                            >  1:2869029/1‑150 (MQ=255)
                                           cAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTc                                                                                                           >  1:3807102/1‑150 (MQ=255)
                                           cAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTc                                                                                                           >  1:4465106/1‑150 (MQ=255)
                                               gATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTa                                                                                                       >  1:3953688/1‑150 (MQ=255)
                                                 ttCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAg                                                                                                     >  2:1052809/1‑150 (MQ=255)
                                                    cccGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAg                                                                                                  <  2:4190181/150‑1 (MQ=255)
                                                      cGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTCAGTCGGa                                                                                                >  1:883879/1‑150 (MQ=255)
                                                      cGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGa                                                                                                <  2:757146/150‑1 (MQ=255)
                                                          aaaCATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGcc                                                                                            <  2:2745515/150‑1 (MQ=255)
                                                          aaaCATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGAACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGcc                                                                                            <  1:2686568/150‑1 (MQ=255)
                                                           aaCATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCt                                                                                           <  2:2813857/150‑1 (MQ=255)
                                                            aCATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTg                                                                                          >  1:2249500/1‑150 (MQ=255)
                                                             cATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGa                                                                                         >  1:1038459/1‑150 (MQ=255)
                                                                cGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTg                                                                                      <  2:1236070/150‑1 (MQ=255)
                                                                cGGGTTTTGCAGCAAGATGGAAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTg                                                                                      <  2:2372156/150‑1 (MQ=255)
                                                                  ggTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGt                                                                                    >  1:2570893/1‑150 (MQ=255)
                                                                   gTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGtt                                                                                   >  2:2990080/1‑150 (MQ=255)
                                                                     tttGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTAt                                                                                 <  2:4021510/150‑1 (MQ=255)
                                                                           gcaAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCa                                                                           <  2:4324351/150‑1 (MQ=255)
                                                                             aaGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAgg                                                                         >  2:1291460/1‑150 (MQ=255)
                                                                              aGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGc                                                                        >  1:1626106/1‑150 (MQ=255)
                                                                                  ggCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTgg                                                                    <  2:2198556/150‑1 (MQ=255)
                                                                                    cAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAg                                                                  >  2:51721/1‑150 (MQ=255)
                                                                                     aaGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGc                                                                 >  1:1986325/1‑150 (MQ=255)
                                                                                     aaGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGc                                                                 <  1:1525124/150‑1 (MQ=255)
                                                                                      aGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCg                                                                >  1:2805849/1‑150 (MQ=255)
                                                                                        aCCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAg                                                              <  2:3099500/150‑1 (MQ=255)
                                                                                        aCCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAg                                                              <  1:600586/150‑1 (MQ=255)
                                                                                          ccaccaCAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGa                                                            >  2:1214820/1‑150 (MQ=255)
                                                                                           caccacAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGAt                                                           <  2:2659783/150‑1 (MQ=255)
                                                                                            accacAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATa                                                          <  2:3959969/150‑1 (MQ=255)
                                                                                            accacAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATa                                                          >  2:54817/1‑150 (MQ=255)
                                                                                              cacaATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAgg                                                        <  1:265568/150‑1 (MQ=255)
                                                                                                caATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGaa                                                      <  2:616997/150‑1 (MQ=255)
                                                                                                 aaTCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAg                                                     <  1:3055038/150‑1 (MQ=255)
                                                                                                  aTCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGa                                                    <  2:3233437/150‑1 (MQ=255)
                                                                                                    cccGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATg                                                  >  1:4203575/1‑150 (MQ=255)
                                                                                                    cccGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATg                                                  <  1:337451/150‑1 (MQ=255)
                                                                                                     ccGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGc                                                 <  2:1103644/150‑1 (MQ=255)
                                                                                                       ggCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGa                                               <  2:1552127/150‑1 (MQ=255)
                                                                                                         cGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATc                                             <  2:3411388/150‑1 (MQ=255)
                                                                                                         cGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATc                                             <  1:2346527/150‑1 (MQ=255)
                                                                                                          gAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCg                                            <  2:2680000/150‑1 (MQ=255)
                                                                                                           aaGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCgg                                           >  2:1849450/1‑150 (MQ=255)
                                                                                                                 aaCAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAAc                                     <  2:1210518/150‑1 (MQ=255)
                                                                                                                 aaCAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAAc                                     >  1:3102201/1‑150 (MQ=255)
                                                                                                                   cagcagCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGATTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACAt                                   <  2:1521611/150‑1 (MQ=255)
                                                                                                                    agcagcACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATg                                  <  1:4266316/150‑1 (MQ=255)
                                                                                                                      cagcaCTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATgac                                <  2:3588328/150‑1 (MQ=255)
                                                                                                                      cagcaCTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGAGAGGAAGATGCGATCGGCAGAACATgac                                >  1:2208887/1‑150 (MQ=255)
                                                                                                                       agcaCTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATgacg                               <  2:1903090/150‑1 (MQ=255)
                                                                                                                          aCTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGAcc                            <  1:2304559/150‑1 (MQ=255)
                                                                                                                          aCTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGAcc                            <  1:39837/150‑1 (MQ=255)
                                                                                                                            tGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCa                          <  2:1744696/150‑1 (MQ=255)
                                                                                                                             ggCGTGAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCaa                         >  2:3695231/1‑150 (MQ=255)
                                                                                                                              gCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAAt                        <  1:1002991/150‑1 (MQ=255)
                                                                                                                              gCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAAt                        <  1:3731923/150‑1 (MQ=255)
                                                                                                                                 tGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATggg                     <  2:3277943/150‑1 (MQ=255)
                                                                                                                                  gAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGa                    >  1:3549925/1‑150 (MQ=255)
                                                                                                                                  gAACATCGCGCCGCCAAAGAGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGa                    <  1:1300105/150‑1 (MQ=255)
                                                                                                                                     cATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTc                 <  2:3868023/150‑1 (MQ=255)
                                                                                                                                     cATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTc                 <  1:863643/150‑1 (MQ=255)
                                                                                                                                     cATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTc                 <  1:51728/150‑1 (MQ=255)
                                                                                                                                     cATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTc                 <  2:1588494/150‑1 (MQ=255)
                                                                                                                                      aTCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCg                <  1:2179825/150‑1 (MQ=255)
                                                                                                                                      aTCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCg                >  1:2427856/1‑150 (MQ=255)
                                                                                                                                      aTCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGACTGATTGGTTATGCGGAAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCg                >  1:2938388/1‑150 (MQ=255)
                                                                                                                                       tCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGt               <  1:316527/150‑1 (MQ=255)
                                                                                                                                        cgcgCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTg              <  1:1839304/150‑1 (MQ=255)
                                                                                                                                        cgcgCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTg              >  1:2189161/1‑150 (MQ=255)
                                                                                                                                          cgccgcCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGAt            <  2:476519/150‑1 (MQ=255)
                                                                                                                                              gccAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCt        >  2:3542642/1‑150 (MQ=255)
                                                                                                                                              gccAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGAGTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCt        >  2:2609013/1‑150 (MQ=255)
                                                                                                                                                   aGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCTTTCAt   >  1:2125959/1‑150 (MQ=255)
                                                                                                                                                   aGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCTTTCAt   >  2:888928/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATAATGTGTACGATTTGCGCGAATAAACTGTTCGGATCGGTCAGTGATTCCCCGACAAACATCGGGTTTTGCAGCAAGATGGCAAGACCCACCACAATCCCGGCGAAGGGAAACAGCAGCACTGGCGTGAACATCGCGCCGCCAAAGCGTTGAATTTGACTGAGCATTTCCACATCCTCTTTTCTCAATTCTGTAAGTAGGAAGCCTGATTGGTTATGCGGCAGGCCTGGAGCGAGGATAGGAAGATGCGATCGGCAGAACATGACGACCCAATGGGATTCGTGATCGCTTTCATG  >  CP009273/3856795‑3857092

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: