Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F218 I0 R1 1506 278.7 9275112 94.5% 8764980 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 619,048 C→A 93.4% A307S (GCC→TCC)  fepB ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273619,0480CA93.4% 822.5 / 47.3 271A307S (GCC→TCC) fepBiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (123/130);  total (130/141)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 274 total aligned reads displayed.

CGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  CP009273/618900‑619195
                                                                                                                                                    |                                                                                                                                                   
 gAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGgagc                                                                                                                                                   <  2:4422484/150‑1 (MQ=255)
  aGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAATCCTATCCAGCACTTGCATGGAGCt                                                                                                                                                  >  1:396559/1‑150 (MQ=255)
   gCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTg                                                                                                                                                 >  1:1285216/1‑150 (MQ=255)
   gCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTg                                                                                                                                                 >  1:53500/1‑150 (MQ=255)
      gTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtag                                                                                                                                              >  2:1754311/1‑150 (MQ=255)
       ttttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtagt                                                                                                                                             >  2:4211320/1‑150 (MQ=255)
         ttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTaa                                                                                                                                           >  1:999794/1‑150 (MQ=255)
           ggTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATc                                                                                                                                         <  2:4132559/150‑1 (MQ=255)
                   gATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGa                                                                                                                                 <  2:817462/150‑1 (MQ=255)
                   gATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGa                                                                                                                                 >  2:4429348/1‑150 (MQ=255)
                    aTTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCATAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGaa                                                                                                                                >  2:935273/1‑150 (MQ=255)
                      tATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACg                                                                                                                              <  2:91893/150‑1 (MQ=255)
                        tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGtc                                                                                                                            <  2:2530427/150‑1 (MQ=255)
                        tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtc                                                                                                                            <  2:2077674/150‑1 (MQ=255)
                        tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtc                                                                                                                            <  2:1950295/150‑1 (MQ=255)
                          ggCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtctc                                                                                                                          <  1:339699/150‑1 (MQ=255)
                           gCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCg                                                                                                                         >  1:2708438/1‑150 (MQ=255)
                            cGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCgg                                                                                                                        <  2:1676305/150‑1 (MQ=255)
                              tgttgtTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGtt                                                                                                                      <  1:4146456/150‑1 (MQ=255)
                               gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTc                                                                                                                     >  2:3611720/1‑150 (MQ=255)
                               gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTc                                                                                                                     >  2:3611681/1‑150 (MQ=255)
                               gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTc                                                                                                                     <  2:3035964/150‑1 (MQ=255)
                                 tgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTccc                                                                                                                   <  2:3182648/150‑1 (MQ=255)
                                  gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                                  <  1:958340/150‑1 (MQ=255)
                                  gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                                  <  1:364831/150‑1 (MQ=255)
                                  gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                                  <  1:4636175/150‑1 (MQ=255)
                                          tggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                          >  2:4024956/1‑150 (MQ=255)
                                              gctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCt                                                                                                      <  1:3781480/150‑1 (MQ=255)
                                               ctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGTAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTg                                                                                                     >  2:1728163/1‑150 (MQ=255)
                                                tggtggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGc                                                                                                    >  2:727722/1‑150 (MQ=255)
                                                    ggAGTTGCGACATTTTCGACAGACGCCGCCGCAGGTGACAGCGTCAGACAGTTAATGCTTAAATCATCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGt                                                                                                <  1:4187253/150‑1 (MQ=255)
                                                       gTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGtttt                                                                                             >  1:1156764/1‑150 (MQ=255)
                                                       gTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGtttt                                                                                             <  2:139984/150‑1 (MQ=255)
                                                         tGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGt                                                                                           <  1:4535916/150‑1 (MQ=255)
                                                         tGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGt                                                                                           <  1:1096363/150‑1 (MQ=255)
                                                           cGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTAc                                                                                         <  2:2688127/150‑1 (MQ=255)
                                                             aCATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTg                                                                                       <  2:3133584/150‑1 (MQ=255)
                                                             aCATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTg                                                                                       >  1:3499794/1‑150 (MQ=255)
                                                             aCATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTg                                                                                       >  1:3499694/1‑150 (MQ=255)
                                                             aCATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTg                                                                                       >  2:1939540/1‑150 (MQ=255)
                                                                ttttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgcag                                                                                    <  1:4177551/150‑1 (MQ=255)
                                                                   tCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgcaggca                                                                                 >  2:2969392/1‑150 (MQ=255)
                                                                       cAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTg                                                                             <  2:188834/150‑1 (MQ=255)
                                                                        aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                            >  2:3768812/1‑150 (MQ=255)
                                                                           cgccgccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCGTGTTTTGTACTGCAGGCAGGTGcgcg                                                                         >  1:2125486/1‑150 (MQ=255)
                                                                            gccgccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGa                                                                        <  1:245525/150‑1 (MQ=255)
                                                                               gccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagca                                                                     >  1:3238967/1‑150 (MQ=255)
                                                                               gccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagca                                                                     <  1:1045520/150‑1 (MQ=255)
                                                                                 cgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagcagc                                                                   >  2:2016589/1‑150 (MQ=255)
                                                                                 cgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagcagc                                                                   <  1:70032/150‑1 (MQ=255)
                                                                                 cgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTCTTTTGTACTGCAGGCAGGTGCGCGagcagc                                                                   <  1:2380485/150‑1 (MQ=255)
                                                                                   cAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCgg                                                                 <  1:2043263/150‑1 (MQ=255)
                                                                                    aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                                >  2:3296187/1‑150 (MQ=255)
                                                                                    aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                                >  2:3862722/1‑150 (MQ=255)
                                                                                    aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                                <  2:1313188/150‑1 (MQ=255)
                                                                                    aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                                <  2:1313258/150‑1 (MQ=255)
                                                                                      gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAtt                                                              <  1:1995683/150‑1 (MQ=255)
                                                                                      gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAtt                                                              <  1:748371/150‑1 (MQ=255)
                                                                                       tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                             >  1:3155151/1‑150 (MQ=255)
                                                                                       tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                             <  2:3595650/150‑1 (MQ=255)
                                                                                          cAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCa                                                          >  1:4366982/1‑150 (MQ=255)
                                                                                           aGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCAt                                                         >  1:1717791/1‑150 (MQ=255)
                                                                                            gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATa                                                        <  1:2951762/150‑1 (MQ=255)
                                                                                             cGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                       <  1:4298085/150‑1 (MQ=255)
                                                                                                 cGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAg                                                   <  2:3327711/150‑1 (MQ=255)
                                                                                                 cGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAg                                                   >  1:2601274/1‑150 (MQ=255)
                                                                                                  gACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                                  <  2:120293/150‑1 (MQ=255)
                                                                                                    cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                                <  1:1408807/150‑1 (MQ=255)
                                                                                                    cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                                <  1:1448779/150‑1 (MQ=255)
                                                                                                     aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                               <  2:2499084/150‑1 (MQ=255)
                                                                                                     aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACGTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                               <  2:2278925/150‑1 (MQ=255)
                                                                                                        tAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGc                                            >  2:109647/1‑150 (MQ=255)
                                                                                                         aaTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                           >  2:355147/1‑150 (MQ=255)
                                                                                                                aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                    >  2:996228/1‑150 (MQ=255)
                                                                                                                aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                    <  1:3796040/150‑1 (MQ=255)
                                                                                                                 aaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTg                                   <  2:3436299/150‑1 (MQ=255)
                                                                                                                 aaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTg                                   >  1:2993110/1‑150 (MQ=255)
                                                                                                                  aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGa                                  <  2:3990414/150‑1 (MQ=255)
                                                                                                                   aCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGAt                                 <  2:3481026/150‑1 (MQ=255)
                                                                                                                    cAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATc                                <  1:3445475/150‑1 (MQ=255)
                                                                                                                    cAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATc                                <  1:4029154/150‑1 (MQ=255)
                                                                                                                    cAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATc                                >  1:2133647/1‑150 (MQ=255)
                                                                                                                        gcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCg                            <  2:2825226/150‑1 (MQ=255)
                                                                                                                          cTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGc                          >  2:2436534/1‑150 (MQ=255)
                                                                                                                              aGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                      <  1:2516074/150‑1 (MQ=255)
                                                                                                                              aGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                      <  2:826664/150‑1 (MQ=255)
                                                                                                                              aGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                      >  1:4576871/1‑150 (MQ=255)
                                                                                                                                  tATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGa                  <  2:2591389/150‑1 (MQ=255)
                                                                                                                                  tATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGa                  >  1:1815537/1‑150 (MQ=255)
                                                                                                                                       aGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGt             >  1:1799402/1‑150 (MQ=255)
                                                                                                                                       aGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGt             >  2:614566/1‑150 (MQ=255)
                                                                                                                                        gCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTg            <  2:53500/150‑1 (MQ=255)
                                                                                                                                         cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa           <  1:2836650/150‑1 (MQ=255)
                                                                                                                                         cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa           <  2:1437486/150‑1 (MQ=255)
                                                                                                                                         cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGAGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa           <  1:2864817/150‑1 (MQ=255)
                                                                                                                                          aCTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAc          <  1:70221/150‑1 (MQ=255)
                                                                                                                                           cTTGGATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACt         <  1:727722/150‑1 (MQ=255)
                                                                                                                                            ttGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTc        >  2:4564020/1‑150 (MQ=255)
                                                                                                                                               cATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc     >  2:3198424/1‑150 (MQ=255)
                                                                                                                                                  ggCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAtt  >  2:3952626/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  CP009273/618900‑619195

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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  CP009273/618902‑619195
                                                                                                                                                  |                                                                                                                                                   
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGATTTGCGACATTTTAGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAATCCTATCCAGCACTTGCATGGAGCT                                                                                                                                                  >  SRR19973073.397329/1‑150 (MQ=60)
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAATCCTATCCAGCACTTGCATGGAGCT                                                                                                                                                  >  SRR19973073.397347/1‑150 (MQ=60)
 GCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTG                                                                                                                                                 >  SRR19973073.53656/1‑150 (MQ=60)
 GCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTG                                                                                                                                                 >  SRR19973073.1287919/1‑150 (MQ=60)
  CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGT                                                                                                                                                <  SRR19973073.1652172/150‑1 (MQ=60)
  CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGT                                                                                                                                                <  SRR19973073.4023573/150‑1 (MQ=60)
      TTTGGTTTGTTGTTTATCGGCGTGTTGTTTTTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTA                                                                                                                                            <  SRR19973073.4475788/150‑1 (MQ=60)
      TTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTA                                                                                                                                            >  SRR19973073.3405338/1‑150 (MQ=60)
       TTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAA                                                                                                                                           >  SRR19973073.1001885/1‑150 (MQ=60)
       TTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAA                                                                                                                                           >  SRR19973073.2309655/1‑150 (MQ=60)
           TTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCA                                                                                                                                       <  SRR19973073.2368885/150‑1 (MQ=60)
              GTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGAC                                                                                                                                    >  SRR19973073.1370327/1‑150 (MQ=60)
                        GGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTC                                                                                                                          <  SRR19973073.340400/150‑1 (MQ=60)
                        GGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTC                                                                                                                          <  SRR19973073.340385/150‑1 (MQ=60)
                         GCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCG                                                                                                                         >  SRR19973073.2622252/1‑150 (MQ=60)
                         GCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCG                                                                                                                         >  SRR19973073.2622259/1‑150 (MQ=60)
                         GCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCG                                                                                                                         >  SRR19973073.2713970/1‑150 (MQ=60)
                          CGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGG                                                                                                                        >  SRR19973073.3792114/1‑150 (MQ=60)
                          CGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGG                                                                                                                        >  SRR19973073.3792138/1‑150 (MQ=60)
                          CGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGG                                                                                                                        >  SRR19973073.4111009/1‑150 (MQ=60)
                           GTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGT                                                                                                                       >  SRR19973073.1939709/1‑150 (MQ=60)
                           GTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGT                                                                                                                       >  SRR19973073.4271308/1‑150 (MQ=60)
                            TGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTT                                                                                                                      <  SRR19973073.4154746/150‑1 (MQ=60)
                              TTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCC                                                                                                                    >  SRR19973073.4333895/1‑150 (MQ=60)
                               TGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCC                                                                                                                   <  SRR19973073.4592821/150‑1 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                  <  SRR19973073.4645416/150‑1 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                  <  SRR19973073.960340/150‑1 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                  <  SRR19973073.3553735/150‑1 (MQ=60)
                                GTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCA                                                                                                                  <  SRR19973073.365557/150‑1 (MQ=60)
                                 TTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAG                                                                                                                 <  SRR19973073.1916463/150‑1 (MQ=60)
                                       CTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATA                                                                                                           <  SRR19973073.3678522/150‑1 (MQ=60)
                                       CTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATA                                                                                                           <  SRR19973073.3678536/150‑1 (MQ=60)
                                        TGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAA                                                                                                          <  SRR19973073.4128305/150‑1 (MQ=60)
                                           TGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACC                                                                                                       >  SRR19973073.544429/1‑150 (MQ=60)
                                            GCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCT                                                                                                      <  SRR19973073.3789052/150‑1 (MQ=60)
                                                 TGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTG                                                                                                 >  SRR19973073.2768335/1‑150 (MQ=60)
                                                  GGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGT                                                                                                <  SRR19973073.4195519/150‑1 (MQ=60)
                                                  GGAGTTGCGACATTTTCGACAGACGCCGCCGCAGGTGACAGCGTCAGACAGTTAATGCTTAAATCATCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGT                                                                                                <  SRR19973073.4195574/150‑1 (MQ=60)
                                                     GTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTT                                                                                             <  SRR19973073.988919/150‑1 (MQ=60)
                                                     GTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTT                                                                                             >  SRR19973073.1159107/1‑150 (MQ=60)
                                                     GTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTGGTTTT                                                                                             >  SRR19973073.1159077/1‑150 (MQ=60)
                                                      TTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTG                                                                                            <  SRR19973073.1205642/150‑1 (MQ=60)
                                                       TGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGT                                                                                           <  SRR19973073.1098615/150‑1 (MQ=60)
                                                       TGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGT                                                                                           <  SRR19973073.4545011/150‑1 (MQ=60)
                                                       TGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGT                                                                                           <  SRR19973073.4545525/150‑1 (MQ=60)
                                                          GACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACT                                                                                        >  SRR19973073.2956358/1‑150 (MQ=60)
                                                           ACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTG                                                                                       >  SRR19973073.3506679/1‑150 (MQ=60)
                                                           ACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTG                                                                                       >  SRR19973073.3506690/1‑150 (MQ=60)
                                                           ACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTG                                                                                       >  SRR19973073.3506779/1‑150 (MQ=60)
                                                            CATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGC                                                                                      >  SRR19973073.1441496/1‑150 (MQ=60)
                                                            CATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGC                                                                                      >  SRR19973073.1441463/1‑150 (MQ=60)
                                                              TTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCNGCGCATAAACCTGCTTGTTTTGTACTGCAG                                                                                    <  SRR19973073.4291806/150‑1 (MQ=60)
                                                              TTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAG                                                                                    <  SRR19973073.4185865/150‑1 (MQ=60)
                                                              TTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAG                                                                                    <  SRR19973073.4185861/150‑1 (MQ=60)
                                                              TTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAG                                                                                    <  SRR19973073.4185828/150‑1 (MQ=60)
                                                                  CGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAG                                                                                >  SRR19973073.4344827/1‑150 (MQ=60)
                                                                   GCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGG                                                                               <  SRR19973073.2501298/150‑1 (MQ=60)
                                                                     CAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTG                                                                             <  SRR19973073.1144482/150‑1 (MQ=60)
                                                                         CGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCG                                                                         >  SRR19973073.2129821/1‑150 (MQ=60)
                                                                         CGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCGTGTTTTGTACTGCAGGCAGGTGCGCG                                                                         >  SRR19973073.2129936/1‑150 (MQ=60)
                                                                          GCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGA                                                                        <  SRR19973073.246084/150‑1 (MQ=60)
                                                                          GCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGA                                                                        >  SRR19973073.86250/1‑150 (MQ=60)
                                                                           CCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAG                                                                       <  SRR19973073.849063/150‑1 (MQ=60)
                                                                           CCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAG                                                                       <  SRR19973073.3171490/150‑1 (MQ=60)
                                                                             GCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCA                                                                     >  SRR19973073.3245507/1‑150 (MQ=60)
                                                                             GCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCA                                                                     <  SRR19973073.1047701/150‑1 (MQ=60)
                                                                             GCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCA                                                                     <  SRR19973073.1047725/150‑1 (MQ=60)
                                                                              CCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAG                                                                    <  SRR19973073.2768186/150‑1 (MQ=60)
                                                                               CGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGC                                                                   <  SRR19973073.70238/150‑1 (MQ=60)
                                                                               CGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGC                                                                   <  SRR19973073.70223/150‑1 (MQ=60)
                                                                               CGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTCTTTTGTACTGCAGGCAGGTGCGCGAGCAGC                                                                   <  SRR19973073.2385362/150‑1 (MQ=60)
                                                                                 CAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGG                                                                 <  SRR19973073.2047510/150‑1 (MQ=60)
                                                                                    GTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATT                                                              <  SRR19973073.749887/150‑1 (MQ=60)
                                                                                    GTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATT                                                              <  SRR19973073.4626918/150‑1 (MQ=60)
                                                                                    GTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATT                                                              <  SRR19973073.1999823/150‑1 (MQ=60)
                                                                                     TGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTA                                                             >  SRR19973073.2114332/1‑150 (MQ=60)
                                                                                     TGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTA                                                             >  SRR19973073.3161506/1‑150 (MQ=60)
                                                                                      GACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAG                                                            >  SRR19973073.2965089/1‑150 (MQ=60)
                                                                                        CAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCA                                                          >  SRR19973073.4375738/1‑150 (MQ=60)
                                                                                         AGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCAT                                                         >  SRR19973073.1782411/1‑150 (MQ=60)
                                                                                         AGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCAT                                                         >  SRR19973073.1721360/1‑150 (MQ=60)
                                                                                         AGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCAT                                                         >  SRR19973073.1782408/1‑150 (MQ=60)
                                                                                          GCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATA                                                        <  SRR19973073.2957783/150‑1 (MQ=60)
                                                                                           CGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAA                                                       <  SRR19973073.4306641/150‑1 (MQ=60)
                                                                                             TCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAAT                                                     >  SRR19973073.2843352/1‑150 (MQ=60)
                                                                                              CCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATA                                                    <  SRR19973073.2450702/150‑1 (MQ=60)
                                                                                              CCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATA                                                    <  SRR19973073.3983569/150‑1 (MQ=60)
                                                                                               CGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAG                                                   >  SRR19973073.2606631/1‑150 (MQ=60)
                                                                                                GACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGC                                                  >  SRR19973073.3068274/1‑150 (MQ=60)
                                                                                                GACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGC                                                  >  SRR19973073.3068307/1‑150 (MQ=60)
                                                                                                  CAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAT                                                <  SRR19973073.1451825/150‑1 (MQ=60)
                                                                                                  CAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAT                                                <  SRR19973073.1411783/150‑1 (MQ=60)
                                                                                                   AGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATC                                               >  SRR19973073.3220194/1‑150 (MQ=60)
                                                                                                    GTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCG                                              <  SRR19973073.3537210/150‑1 (MQ=60)
                                                                                                             TAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTC                                     <  SRR19973073.2863619/150‑1 (MQ=60)
                                                                                                              AAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCT                                    <  SRR19973073.3803626/150‑1 (MQ=60)
                                                                                                               AAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTG                                   >  SRR19973073.680696/1‑150 (MQ=60)
                                                                                                               AAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTG                                   >  SRR19973073.2999190/1‑150 (MQ=60)
                                                                                                                 ACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGAT                                 <  SRR19973073.2515053/150‑1 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                >  SRR19973073.2138116/1‑150 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                <  SRR19973073.3452344/150‑1 (MQ=60)
                                                                                                                  CAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATC                                <  SRR19973073.4037220/150‑1 (MQ=60)
                                                                                                                         TTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCG                         <  SRR19973073.1225682/150‑1 (MQ=60)
                                                                                                                          TAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGA                        >  SRR19973073.2974396/1‑150 (MQ=60)
                                                                                                                           AAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCTCGAGCAGCGGATTAGCATAAATAGCATCTGCATCTTTCTGATCACCGGCGAA                       >  SRR19973073.706888/1‑150 (MQ=60)
                                                                                                                            AGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAC                      >  SRR19973073.4586021/1‑150 (MQ=60)
                                                                                                                            AGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAC                      <  SRR19973073.2521278/150‑1 (MQ=60)
                                                                                                                               CTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGG                   <  SRR19973073.741054/150‑1 (MQ=60)
                                                                                                                                TATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGA                  >  SRR19973073.1819292/1‑150 (MQ=60)
                                                                                                                                TATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGA                  >  SRR19973073.1819173/1‑150 (MQ=60)
                                                                                                                                 ATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAA                 >  SRR19973073.2026074/1‑150 (MQ=60)
                                                                                                                                 ATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAA                 >  SRR19973073.2026083/1‑150 (MQ=60)
                                                                                                                                   CCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATA               >  SRR19973073.242901/1‑150 (MQ=60)
                                                                                                                                    CAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAG              <  SRR19973073.2597321/150‑1 (MQ=60)
                                                                                                                                     AGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCTTCGGCATCTTTCTGATCACCGGCGAACAGGAATAGT             >  SRR19973073.213842/1‑150 (MQ=60)
                                                                                                                                     AGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGT             >  SRR19973073.1803118/1‑150 (MQ=60)
                                                                                                                                      GCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTG            >  SRR19973073.3504684/1‑150 (MQ=60)
                                                                                                                                      GCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGAAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTG            >  SRR19973073.3505131/1‑150 (MQ=60)
                                                                                                                                       CACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGA           <  SRR19973073.2842468/150‑1 (MQ=60)
                                                                                                                                       CACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGAGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGA           <  SRR19973073.2870676/150‑1 (MQ=60)
                                                                                                                                        ACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAC          <  SRR19973073.70412/150‑1 (MQ=60)
                                                                                                                                        ACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAC          <  SRR19973073.70689/150‑1 (MQ=60)
                                                                                                                                        ACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAC          >  SRR19973073.506274/1‑150 (MQ=60)
                                                                                                                                         CTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACT         <  SRR19973073.729225/150‑1 (MQ=60)
                                                                                                                                           TGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCG       <  SRR19973073.4512986/150‑1 (MQ=60)
                                                                                                                                              ATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCA    >  SRR19973073.531606/1‑150 (MQ=60)
                                                                                                                                              ATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCA    <  SRR19973073.1540853/150‑1 (MQ=60)
                                                                                                                                              ATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCA    >  SRR19973073.3261/1‑150 (MQ=60)
                                                                                                                                                GGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  SRR19973073.2847522/1‑150 (MQ=60)
                                                                                                                                                GGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  <  SRR19973073.1944292/150‑1 (MQ=60)
                                                                                                                                                GGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  SRR19973073.2847685/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                   
AGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATT  >  CP009273/618902‑619195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: