Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,978,003 C→A 37.7% E455D (GAG→GAT aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,978,0030CA37.7% 403.1 / 299.8 281E455D (GAG→GATaslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base C (94/81);  new base A (50/56);  total (144/137)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 283 total aligned reads displayed.

ACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAG  >  CP009273/3977855‑3978152
                                                                                                                                                    |                                                                                                                                                     
acacCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCt                                                                                                                                                      <  1:1006056/150‑1 (MQ=255)
 cacCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTc                                                                                                                                                     >  2:1975729/1‑150 (MQ=255)
  acCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCg                                                                                                                                                    >  1:3469248/1‑150 (MQ=255)
   ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                                                                                                                                   <  1:320129/150‑1 (MQ=255)
   ccGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCgg                                                                                                                                                   <  1:2633061/150‑1 (MQ=255)
     gacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGcc                                                                                                                                                 >  2:2286162/1‑150 (MQ=255)
     gacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGcc                                                                                                                                                 >  2:2430172/1‑150 (MQ=255)
     gacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGcc                                                                                                                                                 >  1:3754641/1‑150 (MQ=255)
      acgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                                                                                                <  2:2932207/150‑1 (MQ=255)
         aTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                                                                                                                             <  1:291944/150‑1 (MQ=255)
         aTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTAc                                                                                                                                             >  1:4637359/1‑150 (MQ=255)
           cccGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACgg                                                                                                                                           <  1:974462/150‑1 (MQ=255)
           cccGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACgg                                                                                                                                           <  1:901777/150‑1 (MQ=255)
              gccgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTa                                                                                                                                        >  1:1349318/1‑150 (MQ=255)
                cgTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGa                                                                                                                                      >  1:3273228/1‑150 (MQ=255)
                  tttGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACt                                                                                                                                    <  2:4076416/150‑1 (MQ=255)
                       aTTACTGTTCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCa                                                                                                                               >  2:174312/1‑150 (MQ=255)
                       aTTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCa                                                                                                                               >  1:4137107/1‑150 (MQ=255)
                            tgtgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGt                                                                                                                          >  2:2625868/1‑150 (MQ=255)
                              tgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTc                                                                                                                        >  2:2477362/1‑150 (MQ=255)
                                 cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                                                                                                                     >  1:2682013/1‑150 (MQ=255)
                                  ggTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAg                                                                                                                    <  2:1504880/150‑1 (MQ=255)
                                    tGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGga                                                                                                                  >  2:750747/1‑150 (MQ=255)
                                            cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt                                                                                                          <  2:3441033/150‑1 (MQ=255)
                                            cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGAAGTGATACTTGAACTCATCAATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt                                                                                                          <  1:1336936/150‑1 (MQ=255)
                                             ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATg                                                                                                         <  2:2322823/150‑1 (MQ=255)
                                               tGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc                                                                                                       <  2:450097/150‑1 (MQ=255)
                                                gATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCt                                                                                                      <  2:4524418/150‑1 (MQ=255)
                                                 atatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTg                                                                                                     <  2:4139496/150‑1 (MQ=255)
                                                  tatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCTGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTgg                                                                                                    >  1:130266/1‑150 (MQ=255)
                                                        ctctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCcaca                                                                                              >  1:1061978/1‑150 (MQ=255)
                                                        ctctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCcaca                                                                                              >  1:4636146/1‑150 (MQ=255)
                                                        ctctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCcaca                                                                                              <  2:1136856/150‑1 (MQ=255)
                                                              gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                                        >  2:2186186/1‑150 (MQ=255)
                                                              gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                                        <  1:3259994/150‑1 (MQ=255)
                                                               tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGa                                                                                       >  1:3743594/1‑150 (MQ=255)
                                                                 taAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATa                                                                                     >  1:3127494/1‑150 (MQ=255)
                                                                     cGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAgg                                                                                 <  2:2960327/150‑1 (MQ=255)
                                                                     cGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAgg                                                                                 >  2:1570453/1‑150 (MQ=255)
                                                                      gTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggt                                                                                <  2:522924/150‑1 (MQ=255)
                                                                       tAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtg                                                                               <  1:2974723/150‑1 (MQ=255)
                                                                         aGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtggt                                                                             >  2:3727783/1‑150 (MQ=255)
                                                                          ggTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtt                                                                            <  1:146463/150‑1 (MQ=255)
                                                                          ggTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtt                                                                            <  1:2964201/150‑1 (MQ=255)
                                                                           gTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGttt                                                                           <  1:3149398/150‑1 (MQ=255)
                                                                           gTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGttt                                                                           <  2:1643795/150‑1 (MQ=255)
                                                                             tGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTc                                                                         <  1:3236695/150‑1 (MQ=255)
                                                                                   aTCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACt                                                                   <  2:4527206/150‑1 (MQ=255)
                                                                                     cAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTaa                                                                 >  2:880294/1‑150 (MQ=255)
                                                                                      aGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTaaa                                                                >  2:3426881/1‑150 (MQ=255)
                                                                                      aGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTaaa                                                                >  2:1830356/1‑150 (MQ=255)
                                                                                       ggACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAAt                                                               <  1:1632777/150‑1 (MQ=255)
                                                                                        gACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAAtt                                                              >  2:2058729/1‑150 (MQ=255)
                                                                                         aCGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTc                                                             <  1:2510944/150‑1 (MQ=255)
                                                                                            tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGcc                                                          >  1:3646124/1‑150 (MQ=255)
                                                                                            tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGcc                                                          <  1:4398751/150‑1 (MQ=255)
                                                                                            tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGcc                                                          >  1:2126322/1‑150 (MQ=255)
                                                                                               tACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACt                                                       <  1:4212878/150‑1 (MQ=255)
                                                                                               tACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACt                                                       <  2:3012889/150‑1 (MQ=255)
                                                                                               tACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACt                                                       <  2:328581/150‑1 (MQ=255)
                                                                                                aCTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACtt                                                      <  1:710556/150‑1 (MQ=255)
                                                                                                aCTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACtt                                                      <  1:4486759/150‑1 (MQ=255)
                                                                                                   tGAACTCATCCATCCGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCg                                                   <  2:1224422/150‑1 (MQ=255)
                                                                                                      aCTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTc                                                >  2:3678490/1‑150 (MQ=255)
                                                                                                       cTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTcc                                               >  1:589815/1‑150 (MQ=255)
                                                                                                        tCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCa                                              >  2:331535/1‑150 (MQ=255)
                                                                                                        tCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCa                                              <  2:3351638/150‑1 (MQ=255)
                                                                                                         cATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAATGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAg                                             >  2:1583906/1‑150 (MQ=255)
                                                                                                         cATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAg                                             <  1:1818439/150‑1 (MQ=255)
                                                                                                            ccATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGAt                                          <  1:4617041/150‑1 (MQ=255)
                                                                                                               tACGCACCGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGCTAAAGGGGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGcc                                       >  2:1041903/1‑150 (MQ=255)
                                                                                                               tACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGcc                                       <  2:4356075/150‑1 (MQ=255)
                                                                                                                aCGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGccc                                      >  2:3660829/1‑150 (MQ=255)
                                                                                                                    acaGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGcc                                  >  1:509971/1‑150 (MQ=255)
                                                                                                                     cagcagCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCa                                 >  1:1284638/1‑150 (MQ=255)
                                                                                                                        cagcGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGAt                              >  1:2076619/1‑150 (MQ=255)
                                                                                                                           cGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCa                           <  2:653665/150‑1 (MQ=255)
                                                                                                                           cGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCa                           >  2:3549474/1‑150 (MQ=255)
                                                                                                                             aGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCTCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAgc                         >  1:1531659/1‑150 (MQ=255)
                                                                                                                              gTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAgcg                        <  2:1349318/150‑1 (MQ=255)
                                                                                                                                ttACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGcgcg                      >  1:1372652/1‑150 (MQ=255)
                                                                                                                                 tACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCgg                     <  1:246002/150‑1 (MQ=255)
                                                                                                                                  aCCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATAGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGt                    >  2:4011284/1‑150 (MQ=255)
                                                                                                                                     gTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAgg                 <  1:835292/150‑1 (MQ=255)
                                                                                                                                       tGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATACAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGa               <  2:1181716/150‑1 (MQ=255)
                                                                                                                                       tGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTCAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGa               >  2:2092085/1‑150 (MQ=255)
                                                                                                                                        gAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGaa              <  2:2621007/150‑1 (MQ=255)
                                                                                                                                         aGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAag             <  1:1145350/150‑1 (MQ=255)
                                                                                                                                          ggAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAaga            <  2:1453325/150‑1 (MQ=255)
                                                                                                                                           gAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAagag           <  1:1904992/150‑1 (MQ=255)
                                                                                                                                           gAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAagag           >  1:1688026/1‑150 (MQ=255)
                                                                                                                                            aaGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAgaga          <  1:4001098/150‑1 (MQ=255)
                                                                                                                                            aaGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAgaga          <  2:3926223/150‑1 (MQ=255)
                                                                                                                                            aaGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAgaga          <  2:1061978/150‑1 (MQ=255)
                                                                                                                                            aaGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAgaga          >  2:2731484/1‑150 (MQ=255)
                                                                                                                                             agtagtGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGAt         <  2:4306112/150‑1 (MQ=255)
                                                                                                                                                agtGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTg      <  1:1009051/150‑1 (MQ=255)
                                                                                                                                                  tGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGcc    >  1:926386/1‑150 (MQ=255)
                                                                                                                                                  tGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGcc    >  1:1480262/1‑150 (MQ=255)
                                                                                                                                                    cTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAg  <  2:2885509/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAG  >  CP009273/3977855‑3978152

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: