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 1548 369.5 20095359 96.9% 19472402 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,572,698 G→A R177C (CGC→TGC)  pqqL ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,6980GA99.6% 957.4 / NA 277R177C (CGC→TGC) pqqLputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (128/148);  minor base T (0/1);  total (128/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 284 total aligned reads displayed.

CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787
                                                                                          |                                                                                         
ccGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCAcc                                                                                         >  1:1962138/1‑93 (MQ=255)
  gTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAt                                                                                       <  1:5295402/93‑1 (MQ=255)
   tGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAtt                                                                                      <  2:1276325/93‑1 (MQ=255)
    ggCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCAttt                                                                                     <  2:3554156/93‑1 (MQ=255)
        aCAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCa                                                                                 <  2:825181/93‑1 (MQ=255)
          aGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCCTc                                                                               >  2:622293/1‑93 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              >  1:4138834/1‑93 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATcc                                                                              <  2:1520480/93‑1 (MQ=255)
           gTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCATGAGAGGTGCAGCATTTCGCATcc                                                                              >  2:2379684/1‑93 (MQ=255)
            tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCt                                                                             <  2:5096664/93‑1 (MQ=255)
            tATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCt                                                                             >  1:4091143/1‑93 (MQ=255)
             aTCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTg                                                                            <  2:3373465/93‑1 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           <  1:8072220/93‑1 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  1:6352081/1‑93 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           >  2:3661351/1‑93 (MQ=255)
              tCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGa                                                                           <  1:4292231/93‑1 (MQ=255)
               ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          >  2:9679283/1‑93 (MQ=255)
               ccATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAt                                                                          >  2:8444792/1‑93 (MQ=255)
                cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         >  2:982412/1‑93 (MQ=255)
                cATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATg                                                                         <  1:3716382/93‑1 (MQ=255)
                  tCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATggg                                                                       >  2:853549/1‑93 (MQ=255)
                   cAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGc                                                                      <  2:8558828/93‑1 (MQ=255)
                   cAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCACTTCGCATCCTGATGGGc                                                                      >  1:8989909/1‑93 (MQ=255)
                    aGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCa                                                                     >  1:7381807/1‑93 (MQ=255)
                     ggCCGAGAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAc                                                                    <  2:1465680/93‑1 (MQ=255)
                      gCCGATAGGTTCACGGTCTAAATTACGGGAATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGAGCACCATTTCGCATCCTGATGGGCACg                                                                   >  2:3501614/1‑93 (MQ=255)
                         gATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCa                                                                <  2:9305118/93‑1 (MQ=255)
                         gATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCa                                                                >  1:3970818/1‑93 (MQ=255)
                          aTAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCAt                                                               >  2:1633323/1‑93 (MQ=255)
                           tAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGCGGTGCACCATTTCGCATCCTGATGGGCACGCCAtt                                                              >  2:7562260/1‑93 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACTCCATTcc                                                            <  1:6605823/93‑1 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTcc                                                            >  1:2047137/1‑93 (MQ=255)
                             ggTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTcc                                                            >  2:6022669/1‑93 (MQ=255)
                                tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCa                                                         <  1:9129297/93‑1 (MQ=255)
                                tCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGGGCCCCATTTCGCATCCTGATGGGCACGCCATTCCTcg                                                         >  1:2245331/1‑92 (MQ=255)
                                  aCGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGt                                                       <  2:9524127/93‑1 (MQ=255)
                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGGATCCTCATGGGTACGCCATTCCTCAGTa                                                      >  2:7675092/1‑93 (MQ=255)
                                   cGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTa                                                      <  2:9212161/93‑1 (MQ=255)
                                    cgTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTaa                                                     <  1:2766104/92‑1 (MQ=255)
                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    <  2:3059581/93‑1 (MQ=255)
                                     gTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAt                                                    <  1:3659789/93‑1 (MQ=255)
                                      tCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGGACGCCATTCCTCAGTAAtt                                                   <  2:697706/93‑1 (MQ=255)
                                        tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 <  1:5329600/93‑1 (MQ=255)
                                        tAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAc                                                 >  1:5820674/1‑93 (MQ=255)
                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               <  1:8200401/93‑1 (MQ=255)
                                          aaTTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGc                                               <  2:4471931/93‑1 (MQ=255)
                                             tACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAc                                            <  2:8688527/93‑1 (MQ=255)
                                             tACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAc                                            >  2:5847366/1‑93 (MQ=255)
                                               cGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGt                                          >  2:7922080/1‑93 (MQ=255)
                                                 ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTc                                        >  1:2709320/1‑93 (MQ=255)
                                                 ggTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGAACGCCATTCCTCAGTAATTACGCCACGTTc                                        >  2:1191794/1‑93 (MQ=255)
                                                  gTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCa                                       <  1:2986742/93‑1 (MQ=255)
                                                   tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCATCATTTCGCATCCTGATGGCCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      <  2:5666836/93‑1 (MQ=255)
                                                   tATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAg                                      <  2:250032/93‑1 (MQ=255)
                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGc                                     <  1:515020/93‑1 (MQ=255)
                                                    aTTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGc                                     >  2:1916911/1‑93 (MQ=255)
                                                     tttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCg                                    >  1:7660074/1‑93 (MQ=255)
                                                      ttGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGt                                   >  1:2347887/1‑93 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  <  1:2655479/93‑1 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  >  1:7083018/1‑93 (MQ=255)
                                                       tGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTc                                  >  1:7240307/1‑93 (MQ=255)
                                                          cagcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTAc                               >  2:7062931/1‑93 (MQ=255)
                                                           agcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACt                              <  1:9670871/93‑1 (MQ=255)
                                                            gcagGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACtt                             <  1:10030741/93‑1 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          <  2:4090948/93‑1 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          >  2:5247945/1‑93 (MQ=255)
                                                               ggAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGa                          >  1:8365921/1‑93 (MQ=255)
                                                                gAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAg                         <  2:4052208/93‑1 (MQ=255)
                                                                  aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       >  1:3168419/1‑93 (MQ=255)
                                                                  aaGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtt                       <  1:1578769/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      <  2:7412675/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      <  1:177701/93‑1 (MQ=255)
                                                                   aGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGttt                      >  1:5683622/1‑93 (MQ=255)
                                                                    gggCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGtttt                     <  2:4306548/93‑1 (MQ=255)
                                                                       cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  <  1:1283440/93‑1 (MQ=255)
                                                                       cggcgCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCa                  <  1:9305568/93‑1 (MQ=255)
                                                                         gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                >  1:1661239/1‑93 (MQ=255)
                                                                         gcgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACACCACGTTCAGCGTCTACTTCGAGTTTTTCaaa                >  1:1179048/1‑93 (MQ=255)
                                                                          cgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  1:4225317/93‑1 (MQ=255)
                                                                          cgcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAg               <  2:6065460/93‑1 (MQ=255)
                                                                           gcgcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAgg              <  1:3448881/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:6164735/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            >  1:9743702/1‑93 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  2:3112875/93‑1 (MQ=255)
                                                                             gcCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGtt            <  1:4591171/93‑1 (MQ=255)
                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           >  1:3894791/1‑93 (MQ=255)
                                                                              ccTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTg           >  1:4784564/1‑93 (MQ=255)
                                                                               cTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGc          <  2:2260981/93‑1 (MQ=255)
                                                                                tGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCg         >  1:1872573/1‑93 (MQ=255)
                                                                                 gagaGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCgg        <  2:8117917/93‑1 (MQ=255)
                                                                                 gagaGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCgg        >  1:9385724/1‑93 (MQ=255)
                                                                                  agagGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGc       <  1:9424815/93‑1 (MQ=255)
                                                                                     ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    >  1:2537169/1‑93 (MQ=255)
                                                                                     ggTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAtt    <  2:4919672/93‑1 (MQ=255)
                                                                                      gTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:4996144/93‑1 (MQ=255)
                                                                                      gTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTa   <  1:788351/93‑1 (MQ=255)
                                                                                       tGCACCATTTTGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  >  2:6356796/1‑93 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTaa  >  2:2122859/1‑92 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  <  1:4391575/93‑1 (MQ=255)
                                                                                       tGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAc  >  2:6474034/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC  >  NC_000913/1572608‑1572787

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ 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
TGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCCA  >  NC_000913/1572605‑1572791
                                                                                             |                                                                                             
TGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCA                                                                                               >  SRR4302540.778352/1‑94 (MQ=60)
  ACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACC                                                                                             >  SRR4302540.1983455/1‑94 (MQ=60)
  ACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACC                                                                                             <  SRR4302540.3182209/94‑1 (MQ=60)
     GTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATT                                                                                          <  SRR4302540.5353315/94‑1 (MQ=60)
      TGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTT                                                                                         >  SRR4302540.1665531/1‑94 (MQ=60)
         CGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGC                                                                                      <  SRR4302540.6040397/94‑1 (MQ=60)
          GACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCA                                                                                     >  SRR4302540.3765539/1‑94 (MQ=60)
           ACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCAT                                                                                    >  SRR4302540.1878373/1‑94 (MQ=60)
           ACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCAT                                                                                    <  SRR4302540.2199196/94‑1 (MQ=60)
            cagtgtccgcCAGGCCGATAGCTGCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATC                                                                                   <  SRR4302540.4611548/84‑1 (MQ=60)
             AGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCC                                                                                  >  SRR4302540.3262125/1‑94 (MQ=60)
             AGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCC                                                                                  >  SRR4302540.4183619/1‑94 (MQ=60)
             AGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCC                                                                                  >  SRR4302540.4183620/1‑94 (MQ=60)
             AGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCC                                                                                  >  SRR4302540.6536158/1‑94 (MQ=60)
              GTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCT                                                                                 >  SRR4302540.4135280/1‑94 (MQ=60)
              GTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCT                                                                                 >  SRR4302540.171384/1‑94 (MQ=60)
              GTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCT                                                                                 >  SRR4302540.8187289/1‑94 (MQ=60)
               TATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTG                                                                                <  SRR4302540.9412852/94‑1 (MQ=60)
                ATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGA                                                                               >  SRR4302540.6422821/1‑94 (MQ=60)
                ATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGA                                                                               <  SRR4302540.2106210/94‑1 (MQ=60)
                ATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGA                                                                               >  SRR4302540.292237/1‑94 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              <  SRR4302540.4292421/94‑1 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              <  SRR4302540.9331997/94‑1 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              >  SRR4302540.7422135/1‑94 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              <  SRR4302540.8165882/94‑1 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              <  SRR4302540.4338929/94‑1 (MQ=60)
                 TCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGAT                                                                              <  SRR4302540.1137084/94‑1 (MQ=60)
                  CCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATG                                                                             <  SRR4302540.6526368/94‑1 (MQ=60)
                   CATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGG                                                                            <  SRR4302540.3756932/94‑1 (MQ=60)
                   CATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGG                                                                            >  SRR4302540.8266873/1‑94 (MQ=60)
                   CATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGG                                                                            <  SRR4302540.7253298/94‑1 (MQ=60)
                   CATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATAG                                                                            >  SRR4302540.1287517/1‑94 (MQ=60)
                    ATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGG                                                                           <  SRR4302540.2345640/94‑1 (MQ=60)
                    ATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGG                                                                           >  SRR4302540.3979571/1‑94 (MQ=60)
                    ATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGG                                                                           >  SRR4302540.5765919/1‑94 (MQ=60)
                    ATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGG                                                                           >  SRR4302540.9504950/1‑94 (MQ=60)
                    ATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGG                                                                           <  SRR4302540.1470882/94‑1 (MQ=60)
                     TCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGC                                                                          <  SRR4302540.2398831/94‑1 (MQ=60)
                     TCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCACTTCGCATCCTGATGGGC                                                                          >  SRR4302540.9097522/1‑94 (MQ=60)
                      CAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCA                                                                         >  SRR4302540.7464062/1‑94 (MQ=60)
                       AGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAC                                                                        >  SRR4302540.5378537/1‑94 (MQ=60)
                       AGGCCGATAGGTTCACGGTCTAAATTACGGGTAGTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCAC                                                                        <  SRR4302540.8365686/94‑1 (MQ=60)
                          CCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCC                                                                     >  SRR4302540.8154379/1‑94 (MQ=60)
                           CGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCA                                                                    >  SRR4302540.9610093/1‑94 (MQ=60)
                           CGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCA                                                                    >  SRR4302540.4013622/1‑94 (MQ=60)
                           CGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCA                                                                    >  SRR4302540.4329982/1‑94 (MQ=60)
                           CGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCA                                                                    >  SRR4302540.548734/1‑94 (MQ=60)
                             ATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATT                                                                  <  SRR4302540.5272689/94‑1 (MQ=60)
                              TAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTC                                                                 >  SRR4302540.7569789/1‑94 (MQ=60)
                               AGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCC                                                                <  SRR4302540.153369/94‑1 (MQ=60)
                               AGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCC                                                                >  SRR4302540.2069191/1‑94 (MQ=60)
                                GGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACTCCATTCCT                                                               <  SRR4302540.6678970/94‑1 (MQ=60)
                                 GTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTC                                                              >  SRR4302540.4164756/1‑94 (MQ=60)
                                  TTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCA                                                             >  SRR4302540.9477487/1‑94 (MQ=60)
                                  TTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGGGCCCCATTTCGCATCCTGATGGGCACGCCATTCCTCG                                                             >  SRR4302540.2269456/1‑94 (MQ=60)
                                   TCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAG                                                            <  SRR4302540.9239023/94‑1 (MQ=60)
                                    CACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGT                                                           >  SRR4302540.3151755/1‑94 (MQ=60)
                                       GGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAT                                                        >  SRR4302540.165817/1‑94 (MQ=60)
                                       GGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAT                                                        <  SRR4302540.2795989/94‑1 (MQ=60)
                                       GGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAT                                                        >  SRR4302540.6074342/1‑94 (MQ=60)
                                       CGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAAT                                                        <  SRR4302540.2795990/94‑1 (MQ=60)
                                        GTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATT                                                       <  SRR4302540.3699856/94‑1 (MQ=60)
                                        GTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATT                                                       <  SRR4302540.5135761/94‑1 (MQ=60)
                                         TCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTTTGAGGTGCACCAATTCGCATCCTGATGGTCACGCCATTCCTCAGCTATTA                                                      >  SRR4302540.9077422/1‑94 (MQ=60)
                                          CTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTAC                                                     >  SRR4302540.5885008/1‑94 (MQ=60)
                                           TAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACG                                                    <  SRR4302540.5387984/94‑1 (MQ=60)
                                             AATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCC                                                  <  SRR4302540.8295696/94‑1 (MQ=60)
                                               atatGGGTATTTTCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCAC                                                <  SRR4302540.2714866/90‑1 (MQ=60)
                                                TACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACG                                               <  SRR4302540.2308140/94‑1 (MQ=60)
                                                   GGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTC                                            >  SRR4302540.9533296/1‑94 (MQ=60)
                                                   GGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTC                                            >  SRR4302540.2738540/1‑94 (MQ=60)
                                                   GGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTC                                            >  SRR4302540.1823562/1‑94 (MQ=60)
                                                    GGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCA                                           <  SRR4302540.5859691/94‑1 (MQ=60)
                                                     GTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAG                                          <  SRR4302540.3019690/94‑1 (MQ=60)
                                                     GTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAG                                          >  SRR4302540.3815030/1‑94 (MQ=60)
                                                       ATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCG                                        <  SRR4302540.520546/94‑1 (MQ=60)
                                                       ATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCG                                        >  SRR4302540.7747554/1‑94 (MQ=60)
                                                       ATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCG                                        <  SRR4302540.2623494/94‑1 (MQ=60)
                                                       ATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCG                                        >  SRR4302540.7746332/1‑94 (MQ=60)
                                                        TTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGT                                       >  SRR4302540.2373174/1‑94 (MQ=60)
                                                         TTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTC                                      >  SRR4302540.7161850/1‑94 (MQ=60)
                                                         TTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTC                                      >  SRR4302540.7321065/1‑94 (MQ=60)
                                                         TTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTC                                      >  SRR4302540.1256516/1‑94 (MQ=60)
                                                          TGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCT                                     >  SRR4302540.1982356/1‑94 (MQ=60)
                                                          TGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCT                                     <  SRR4302540.2684299/94‑1 (MQ=60)
                                                          TGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCT                                     >  SRR4302540.4672002/1‑94 (MQ=60)
                                                          TGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCT                                     <  SRR4302540.2047780/94‑1 (MQ=60)
                                                           GCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTA                                    <  SRR4302540.4104936/94‑1 (MQ=60)
                                                           GCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTA                                    >  SRR4302540.28729/1‑94 (MQ=60)
                                                             CAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACT                                  <  SRR4302540.3524582/94‑1 (MQ=60)
                                                             CAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACT                                  >  SRR4302540.374689/1‑94 (MQ=60)
                                                              AGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTT                                 <  SRR4302540.9784681/94‑1 (MQ=60)
                                                              AGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTT                                 <  SRR4302540.9784682/94‑1 (MQ=60)
                                                              AGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTT                                 >  SRR4302540.10064097/1‑94 (MQ=60)
                                                              AGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTT                                 >  SRR4302540.5344574/1‑94 (MQ=60)
                                                               GCAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTC                                <  SRR4302540.10147425/94‑1 (MQ=60)
                                                                CAGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCG                               <  SRR4302540.3082792/94‑1 (MQ=60)
                                                                 AGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGA                              <  SRR4302540.2749361/94‑1 (MQ=60)
                                                                 AGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGA                              >  SRR4302540.8463002/1‑94 (MQ=60)
                                                                 AGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGA                              <  SRR4302540.7009544/94‑1 (MQ=60)
                                                                 AGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGA                              >  SRR4302540.4013377/1‑94 (MQ=60)
                                                                 AGGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGA                              >  SRR4302540.7783695/1‑94 (MQ=60)
                                                                  GGAAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAG                             <  SRR4302540.3544737/94‑1 (MQ=60)
                                                                    AAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTT                           >  SRR4302540.3202866/1‑94 (MQ=60)
                                                                    AAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTT                           >  SRR4302540.4377315/1‑94 (MQ=60)
                                                                    AAAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTT                           <  SRR4302540.1999911/94‑1 (MQ=60)
                                                                     AAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTT                          <  SRR4302540.1596028/94‑1 (MQ=60)
                                                                     AAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTT                          >  SRR4302540.5746197/1‑94 (MQ=60)
                                                                     AAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTT                          >  SRR4302540.6710404/1‑94 (MQ=60)
                                                                     AAGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTT                          <  SRR4302540.3336435/94‑1 (MQ=60)
                                                                      AGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTT                         <  SRR4302540.2463890/94‑1 (MQ=60)
                                                                      AGGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTT                         <  SRR4302540.179465/94‑1 (MQ=60)
                                                                       GGGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTT                        <  SRR4302540.5088580/94‑1 (MQ=60)
                                                                        GGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTC                       <  SRR4302540.8648740/94‑1 (MQ=60)
                                                                        GGCGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTC                       <  SRR4302540.760631/94‑1 (MQ=60)
                                                                          CGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAA                     <  SRR4302540.942921/94‑1 (MQ=60)
                                                                          CGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAA                     <  SRR4302540.1296669/94‑1 (MQ=60)
                                                                          CGGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAA                     <  SRR4302540.9416699/94‑1 (MQ=60)
                                                                           GGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAA                    >  SRR4302540.1679332/1‑94 (MQ=60)
                                                                           GGCGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACACCACGTTCAGCGTCTACTTCGAGTTTTTCAAA                    >  SRR4302540.1191230/1‑94 (MQ=60)
                                                                             CGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGG                  <  SRR4302540.1720829/94‑1 (MQ=60)
                                                                             CGCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGG                  <  SRR4302540.4271314/94‑1 (MQ=60)
                                                                              GCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGT                 <  SRR4302540.3486397/94‑1 (MQ=60)
                                                                              GCGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGT                 >  SRR4302540.9199290/1‑94 (MQ=60)
                                                                               CGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTT                >  SRR4302540.9858062/1‑94 (MQ=60)
                                                                               CGCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTT                <  SRR4302540.2677335/94‑1 (MQ=60)
                                                                                GCCTGAGAGGTTCCCCATTTCGCATTCTGATGTGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               <  SRR4302540.5645249/94‑1 (MQ=60)
                                                                                GCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               >  SRR4302540.3937067/1‑94 (MQ=60)
                                                                                GCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               >  SRR4302540.4523342/1‑94 (MQ=60)
                                                                                GCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               <  SRR4302540.4640588/94‑1 (MQ=60)
                                                                                GCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               <  SRR4302540.9125228/94‑1 (MQ=60)
                                                                                GCCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTG               >  SRR4302540.4835712/1‑94 (MQ=60)
                                                                                 CCTGAGAGGTGCACCATTTCGCATCCTGATGGGTACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGC              <  SRR4302540.9796411/94‑1 (MQ=60)
                                                                                 CCTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGC              >  SRR4302540.4756341/1‑94 (MQ=60)
                                                                                  CTGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCG             >  SRR4302540.1893090/1‑94 (MQ=60)
                                                                                   TGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGG            >  SRR4302540.5018218/1‑94 (MQ=60)
                                                                                   TGAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGG            >  SRR4302540.9497436/1‑94 (MQ=60)
                                                                                    GAGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGC           <  SRR4302540.4074807/94‑1 (MQ=60)
                                                                                     AGAGGTTCGCCATTTCGCATCCGGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCA          <  SRR4302540.8040321/94‑1 (MQ=60)
                                                                                     AGAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCA          <  SRR4302540.9536810/94‑1 (MQ=60)
                                                                                      GAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAT         >  SRR4302540.4209682/1‑94 (MQ=60)
                                                                                      GAGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCAACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCAT         <  SRR4302540.4170402/94‑1 (MQ=60)
                                                                                       AGGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATT        >  SRR4302540.2564530/1‑94 (MQ=60)
                                                                                        GGTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTA       <  SRR4302540.7453508/94‑1 (MQ=60)
                                                                                         GTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC      >  SRR4302540.4943461/1‑94 (MQ=60)
                                                                                         GTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC      <  SRR4302540.796734/94‑1 (MQ=60)
                                                                                         GTGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTAC      <  SRR4302540.5050231/94‑1 (MQ=60)
                                                                                          TGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACT     <  SRR4302540.4439232/94‑1 (MQ=60)
                                                                                          TGCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACT     >  SRR4302540.9073427/1‑94 (MQ=60)
                                                                                           GCACCATTTCGCATCCTGGTGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCTTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTC    <  SRR4302540.7034487/94‑1 (MQ=60)
                                                                                           GCACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTC    <  SRR4302540.7034483/94‑1 (MQ=60)
                                                                                            CACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCC   <  SRR4302540.6119559/94‑1 (MQ=60)
                                                                                             ACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCCA  <  SRR4302540.4180125/94‑1 (MQ=60)
                                                                                             ACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCCA  >  SRR4302540.4621745/1‑94 (MQ=60)
                                                                                             ACCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCCA  >  SRR4302540.8928925/1‑94 (MQ=60)
                                                                                             |                                                                                             
TGACCGTGGCGACAGTATCCATCAGGCCGATAGGTTCACGGTCTAAATTACGGGTATTTGCCAGCAGGAAAGGGCGGCGCGCCTGAGAGGTGCGCCATTTCGCATCCTGATGGGCACGCCATTCCTCAGTAATTACGCCACGTTCAGCGTCTACTTCGAGTTTTTCAAAGGTTGCGGCATTACTCCA  >  NC_000913/1572605‑1572791

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