Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,323,370 G→T 5.6% pseudogene (4512/4605 nt) yfaS ← pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,323,3700GT 5.6% 1312.7 / 32.0 336pseudogene (4512/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (171/146);  new base T (8/11);  total (179/157)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 345 total aligned reads displayed.

GCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATT  >  CP009273/2323221‑2323519
                                                                                                                                                     |                                                                                                                                                     
gCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACg                                                                                                                                                       >  2:5042339/1‑150 (MQ=255)
   aGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAAt                                                                                                                                                    >  2:5016900/1‑150 (MQ=255)
   aGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAAt                                                                                                                                                    >  2:1752230/1‑150 (MQ=255)
     tgtTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTg                                                                                                                                                  >  1:2855142/1‑150 (MQ=255)
      gtTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGc                                                                                                                                                 >  2:5124959/1‑150 (MQ=255)
      gtTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGc                                                                                                                                                 <  1:698782/150‑1 (MQ=255)
       ttGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGcc                                                                                                                                                >  2:4112333/1‑150 (MQ=255)
        tGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGccc                                                                                                                                               >  2:3156077/1‑150 (MQ=255)
            agagTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCtttt                                                                                                                                           <  2:4921318/150‑1 (MQ=255)
            agagTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCtttt                                                                                                                                           <  2:4689058/150‑1 (MQ=255)
             gagTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCttttt                                                                                                                                          <  1:1834175/150‑1 (MQ=255)
              agTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTg                                                                                                                                         >  1:4388402/1‑150 (MQ=255)
               gTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGc                                                                                                                                        <  1:569974/150‑1 (MQ=255)
                    cAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGaa                                                                                                                                   >  2:2861897/1‑150 (MQ=255)
                    cAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGaa                                                                                                                                   >  2:3916296/1‑150 (MQ=255)
                        gCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGACgcgc                                                                                                                               >  1:5053794/1‑150 (MQ=255)
                        gCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAgcgc                                                                                                                               <  1:258909/150‑1 (MQ=255)
                         ccAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCgca                                                                                                                              <  2:3690852/150‑1 (MQ=255)
                          cAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCgcag                                                                                                                             <  2:2585128/150‑1 (MQ=255)
                          cAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCgcag                                                                                                                             >  2:605105/1‑150 (MQ=255)
                             tAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTTCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCaa                                                                                                                          <  2:378277/150‑1 (MQ=255)
                                ccAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATg                                                                                                                       >  2:766468/1‑150 (MQ=255)
                                 cAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATgg                                                                                                                      >  2:4270944/1‑150 (MQ=255)
                                  aGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGc                                                                                                                     <  2:4586290/150‑1 (MQ=255)
                                     cGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGa                                                                                                                  <  2:907557/150‑1 (MQ=255)
                                      gATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGaa                                                                                                                 <  2:2081130/150‑1 (MQ=255)
                                      gATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGaa                                                                                                                 >  1:4017033/1‑150 (MQ=255)
                                           tGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGa                                                                                                            <  2:1356663/150‑1 (MQ=255)
                                              gcCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCg                                                                                                         >  1:4167303/1‑150 (MQ=255)
                                                  gTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTcc                                                                                                     >  1:969802/1‑150 (MQ=255)
                                                   ttCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTccc                                                                                                    <  1:3514730/150‑1 (MQ=255)
                                                   ttCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTccc                                                                                                    <  1:3825464/150‑1 (MQ=255)
                                                     cACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGt                                                                                                  >  2:3975162/1‑150 (MQ=255)
                                                       cTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCa                                                                                                <  2:3353806/150‑1 (MQ=255)
                                                        ttATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAg                                                                                               >  1:4181151/1‑150 (MQ=255)
                                                        ttATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAg                                                                                               >  1:2472880/1‑150 (MQ=255)
                                                         tATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGt                                                                                              <  1:3939537/150‑1 (MQ=255)
                                                          aTTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGtt                                                                                             <  1:2193720/150‑1 (MQ=255)
                                                            ttCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCt                                                                                           <  2:4542423/150‑1 (MQ=255)
                                                            ttCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCt                                                                                           >  1:1916570/1‑150 (MQ=255)
                                                               aCCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTc                                                                                        >  1:3128282/1‑150 (MQ=255)
                                                               aCCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTc                                                                                        >  2:1916981/1‑150 (MQ=255)
                                                               aCCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTc                                                                                        >  2:1916414/1‑150 (MQ=255)
                                                                ccTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCa                                                                                       >  2:1819287/1‑150 (MQ=255)
                                                                 cTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCAc                                                                                      >  1:1730193/1‑150 (MQ=255)
                                                                  tGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGATCCCGTCAGTTCTTTCAcc                                                                                     >  1:3332907/1‑150 (MQ=255)
                                                                    cATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCgg                                                                                   >  2:82533/1‑150 (MQ=255)
                                                                    cATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCgg                                                                                   <  1:4998940/150‑1 (MQ=255)
                                                                      tCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCa                                                                                 >  2:2463690/1‑150 (MQ=255)
                                                                      tCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCa                                                                                 >  2:4877006/1‑150 (MQ=255)
                                                                          ggTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCat                                                                             >  1:2618735/1‑150 (MQ=255)
                                                                          ggTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCat                                                                             >  1:2962882/1‑150 (MQ=255)
                                                                             ccATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCAtata                                                                          >  1:2333329/1‑150 (MQ=255)
                                                                              cATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAg                                                                         >  2:3443642/1‑150 (MQ=255)
                                                                                 tCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGcc                                                                      >  2:2291547/1‑150 (MQ=255)
                                                                                     tCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGtt                                                                  <  1:2926446/150‑1 (MQ=255)
                                                                                      cccTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTc                                                                 >  2:2062930/1‑150 (MQ=255)
                                                                                      cccTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTc                                                                 <  2:2179866/150‑1 (MQ=255)
                                                                                      cccTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTc                                                                 <  2:2314495/150‑1 (MQ=255)
                                                                                        cTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGc                                                               >  1:1562042/1‑150 (MQ=255)
                                                                                         tGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGcc                                                              <  1:4823262/150‑1 (MQ=255)
                                                                                          gCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGccc                                                             <  2:1770055/150‑1 (MQ=255)
                                                                                            cGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCcat                                                           >  2:1106419/1‑150 (MQ=255)
                                                                                            cGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCcat                                                           <  2:3296003/150‑1 (MQ=255)
                                                                                              cAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCcattt                                                         <  1:391671/150‑1 (MQ=255)
                                                                                                   cTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCcatttcattt                                                    >  2:4765711/1‑150 (MQ=255)
                                                                                                        ctgcGAAGGTGCATAGGAACGCACATAACGAGCATGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTTACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGTCCAGTTCGCCCATTTCATTTCgcgc                                               >  1:664305/1‑150 (MQ=255)
                                                                                                         tgcGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCt                                              >  1:4556294/1‑150 (MQ=255)
                                                                                                          gcgcAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCtt                                             <  1:770949/150‑1 (MQ=255)
                                                                                                          gcgcAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCtt                                             >  2:16063/1‑150 (MQ=255)
                                                                                                            gcAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCtttt                                           <  2:4359282/150‑1 (MQ=255)
                                                                                                               ggTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTcc                                        <  1:3684167/150‑1 (MQ=255)
                                                                                                               ggTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTcc                                        <  1:1111795/150‑1 (MQ=255)
                                                                                                                  gCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGc                                     >  2:1309824/1‑150 (MQ=255)
                                                                                                                  gCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATGTCGCGCTTTTTCCAGc                                     >  2:386908/1‑150 (MQ=255)
                                                                                                                   cATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCa                                    <  1:2839429/150‑1 (MQ=255)
                                                                                                                    aTAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCaa                                   >  1:3755829/1‑150 (MQ=255)
                                                                                                                     tAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAAt                                  >  1:3212944/1‑150 (MQ=255)
                                                                                                                      aGGAACGCACCTAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAAtt                                 <  2:274261/150‑1 (MQ=255)
                                                                                                                        gAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGc                               >  1:1745220/1‑150 (MQ=255)
                                                                                                                         aaCGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGcc                              <  2:4939893/150‑1 (MQ=255)
                                                                                                                            gCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCct                           >  1:1784978/1‑150 (MQ=255)
                                                                                                                               caTAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCctgct                        >  2:1219872/1‑150 (MQ=255)
                                                                                                                                 tAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGt                      >  2:4388035/1‑150 (MQ=255)
                                                                                                                                  aaCGAGCAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGtt                     <  1:3160058/150‑1 (MQ=255)
                                                                                                                                  aaCGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGtt                     >  2:2985373/1‑150 (MQ=255)
                                                                                                                                   aCGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGttt                    >  1:4507991/1‑150 (MQ=255)
                                                                                                                                    cGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTc                   <  2:3556286/150‑1 (MQ=255)
                                                                                                                                    cGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTc                   <  1:881651/150‑1 (MQ=255)
                                                                                                                                     gAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTcg                  >  2:672528/1‑150 (MQ=255)
                                                                                                                                        cAGGAGGCAGAACTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCgcgg               >  1:3600806/1‑150 (MQ=255)
                                                                                                                                         aggaggCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCgcggc              >  2:1155434/1‑150 (MQ=255)
                                                                                                                                         aggaggCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCgcggc              >  1:2128330/1‑150 (MQ=255)
                                                                                                                                            aggCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGtt           >  2:4195840/1‑150 (MQ=255)
                                                                                                                                                aGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGt       >  1:1007081/1‑150 (MQ=255)
                                                                                                                                                 gAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGtt      >  2:1707068/1‑150 (MQ=255)
                                                                                                                                                  aaCTAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGttt     <  1:5097089/150‑1 (MQ=255)
                                                                                                                                                   aCGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTa    <  1:1072820/150‑1 (MQ=255)
                                                                                                                                                     gAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTAtt  >  1:3432714/1‑150 (MQ=255)
                                                                                                                                                     gAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTAtt  >  1:1369791/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCTGCGCAGGTGCATAGGAACGCACATAACGAGCAGGAGGCAGAACGAATTGCCCTTTTTGCGAGAAGCGCAGCAAATGGCGGAAAGTGACCGTTCCCGTCAGTTCTTTCACCGGCACCATATAGGCCAGTTCGCCCATTTCATTTCGCGCTTTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATT  >  CP009273/2323221‑2323519

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: