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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 294,929 C→A 87.1% E573* (GAA→TAA)  paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273294,9290CA87.1% 349.7 / 48.7 147E573* (GAA→TAA) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base C (7/12);  new base A (57/71);  total (64/83)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 150 total aligned reads displayed.

GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTT  >  CP009273/294784‑295073
                                                                                                                                                 |                                                                                                                                                
gTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTaaa                                                                                                                                                <  1:1842495/148‑1 (MQ=255)
gTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTaaa                                                                                                                                                <  2:3282959/148‑1 (MQ=255)
  tttCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATg                                                                                                                                              <  1:3039589/148‑1 (MQ=255)
    tCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGct                                                                                                                                            <  2:591421/148‑1 (MQ=255)
      ggATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGctct                                                                                                                                          <  1:1169011/148‑1 (MQ=255)
          tCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTc                                                                                                                                      <  2:1177139/148‑1 (MQ=255)
             ggATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGc                                                                                                                                   <  1:2495301/148‑1 (MQ=255)
              gATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCt                                                                                                                                  <  1:840344/148‑1 (MQ=255)
              gATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCt                                                                                                                                  <  1:2523770/148‑1 (MQ=255)
              gATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCt                                                                                                                                  >  1:1158994/1‑148 (MQ=255)
              gATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCt                                                                                                                                  >  2:1363031/1‑148 (MQ=255)
               aTGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTg                                                                                                                                 <  2:2935166/148‑1 (MQ=255)
               aTGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTg                                                                                                                                 >  2:1787917/1‑148 (MQ=255)
                  cGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCa                                                                                                                              >  1:3210770/1‑148 (MQ=255)
                  cGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCa                                                                                                                              <  1:559245/148‑1 (MQ=255)
                  cGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCa                                                                                                                              >  1:193238/1‑148 (MQ=255)
                    tCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGt                                                                                                                            >  1:733156/1‑148 (MQ=255)
                         cagcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgcc                                                                                                                       >  2:2525104/1‑148 (MQ=255)
                         cagcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgcc                                                                                                                       >  2:1657643/1‑148 (MQ=255)
                          agcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgccg                                                                                                                      >  1:296083/1‑148 (MQ=255)
                          agcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgccg                                                                                                                      >  1:296092/1‑148 (MQ=255)
                           gcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgccgc                                                                                                                     <  2:1882702/148‑1 (MQ=255)
                           gcaCACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgccgc                                                                                                                     <  2:310520/148‑1 (MQ=255)
                            cacacacaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTgccgcc                                                                                                                    <  1:1354148/148‑1 (MQ=255)
                               acacaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCAACAAAATGCGCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGGCTGc                                                                                                                 <  2:17815/148‑1 (MQ=255)
                                 acaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCgg                                                                                                               <  1:1817244/148‑1 (MQ=255)
                                 acaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCgg                                                                                                               >  2:2010494/1‑148 (MQ=255)
                                 acaGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCgg                                                                                                               >  2:2011396/1‑148 (MQ=255)
                                        gCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTc                                                                                                        <  2:831705/148‑1 (MQ=255)
                                           tACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATg                                                                                                     <  2:2059039/148‑1 (MQ=255)
                                           tACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATg                                                                                                     <  2:2517635/148‑1 (MQ=255)
                                            aCGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGt                                                                                                    >  1:518925/1‑148 (MQ=255)
                                             cGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTa                                                                                                   <  2:3218415/148‑1 (MQ=255)
                                                  gACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTg                                                                                              >  2:1100800/1‑148 (MQ=255)
                                                   aCCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTgg                                                                                             <  2:2614048/148‑1 (MQ=255)
                                                   aCCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTgg                                                                                             >  2:1648410/1‑148 (MQ=255)
                                                      cGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGgcgc                                                                                          <  1:27916/148‑1 (MQ=255)
                                                      cGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGgcgc                                                                                          >  2:1636332/1‑148 (MQ=255)
                                                        aaCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAACTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        <  2:850484/148‑1 (MQ=255)
                                                        aaCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        >  1:560087/1‑148 (MQ=255)
                                                        aaCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        >  1:560083/1‑148 (MQ=255)
                                                        aaCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        <  2:2700928/148‑1 (MQ=255)
                                                        aaCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        <  2:2700934/148‑1 (MQ=255)
                                                        aaCGTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCtt                                                                                        <  2:3075774/148‑1 (MQ=255)
                                                         aCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTc                                                                                       >  1:3245877/1‑148 (MQ=255)
                                                         aCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTc                                                                                       >  2:747936/1‑148 (MQ=255)
                                                         aCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTc                                                                                       >  1:2079617/1‑148 (MQ=255)
                                                         aCTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCGTc                                                                                       >  2:2316723/1‑148 (MQ=255)
                                                           ttCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCgg                                                                                     <  2:2015229/148‑1 (MQ=255)
                                                           gtctcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCACAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCgg                                                                                     <  2:1254813/147‑1 (MQ=255)
                                                            tctcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCggg                                                                                    >  1:2265882/1‑148 (MQ=255)
                                                             ctcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGt                                                                                   <  1:3052138/148‑1 (MQ=255)
                                                             ctcCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGt                                                                                   >  1:734099/1‑148 (MQ=255)
                                                                 cGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTCCCg                                                                               >  2:2098912/1‑148 (MQ=255)
                                                                    cgcgCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTg                                                                            <  2:256854/148‑1 (MQ=255)
                                                                    cgcgCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTg                                                                            <  1:955158/148‑1 (MQ=255)
                                                                      cgcTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGt                                                                          >  1:2890256/1‑148 (MQ=255)
                                                                       gcTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTa                                                                         <  2:312500/148‑1 (MQ=255)
                                                                         tATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAAt                                                                       <  1:3184230/148‑1 (MQ=255)
                                                                             cACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTg                                                                   <  2:193238/148‑1 (MQ=255)
                                                                              aCGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGc                                                                  <  1:1348501/148‑1 (MQ=255)
                                                                               cGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGcc                                                                 >  2:609234/1‑148 (MQ=255)
                                                                               cGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGcc                                                                 >  1:1753729/1‑148 (MQ=255)
                                                                                       tCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCaa                                                         <  2:1093496/148‑1 (MQ=255)
                                                                                        ccACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTTCCGTCGGCaaa                                                        >  1:289122/1‑148 (MQ=255)
                                                                                        ccACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAAGGCTCTCTTCCGCTGTCAGTCTGCCGCCTTCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCaaa                                                        <  2:689971/148‑1 (MQ=255)
                                                                                         cacaAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAAc                                                       <  2:89985/148‑1 (MQ=255)
                                                                                         cacaAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAAc                                                       <  2:90065/148‑1 (MQ=255)
                                                                                             aaaTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCg                                                   <  1:1514671/148‑1 (MQ=255)
                                                                                              aaTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGa                                                  <  2:2227995/148‑1 (MQ=255)
                                                                                               aTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGAc                                                 >  1:3133411/1‑148 (MQ=255)
                                                                                               aTGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGAc                                                 >  1:2745199/1‑148 (MQ=255)
                                                                                                tGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACt                                                >  2:1134801/1‑148 (MQ=255)
                                                                                                 gCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTg                                               <  2:2913805/148‑1 (MQ=255)
                                                                                                 gCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTg                                               >  1:819863/1‑148 (MQ=255)
                                                                                                  cccGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGc                                              <  2:1541062/148‑1 (MQ=255)
                                                                                                  cccGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGc                                              <  2:1541112/148‑1 (MQ=255)
                                                                                                       cAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAgg                                         <  1:389217/148‑1 (MQ=255)
                                                                                                       cAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAgg                                         <  1:2159897/148‑1 (MQ=255)
                                                                                                         aaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAt                                       >  2:3290657/1‑148 (MQ=255)
                                                                                                         aaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAt                                       <  2:2240922/148‑1 (MQ=255)
                                                                                                         aaGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAt                                       <  2:1753729/148‑1 (MQ=255)
                                                                                                            gTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCaa                                    >  1:191699/1‑148 (MQ=255)
                                                                                                                  tgctgGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGa                              >  2:1652769/1‑148 (MQ=255)
                                                                                                                        tACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCg                        <  2:2990032/148‑1 (MQ=255)
                                                                                                                        tACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCg                        >  1:2423759/1‑148 (MQ=255)
                                                                                                                         aCTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGa                       >  2:521208/1‑148 (MQ=255)
                                                                                                                           tctTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAgg                     >  1:2197729/1‑148 (MQ=255)
                                                                                                                           tctTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAgg                     >  1:2197752/1‑148 (MQ=255)
                                                                                                                            ctTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGc                    <  2:24247/148‑1 (MQ=255)
                                                                                                                               tGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAAt                 <  2:296083/148‑1 (MQ=255)
                                                                                                                               tGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAAt                 <  1:1636332/148‑1 (MQ=255)
                                                                                                                                gCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATc                <  2:518925/148‑1 (MQ=255)
                                                                                                                                 cTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTATTCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCa               <  1:451256/148‑1 (MQ=255)
                                                                                                                                     gtgtTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTc           >  1:1474342/1‑148 (MQ=255)
                                                                                                                                         tCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGa       <  1:2112864/148‑1 (MQ=255)
                                                                                                                                         tCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGa       <  2:1638411/148‑1 (MQ=255)
                                                                                                                                         tCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGa       <  1:2112933/148‑1 (MQ=255)
                                                                                                                                         tCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGa       >  2:1176135/1‑148 (MQ=255)
                                                                                                                                              aTTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCtt  >  2:2099181/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTT  >  CP009273/294784‑295073

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAC  >  CP009273/294784‑295078
                                                                                                                                                 |                                                                                                                                                     
GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATG                                                                                                                                                   <  SRR19973070.1243805/150‑1 (MQ=60)
GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATG                                                                                                                                                   <  SRR19973070.1855139/150‑1 (MQ=60)
  TTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCT                                                                                                                                                 <  SRR19973070.3058914/150‑1 (MQ=60)
     CGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTC                                                                                                                                              >  SRR19973070.2932633/1‑150 (MQ=60)
      GGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCT                                                                                                                                             <  SRR19973070.1177197/150‑1 (MQ=60)
            AGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCT                                                                                                                                       >  SRR19973070.1167098/1‑150 (MQ=60)
             GGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTG                                                                                                                                      <  SRR19973070.2511483/150‑1 (MQ=60)
             GGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTG                                                                                                                                      <  SRR19973070.2511473/150‑1 (MQ=60)
              GATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGT                                                                                                                                     <  SRR19973070.2540235/150‑1 (MQ=60)
              GATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGT                                                                                                                                     <  SRR19973070.846030/150‑1 (MQ=60)
               ATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTC                                                                                                                                    <  SRR19973070.590824/150‑1 (MQ=60)
                TGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCA                                                                                                                                   >  SRR19973070.3231536/1‑150 (MQ=60)
                TGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCA                                                                                                                                   >  SRR19973070.194727/1‑150 (MQ=60)
                  CGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGT                                                                                                                                 >  SRR19973070.738347/1‑150 (MQ=60)
                  CGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGT                                                                                                                                 <  SRR19973070.562923/150‑1 (MQ=60)
                     CCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCNG                                                                                                                              <  SRR19973070.2066627/150‑1 (MQ=60)
                        GCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCG                                                                                                                           >  SRR19973070.298129/1‑150 (MQ=60)
                        GCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCG                                                                                                                           >  SRR19973070.298138/1‑150 (MQ=60)
                           GCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCT                                                                                                                        <  SRR19973070.2193552/150‑1 (MQ=60)
                            CACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTG                                                                                                                       <  SRR19973070.1363409/150‑1 (MQ=60)
                             ACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGC                                                                                                                      >  SRR19973070.94823/1‑150 (MQ=60)
                                 ACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTG                                                                                                                  <  SRR19973070.1829767/150‑1 (MQ=60)
                                   AGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGC                                                                                                                <  SRR19973070.2554171/150‑1 (MQ=60)
                                          ATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGT                                                                                                         >  SRR19973070.522465/1‑150 (MQ=60)
                                           TACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTA                                                                                                        >  SRR19973070.1445342/1‑150 (MQ=60)
                                                      CGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTT                                                                                             >  SRR19973070.563779/1‑150 (MQ=60)
                                                      CGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTT                                                                                             >  SRR19973070.563774/1‑150 (MQ=60)
                                                      CGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTT                                                                                             <  SRR19973070.28294/150‑1 (MQ=60)
                                                      CGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCGT                                                                                             >  SRR19973070.563767/1‑150 (MQ=60)
                                                       GAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTC                                                                                            <  SRR19973070.2990079/150‑1 (MQ=60)
                                                       GAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTC                                                                                            >  SRR19973070.3266795/1‑150 (MQ=60)
                                                       GAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTC                                                                                            >  SRR19973070.2093611/1‑150 (MQ=60)
                                                          CTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGG                                                                                         >  SRR19973070.2280922/1‑150 (MQ=60)
                                                           TTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGT                                                                                        >  SRR19973070.739295/1‑150 (MQ=60)
                                                            TCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTAGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTA                                                                                       <  SRR19973070.2827589/150‑1 (MQ=60)
                                                             CTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTAC                                                                                      <  SRR19973070.3071591/150‑1 (MQ=60)
                                                                    CGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGT                                                                               <  SRR19973070.961976/150‑1 (MQ=60)
                                                                    CGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGT                                                                               >  SRR19973070.2908962/1‑150 (MQ=60)
                                                                        CTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATC                                                                           <  SRR19973070.179598/150‑1 (MQ=60)
                                                                         TATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCT                                                                          >  SRR19973070.313681/1‑150 (MQ=60)
                                                                         TATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCT                                                                          <  SRR19973070.3204831/150‑1 (MQ=60)
                                                                             CACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGGTGGTAATCTTGCC                                                                      >  SRR19973070.1765880/1‑150 (MQ=60)
                                                                              ACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCG                                                                     <  SRR19973070.1357747/150‑1 (MQ=60)
                                                                                      CTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTTCCGTCGGCAAA                                                             >  SRR19973070.291077/1‑150 (MQ=60)
                                                                                             AAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGAC                                                      <  SRR19973070.1525278/150‑1 (MQ=60)
                                                                                             AAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGAC                                                      >  SRR19973070.2762849/1‑150 (MQ=60)
                                                                                             AAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGAC                                                      >  SRR19973070.3153650/1‑150 (MQ=60)
                                                                                              AATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACT                                                     <  SRR19973070.2233593/150‑1 (MQ=60)
                                                                                               ATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTG                                                    >  SRR19973070.825458/1‑150 (MQ=60)
                                                                                                       CAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAT                                            <  SRR19973070.392068/150‑1 (MQ=60)
                                                                                                       CAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGAT                                            <  SRR19973070.2174254/150‑1 (MQ=60)
                                                                                                          AGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAA                                         >  SRR19973070.193174/1‑150 (MQ=60)
                                                                                                                      GGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCG                             >  SRR19973070.770150/1‑150 (MQ=60)
                                                                                                                      GGTACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCG                             >  SRR19973070.2439693/1‑150 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          <  SRR19973070.3224218/150‑1 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          >  SRR19973070.3282466/1‑150 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          >  SRR19973070.2212226/1‑150 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          >  SRR19973070.1920111/1‑150 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          >  SRR19973070.1635449/1‑150 (MQ=60)
                                                                                                                         ACTCTTTGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGG                          >  SRR19973070.2212249/1‑150 (MQ=60)
                                                                                                                               TGCTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCA                    <  SRR19973070.1647672/150‑1 (MQ=60)
                                                                                                                                 CTCAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTATTCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATT                  <  SRR19973070.454508/150‑1 (MQ=60)
                                                                                                                                   CAGTGTTCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTC                >  SRR19973070.1484655/1‑150 (MQ=60)
                                                                                                                                         TCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAG          <  SRR19973070.2127030/150‑1 (MQ=60)
                                                                                                                                         TCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAG          <  SRR19973070.2127069/150‑1 (MQ=60)
                                                                                                                                         TCCGAATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAG          <  SRR19973070.2127100/150‑1 (MQ=60)
                                                                                                                                           CGAATTATATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCT        <  SRR19973070.522097/150‑1 (MQ=60)
                                                                                                                                              ATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCA     >  SRR19973070.687620/1‑150 (MQ=60)
                                                                                                                                              ATTAAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCA     >  SRR19973070.1846031/1‑150 (MQ=60)
                                                                                                                                                 AAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAC  >  SRR19973070.1014881/1‑150 (MQ=60)
                                                                                                                                                 |                                                                                                                                                     
GTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCGAACTTCTCCCGTCGCGCTATGCACGCCGACCTCCACAAAATGCCCGGCAAAGGTCGACTGCTGGTACTCTTTGCTCAGTGTTCCGAATTCAATGCTCTCTTCCGCTGTCAGTCTGCCGCCTGCGGTGGCTTCATGTAGCGTGGCGCTTCGGGTACCGTTGGTAATCTTGCCGTCGGCAAACTGCGACTGCTCAGGATCAAACCCGACTGCCGAGGCAATCATTTCGCGAAGCTTCATAC  >  CP009273/294784‑295078

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