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 F1 I1 R1 744 338.8 17650688 97.8% 17262372 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,349 T→C P77P (CCA→CCG)  mobA ← MobA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,3490TC99.0% 973.7 / ‑2.8 287P77P (CCA→CCG) mobAMobA
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (149/135);  total (151/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 297 total aligned reads displayed.

AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAG  >  NC_000913/4039259‑4039441
                                                                                          |                                                                                            
aGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCgg                                                                                            <  2:7036676/93‑1 (MQ=255)
 gggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCggg                                                                                           >  1:6658266/1‑93 (MQ=255)
 gggggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCggg                                                                                           <  1:3750909/93‑1 (MQ=255)
    ggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTaa                                                                                        <  1:7461401/93‑1 (MQ=255)
    ggAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTaa                                                                                        >  1:2325045/1‑93 (MQ=255)
     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  1:4935349/93‑1 (MQ=255)
     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       >  2:3467242/1‑93 (MQ=255)
     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  1:2850432/93‑1 (MQ=255)
     gAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAt                                                                                       <  2:5838260/93‑1 (MQ=255)
       aTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATcc                                                                                     <  2:8634840/93‑1 (MQ=255)
       aTGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATcc                                                                                     <  1:682459/93‑1 (MQ=255)
         gTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGc                                                                                   <  1:264549/93‑1 (MQ=255)
            aGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAg                                                                                >  2:7149084/1‑93 (MQ=255)
             ggCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGt                                                                               >  1:5355558/1‑93 (MQ=255)
              gCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTg                                                                              >  2:2633919/1‑93 (MQ=255)
               cGTATCGCGCGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  1:8729962/93‑1 (MQ=255)
               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  2:781434/93‑1 (MQ=255)
               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             >  1:8348982/1‑93 (MQ=255)
               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  2:6599763/93‑1 (MQ=255)
               cGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGa                                                                             <  1:3968036/93‑1 (MQ=255)
                 tATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCTCCAGTGAAt                                                                           >  2:7885109/1‑93 (MQ=255)
                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCt                                                                         <  1:4335667/93‑1 (MQ=255)
                   tCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCt                                                                         <  2:8210490/93‑1 (MQ=255)
                    cGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCtt                                                                        <  1:7745590/93‑1 (MQ=255)
                         gggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTACTCCGACCGTGAATCTTCAATc                                                                   >  2:2958365/1‑93 (MQ=255)
                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  >  2:1281015/1‑93 (MQ=255)
                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  <  1:5177878/93‑1 (MQ=255)
                          ggCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCa                                                                  >  1:3814012/1‑93 (MQ=255)
                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                 <  1:1326492/93‑1 (MQ=255)
                           gCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCAc                                                                 >  2:2391055/1‑93 (MQ=255)
                            caaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACt                                                                >  1:1103042/1‑93 (MQ=255)
                            caaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACt                                                                <  1:7853640/93‑1 (MQ=255)
                             aaaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACtt                                                               <  1:6482777/93‑1 (MQ=255)
                              aaacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACttt                                                              >  1:7766206/1‑93 (MQ=255)
                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             >  2:3708594/1‑93 (MQ=255)
                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             >  2:3541560/1‑93 (MQ=255)
                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             <  2:3441342/93‑1 (MQ=255)
                               aacaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTc                                                             <  1:4549884/93‑1 (MQ=255)
                                acaaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCa                                                            >  1:5877426/1‑93 (MQ=255)
                                 caaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAg                                                           >  1:7057383/1‑93 (MQ=255)
                                 caaaaACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAg                                                           >  1:5812649/1‑93 (MQ=255)
                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCCGAGGGCCCGGGTAATCCGCCAGTGCATCTTCAATCCCTTTCAgcc                                                         >  2:4314356/1‑91 (MQ=255)
                                   aaaaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAc                                                         >  2:4664861/1‑93 (MQ=255)
                                     aaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCg                                                       >  2:6345071/1‑93 (MQ=255)
                                     aaCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCg                                                       <  1:6921858/93‑1 (MQ=255)
                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                      <  2:7283322/93‑1 (MQ=255)
                                      aCCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGc                                                      <  1:6246786/93‑1 (MQ=255)
                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                    >  2:1443404/1‑93 (MQ=255)
                                        cACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCtt                                                    <  2:5099433/93‑1 (MQ=255)
                                           tCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 >  1:8595453/1‑93 (MQ=255)
                                           ccACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 <  1:6581283/92‑1 (MQ=255)
                                           ccACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGcc                                                 <  2:2092672/92‑1 (MQ=255)
                                            cACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCt                                                >  1:8165343/1‑93 (MQ=255)
                                            cACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCt                                                >  1:2521956/1‑93 (MQ=255)
                                             aCCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTg                                               >  1:8635700/1‑93 (MQ=255)
                                              ccAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGa                                              >  2:2789230/1‑93 (MQ=255)
                                                aGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATa                                            >  1:4832402/1‑93 (MQ=255)
                                                 gCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAg                                           <  2:2009661/93‑1 (MQ=255)
                                                 gCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAg                                           >  2:3961820/1‑93 (MQ=255)
                                                     ccTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGAttt                                       >  1:5493511/1‑93 (MQ=255)
                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      >  2:6966333/1‑93 (MQ=255)
                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      >  1:6434502/1‑93 (MQ=255)
                                                      ctgctgCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTc                                      <  1:7174667/93‑1 (MQ=255)
                                                        gctgcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCt                                    <  2:3032950/93‑1 (MQ=255)
                                                           gcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgat                                 >  2:1658567/1‑93 (MQ=255)
                                                           gcATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgat                                 >  2:3599616/1‑93 (MQ=255)
                                                             aTTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCtgatga                               <  2:1510950/93‑1 (MQ=255)
                                                             aTTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTACCtgatga                               >  2:8380482/1‑93 (MQ=255)
                                                               tACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACg                             >  1:6335753/1‑93 (MQ=255)
                                                                 cTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAt                           <  1:670030/93‑1 (MQ=255)
                                                                  tGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAtt                          <  2:4942547/93‑1 (MQ=255)
                                                                   gAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTa                         >  1:3172266/1‑93 (MQ=255)
                                                                    aaaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                        >  2:8026406/1‑93 (MQ=255)
                                                                    aaaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAg                        <  2:7496775/93‑1 (MQ=255)
                                                                     aaGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGc                       >  1:4648564/1‑93 (MQ=255)
                                                                      aGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCa                      <  1:4088627/93‑1 (MQ=255)
                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     <  2:770257/93‑1 (MQ=255)
                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     <  2:557636/93‑1 (MQ=255)
                                                                       gCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAt                     >  2:5780051/1‑93 (MQ=255)
                                                                        cATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAtt                    <  1:7500826/93‑1 (MQ=255)
                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:4525059/1‑93 (MQ=255)
                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:8229261/1‑93 (MQ=255)
                                                                          ttCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTaa                  >  2:3812583/1‑93 (MQ=255)
                                                                            ccTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAcc                >  2:5283520/1‑93 (MQ=255)
                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa               >  1:3057080/1‑93 (MQ=255)
                                                                             cTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCa               <  2:4474465/93‑1 (MQ=255)
                                                                              tgtgagagGGCCCGGGTAAACCGCCAGTGAATCTTCAATCCCTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              <  2:2224180/89‑1 (MQ=255)
                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              >  2:6754833/1‑93 (MQ=255)
                                                                              tGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCac              <  1:806620/93‑1 (MQ=255)
                                                                               gCCAGAGGGCCCGGGTAATACGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacg             >  1:5946344/1‑93 (MQ=255)
                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           <  2:181864/93‑1 (MQ=255)
                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           >  2:1022578/1‑93 (MQ=255)
                                                                                 cAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgac           >  2:4698741/1‑93 (MQ=255)
                                                                                  agagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg          >  1:1792776/1‑93 (MQ=255)
                                                                                  agagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCacgacg          >  1:1693084/1‑93 (MQ=255)
                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt         >  2:6175482/1‑93 (MQ=255)
                                                                                   gagGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGt         >  1:6515110/1‑93 (MQ=255)
                                                                                    agGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTg        <  2:1547409/93‑1 (MQ=255)
                                                                                     gggCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTga       <  1:1502206/93‑1 (MQ=255)
                                                                                          tGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGagagag  <  1:1462529/93‑1 (MQ=255)
                                                                                          |                                                                                            
AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAG  >  NC_000913/4039259‑4039441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAGC  >  NC_000913/4039256‑4039442
                                                                                             |                                                                                             
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCC                                                                                               <  SRR4302541.6639253/94‑1 (MQ=60)
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCC                                                                                               <  SRR4302541.4702462/94‑1 (MQ=60)
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCC                                                                                               >  SRR4302541.3489311/1‑94 (MQ=60)
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCC                                                                                               <  SRR4302541.6378911/94‑1 (MQ=60)
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCC                                                                                               <  SRR4302541.1615094/94‑1 (MQ=60)
  CAGGGGGAATGTAATGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGG                                                                                             <  SRR4302541.8749869/94‑1 (MQ=60)
  CAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGG                                                                                             <  SRR4302541.621883/94‑1 (MQ=60)
  CAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGG                                                                                             >  SRR4302541.1509199/1‑94 (MQ=60)
  CAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGG                                                                                             >  SRR4302541.2773543/1‑94 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            >  SRR4302541.6726573/1‑94 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            <  SRR4302541.3045460/94‑1 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            <  SRR4302541.8187356/94‑1 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            >  SRR4302541.8427947/1‑94 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            >  SRR4302541.1756137/1‑94 (MQ=60)
   AGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGG                                                                                            <  SRR4302541.1709388/94‑1 (MQ=60)
    GGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGT                                                                                           <  SRR4302541.3788645/94‑1 (MQ=60)
     GGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTA                                                                                          <  SRR4302541.7834967/94‑1 (MQ=60)
     GGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTA                                                                                          >  SRR4302541.5537582/1‑94 (MQ=60)
     GGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTA                                                                                          >  SRR4302541.5874979/1‑94 (MQ=60)
      GGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAA                                                                                         >  SRR4302541.2347605/1‑94 (MQ=60)
       GGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAT                                                                                        <  SRR4302541.5910502/94‑1 (MQ=60)
       GGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAAT                                                                                        <  SRR4302541.7541115/94‑1 (MQ=60)
        GAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATC                                                                                       <  SRR4302541.4985719/94‑1 (MQ=60)
        GAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATC                                                                                       <  SRR4302541.2879098/94‑1 (MQ=60)
          ATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCG                                                                                     >  SRR4302541.5828730/1‑94 (MQ=60)
          ATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCG                                                                                     >  SRR4302541.6551397/1‑94 (MQ=60)
          ATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCG                                                                                     <  SRR4302541.707940/94‑1 (MQ=60)
          ATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCG                                                                                     <  SRR4302541.688896/94‑1 (MQ=60)
           TGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGC                                                                                    <  SRR4302541.3927487/94‑1 (MQ=60)
            GTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCC                                                                                   <  SRR4302541.266720/94‑1 (MQ=60)
              AAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAG                                                                                 >  SRR4302541.8509519/1‑94 (MQ=60)
              AAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAG                                                                                 >  SRR4302541.6568768/1‑94 (MQ=60)
               AGGCGTATCGCGCGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGT                                                                                >  SRR4302541.6481874/1‑94 (MQ=60)
               AGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGT                                                                                >  SRR4302541.5410580/1‑94 (MQ=60)
               AGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGT                                                                                <  SRR4302541.6344084/94‑1 (MQ=60)
                GGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTG                                                                               >  SRR4302541.7314050/1‑94 (MQ=60)
                 GCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGA                                                                              >  SRR4302541.8439602/1‑94 (MQ=60)
                  CGTATCGCGCGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAA                                                                             <  SRR4302541.8823435/94‑1 (MQ=60)
                  CGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAA                                                                             <  SRR4302541.4007430/94‑1 (MQ=60)
                   GTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAAT                                                                            <  SRR4302541.6816879/94‑1 (MQ=60)
                     ATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCT                                                                          <  SRR4302541.4617749/94‑1 (MQ=60)
                      TCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTT                                                                         <  SRR4302541.4378961/94‑1 (MQ=60)
                       CGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTC                                                                        >  SRR4302541.8155885/1‑94 (MQ=60)
                       CGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTC                                                                        >  SRR4302541.7668465/1‑94 (MQ=60)
                       CGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTC                                                                        <  SRR4302541.7828008/94‑1 (MQ=60)
                          ACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAAT                                                                     <  SRR4302541.7569203/94‑1 (MQ=60)
                            GGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCTGGTAATCCGCCAGTGAATCTTCAATCA                                                                   >  SRR4302541.2893611/1‑94 (MQ=60)
                            GGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCA                                                                   >  SRR4302541.3852430/1‑94 (MQ=60)
                             GGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCAC                                                                  <  SRR4302541.5231144/94‑1 (MQ=60)
                              GCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACT                                                                 <  SRR4302541.1339536/94‑1 (MQ=60)
                              GCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACT                                                                 >  SRR4302541.1113374/1‑94 (MQ=60)
                              GCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACT                                                                 >  SRR4302541.6033129/1‑94 (MQ=60)
                               CAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTT                                                                <  SRR4302541.7939854/94‑1 (MQ=60)
                               CAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTT                                                                <  SRR4302541.5968924/94‑1 (MQ=60)
                               CAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTT                                                                <  SRR4302541.4969136/94‑1 (MQ=60)
                                AAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTT                                                               <  SRR4302541.6548852/94‑1 (MQ=60)
                                AAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTT                                                               <  SRR4302541.6388599/94‑1 (MQ=60)
                                AAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTT                                                               >  SRR4302541.1502702/1‑94 (MQ=60)
                                AAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTT                                                               >  SRR4302541.7848826/1‑94 (MQ=60)
                                 AAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTC                                                              <  SRR4302541.8151154/94‑1 (MQ=60)
                                 AAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTC                                                              <  SRR4302541.6861728/94‑1 (MQ=60)
                                  AACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCA                                                             >  SRR4302541.5937238/1‑94 (MQ=60)
                                  AACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCA                                                             <  SRR4302541.8764158/94‑1 (MQ=60)
                                  AACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCA                                                             >  SRR4302541.5937239/1‑94 (MQ=60)
                                  AACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCA                                                             <  SRR4302541.4595908/94‑1 (MQ=60)
                                   ACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAG                                                            >  SRR4302541.7132437/1‑94 (MQ=60)
                                   ACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAG                                                            <  SRR4302541.7046638/94‑1 (MQ=60)
                                   ACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAG                                                            <  SRR4302541.5430827/94‑1 (MQ=60)
                                   ACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAG                                                            >  SRR4302541.5871711/1‑94 (MQ=60)
                                    CAAAAACCACTCACCAGCTTCCTGCTGCATGACCGAAAGCATTCCTGCCAGAGGGCCCGGGGAATCCGCCAGTGAATCGTCAATCACTTTCAGA                                                           >  SRR4302541.7929274/1‑94 (MQ=60)
                                       AAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCG                                                        <  SRR4302541.478309/94‑1 (MQ=60)
                                        AACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGC                                                       <  SRR4302541.6995723/94‑1 (MQ=60)
                                         ACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCT                                                      <  SRR4302541.7173963/94‑1 (MQ=60)
                                         ACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCT                                                      <  SRR4302541.6310604/94‑1 (MQ=60)
                                          CCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTT                                                     >  SRR4302541.239471/1‑94 (MQ=60)
                                          CCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTT                                                     <  SRR4302541.869450/94‑1 (MQ=60)
                                          CCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTT                                                     <  SRR4302541.3906037/94‑1 (MQ=60)
                                          CCACTCACCAGCTTCCTGCTGCATTACTGAAAGCACTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGACTCTTCAATCACTTTCACACGGCTT                                                     >  SRR4302541.1380484/1‑94 (MQ=60)
                                          CCACTCACCAGCTTCCTACTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTT                                                     <  SRR4302541.7895029/94‑1 (MQ=60)
                                           CACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTG                                                    >  SRR4302541.1404459/1‑94 (MQ=60)
                                             cccccccgccTCCTGCGGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCC                                                  <  SRR4302541.6728762/84‑1 (MQ=60)
                                             CTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCC                                                  >  SRR4302541.8687670/1‑94 (MQ=60)
                                             CTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCC                                                  <  SRR4302541.6728499/94‑1 (MQ=60)
                                              TCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCT                                                 >  SRR4302541.8254528/1‑94 (MQ=60)
                                              TCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCT                                                 >  SRR4302541.2547330/1‑94 (MQ=60)
                                              CCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCT                                                 <  SRR4302541.6648637/94‑1 (MQ=60)
                                              CCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCT                                                 <  SRR4302541.1301872/94‑1 (MQ=60)
                                               CACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGAGCGCTTGCCTG                                                >  SRR4302541.8728220/1‑94 (MQ=60)
                                               CACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTG                                                >  SRR4302541.8728221/1‑94 (MQ=60)
                                                 CCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGAT                                              >  SRR4302541.6439603/1‑94 (MQ=60)
                                                  CAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             >  SRR4302541.2798903/1‑94 (MQ=60)
                                                  CAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             <  SRR4302541.6924876/94‑1 (MQ=60)
                                                  CAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             >  SRR4302541.6704620/1‑94 (MQ=60)
                                                  CAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             >  SRR4302541.4881662/1‑94 (MQ=60)
                                                  CAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             >  SRR4302541.4881663/1‑94 (MQ=60)
                                                  CAGCTTCCGGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATA                                             <  SRR4302541.5374976/94‑1 (MQ=60)
                                                   AGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCATACCGCTTGCCTGATAG                                            >  SRR4302541.7836833/1‑94 (MQ=60)
                                                   AGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAG                                            >  SRR4302541.8764851/1‑94 (MQ=60)
                                                      TTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATT                                         >  SRR4302541.3693820/1‑94 (MQ=60)
                                                      TTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATT                                         <  SRR4302541.4149879/94‑1 (MQ=60)
                                                       TCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTT                                        >  SRR4302541.5549796/1‑94 (MQ=60)
                                                        CCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTC                                       >  SRR4302541.6500085/1‑94 (MQ=60)
                                                         CTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCC                                      <  SRR4302541.7251175/94‑1 (MQ=60)
                                                         CTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCC                                      >  SRR4302541.2586789/1‑94 (MQ=60)
                                                          TGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCTGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCT                                     <  SRR4302541.6593924/94‑1 (MQ=60)
                                                          TGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCT                                     <  SRR4302541.6510291/94‑1 (MQ=60)
                                                          CGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCT                                     <  SRR4302541.6593925/94‑1 (MQ=60)
                                                            CTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGA                                   >  SRR4302541.6901873/1‑94 (MQ=60)
                                                            CTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGA                                   >  SRR4302541.5010757/1‑94 (MQ=60)
                                                            CTGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGA                                   >  SRR4302541.2096157/1‑94 (MQ=60)
                                                             TGCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGAT                                  <  SRR4302541.5117952/94‑1 (MQ=60)
                                                              GCATTACTGAAAGCATTGCTGGCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGTTTGCCTGATAGATTTCCTGATG                                 <  SRR4302541.7686384/94‑1 (MQ=60)
                                                              GCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATG                                 >  SRR4302541.4453707/1‑94 (MQ=60)
                                                              GCATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATG                                 >  SRR4302541.8096120/1‑94 (MQ=60)
                                                               CATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGA                                >  SRR4302541.4566700/1‑94 (MQ=60)
                                                               CATTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGA                                >  SRR4302541.7425273/1‑94 (MQ=60)
                                                                GTTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGAC                               <  SRR4302541.1947791/94‑1 (MQ=60)
                                                                 TTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACG                              <  SRR4302541.3124861/94‑1 (MQ=60)
                                                                 TTACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACG                              >  SRR4302541.6400523/1‑94 (MQ=60)
                                                                   ACTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGAT                            >  SRR4302541.8374470/1‑94 (MQ=60)
                                                                    CTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATT                           <  SRR4302541.676351/94‑1 (MQ=60)
                                                                    CTGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATT                           >  SRR4302541.1655604/1‑94 (MQ=60)
                                                                     TGAAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTA                          >  SRR4302541.3204321/1‑94 (MQ=60)
                                                                       AAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGC                        <  SRR4302541.2387057/94‑1 (MQ=60)
                                                                       AAAGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGC                        >  SRR4302541.4695902/1‑94 (MQ=60)
                                                                         AGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAT                      <  SRR4302541.4128996/94‑1 (MQ=60)
                                                                         AGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAT                      <  SRR4302541.664617/94‑1 (MQ=60)
                                                                         AGCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCAT                      >  SRR4302541.3272139/1‑94 (MQ=60)
                                                                          GCATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATT                     <  SRR4302541.714480/94‑1 (MQ=60)
                                                                           CATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTA                    <  SRR4302541.7581010/94‑1 (MQ=60)
                                                                           CATTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTA                    <  SRR4302541.7581011/94‑1 (MQ=60)
                                                                             TTCCTGCCAGAGGGGCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAC                  <  SRR4302541.6177642/94‑1 (MQ=60)
                                                                             TTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAC                  <  SRR4302541.3818131/94‑1 (MQ=60)
                                                                             TTCCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAAC                  <  SRR4302541.8883324/94‑1 (MQ=60)
                                                                               CCTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCA                >  SRR4302541.3088048/1‑94 (MQ=60)
                                                                                CTGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCAC               <  SRR4302541.4252224/94‑1 (MQ=60)
                                                                                 TGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACG              >  SRR4302541.2650427/1‑94 (MQ=60)
                                                                                 TGCCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACG              <  SRR4302541.814022/94‑1 (MQ=60)
                                                                                 TGCCAGAGGGCCCGGGTAATACGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACG              >  SRR4302541.6006937/1‑94 (MQ=60)
                                                                                   CCAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGAC            >  SRR4302541.8538731/1‑94 (MQ=60)
                                                                                    CAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACG           >  SRR4302541.3837754/1‑94 (MQ=60)
                                                                                    CAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACG           >  SRR4302541.1709702/1‑94 (MQ=60)
                                                                                    CAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACG           <  SRR4302541.1713164/94‑1 (MQ=60)
                                                                                    CAGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACG           >  SRR4302541.1810186/1‑94 (MQ=60)
                                                                                     AGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGT          <  SRR4302541.4853131/94‑1 (MQ=60)
                                                                                     AGAGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGT          >  SRR4302541.6581592/1‑94 (MQ=60)
                                                                                       AGGGCGCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGA        <  SRR4302541.4638872/94‑1 (MQ=60)
                                                                                       AGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGA        <  SRR4302541.4073200/94‑1 (MQ=60)
                                                                                       AGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGA        >  SRR4302541.3123484/1‑94 (MQ=60)
                                                                                       AGGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGA        >  SRR4302541.6529384/1‑94 (MQ=60)
                                                                                        GGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAG       >  SRR4302541.40462/1‑94 (MQ=60)
                                                                                        GGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAG       <  SRR4302541.5729051/94‑1 (MQ=60)
                                                                                        GGGCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAG       <  SRR4302541.1517080/94‑1 (MQ=60)
                                                                                          GCCCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTGGCCTGATAGATTTCCTGATGACGATTAGCATTATCCACGACGTGAGAG     <  SRR4302541.5246187/94‑1 (MQ=60)
                                                                                            CCGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAG   <  SRR4302541.5854979/94‑1 (MQ=60)
                                                                                             TGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAGC  <  SRR4302541.1476839/94‑1 (MQ=60)
                                                                                             CGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAGC  >  SRR4302541.5873806/1‑94 (MQ=60)
                                                                                             |                                                                                             
ATCAGGGGGAATGTAAGGCGTATCGCACGGGCAAAACAAAAACCACTCACCAGCTTCCTGCTGCATTACTGAAAGCATTCCTGCCAGAGGGCCTGGGTAATCCGCCAGTGAATCTTCAATCACTTTCAGACCGCTTGCCTGATAGATTTCCTGATGACGATTAGCATTAACCACGACGTGAGAGAGC  >  NC_000913/4039256‑4039442

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑