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 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,822 2 bp→CG coding (578‑579/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,8220GC100.0% 793.5 / NA 248G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (136/112);  total (136/112)
*NC_0009133,987,8230CG100.0% 788.7 / NA 248G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (136/112);  total (136/112)

Only 100 of 256 total aligned reads displayed.

CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909
                                                                                       ||                                                                                      
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                                  >  2:644175/1‑93 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                                  >  2:1656165/1‑93 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                                  >  2:1995736/1‑93 (MQ=255)
    tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                                >  2:4852752/1‑93 (MQ=255)
     cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                               >  2:4306247/1‑93 (MQ=255)
     cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                               <  1:6049407/93‑1 (MQ=255)
      gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAt                                                                              <  1:787089/93‑1 (MQ=255)
      gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAt                                                                              <  1:4924000/93‑1 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            >  2:563186/1‑93 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            <  1:3999724/93‑1 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            <  1:1916178/93‑1 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAAAGCGGCGCGAATAc                                                                            >  2:4935404/1‑93 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTATCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           >  1:5965540/1‑93 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           <  2:690104/93‑1 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           <  1:4345101/93‑1 (MQ=255)
          gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                                          >  2:5391188/1‑93 (MQ=255)
          gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                                          >  2:4466365/1‑93 (MQ=255)
             cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                       <  2:2598591/93‑1 (MQ=255)
             cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                       >  2:1630495/1‑93 (MQ=255)
             cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                       >  2:1179805/1‑93 (MQ=255)
               tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTc                                                                     <  1:5143393/93‑1 (MQ=255)
               tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTc                                                                     >  2:5415010/1‑93 (MQ=255)
                gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                    >  1:4040457/1‑93 (MQ=255)
                 aaTCATCAGGGCGTCATTCAATACCCACCGCACCTTGTCCTGCCGCCGAAAGAGAGATCTGGGGTGGCAACGCGGCGCTAATACGTGACTCCa                                                                   <  2:1450659/93‑1 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   >  1:5033107/1‑93 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   <  2:933107/93‑1 (MQ=255)
                  atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                                  >  1:6165459/1‑93 (MQ=255)
                  atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                                  >  1:4783697/1‑93 (MQ=255)
                   tcatcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                                 >  2:5444299/1‑93 (MQ=255)
                     atcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAcc                                                               >  1:4765657/1‑93 (MQ=255)
                      tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                              >  1:3846189/1‑93 (MQ=255)
                      tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                              >  1:3419784/1‑93 (MQ=255)
                      tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                              >  2:5835926/1‑93 (MQ=255)
                       caAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTa                                                             <  1:1327092/93‑1 (MQ=255)
                        aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                            <  1:4594173/93‑1 (MQ=255)
                        aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                            <  1:5630042/93‑1 (MQ=255)
                         aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                           >  2:829827/1‑93 (MQ=255)
                           gcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                         >  2:5525436/1‑93 (MQ=255)
                            cggcATTCAATACCCACCGCCTCTTGTCCTACCGCCGGAGGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                                        <  2:2360984/93‑1 (MQ=255)
                             ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                       >  2:255008/1‑93 (MQ=255)
                              gcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGttt                                                      >  2:6267732/1‑93 (MQ=255)
                              gcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGttt                                                      <  1:1543623/93‑1 (MQ=255)
                               cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                     <  2:2905454/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  2:3007124/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  1:5344881/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  2:3340772/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    >  1:233232/1‑93 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  2:2167325/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  1:3322776/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  1:866275/93‑1 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   >  1:1022167/1‑93 (MQ=255)
                                  tCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                  >  2:3535723/1‑93 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:1351006/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 >  1:4635813/1‑93 (MQ=255)
                                    aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                                >  2:5989731/1‑93 (MQ=255)
                                    aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                                <  2:3790570/93‑1 (MQ=255)
                                      tACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                              >  1:197812/1‑93 (MQ=255)
                                       aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                             <  1:71583/93‑1 (MQ=255)
                                        cccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                            >  2:3281580/1‑93 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           <  1:1155048/93‑1 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           <  1:6008970/93‑1 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           >  1:4484745/1‑93 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           >  2:5166744/1‑93 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           >  1:2551699/1‑93 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           <  1:4123132/93‑1 (MQ=255)
                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                                         >  2:183507/1‑93 (MQ=255)
                                             cGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                       <  2:1734547/93‑1 (MQ=255)
                                             cGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                       >  2:5666008/1‑93 (MQ=255)
                                              gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                      <  2:2103784/93‑1 (MQ=255)
                                              gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                      >  1:944759/1‑93 (MQ=255)
                                               cACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                                     <  2:3091108/93‑1 (MQ=255)
                                                aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                                    >  1:4195366/1‑93 (MQ=255)
                                                aCCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                                    >  1:1643520/1‑93 (MQ=255)
                                                 ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAACCGTTTTAGTCCGGCTGCGGCAAGa                                   >  1:2918052/1‑93 (MQ=255)
                                                 ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                   >  1:2375626/1‑93 (MQ=255)
                                                  cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                  >  2:388475/1‑93 (MQ=255)
                                                   ttGTCCTACCGCGGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTAAGTCCGGCTACGGCAAGAat                                 <  1:3001061/93‑1 (MQ=255)
                                                   ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                 >  1:10076/1‑93 (MQ=255)
                                                   ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                 <  1:3756248/93‑1 (MQ=255)
                                                   ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                 >  2:4241126/1‑93 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                <  1:3877389/93‑1 (MQ=255)
                                                     gTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                               <  1:6247346/93‑1 (MQ=255)
                                                         tACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           <  1:2187817/93‑1 (MQ=255)
                                                            cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        <  1:1704210/93‑1 (MQ=255)
                                                            cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        <  2:4195733/93‑1 (MQ=255)
                                                             gccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg                       >  1:5767039/1‑93 (MQ=255)
                                                              ccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                      <  1:3256025/93‑1 (MQ=255)
                                                               cgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                     >  2:2409060/1‑93 (MQ=255)
                                                                ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    <  1:4239295/93‑1 (MQ=255)
                                                                 gAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   <  2:4395944/93‑1 (MQ=255)
                                                                 gAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   >  2:5032332/1‑93 (MQ=255)
                                                                   agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                 >  2:3647958/1‑93 (MQ=255)
                                                                   agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                 >  2:4503793/1‑93 (MQ=255)
                                                                    gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTctc                <  1:2688392/93‑3 (MQ=255)
                                                                        gaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGtt            >  1:4313398/2‑93 (MQ=255)
                                                                          tctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAt          >  2:4111941/1‑93 (MQ=255)
                                                                          tctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAt          >  1:2747925/1‑93 (MQ=255)
                                                                             cGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa       <  2:4875997/93‑1 (MQ=255)
                                                                              gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg      <  2:5945761/93‑1 (MQ=255)
                                                                              gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg      <  1:3531363/93‑1 (MQ=255)
                                                                                       ||                                                                                      
CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909

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
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGC  >  NC_000913/3987729‑3987910
                                                                                             |                                                                                        
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTTCTACCGCCGGAAGAGAGATCTCGGGTGGCAAC                                                                                          >  SRR4302533.5990964/1‑94 (MQ=60)
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAAC                                                                                          >  SRR4302533.206949/1‑94 (MQ=60)
     GCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGC                                                                                     <  SRR4302533.3854472/94‑1 (MQ=60)
     GCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGC                                                                                     <  SRR4302533.5166238/94‑1 (MQ=60)
     GCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGC                                                                                     <  SRR4302533.847628/94‑1 (MQ=60)
      CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCG                                                                                    <  SRR4302533.4038193/94‑1 (MQ=60)
      CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCG                                                                                    <  SRR4302533.5570291/94‑1 (MQ=60)
       AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGC                                                                                   >  SRR4302533.5142570/1‑94 (MQ=60)
       AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGC                                                                                   <  SRR4302533.1072320/94‑1 (MQ=60)
        GCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCG                                                                                  <  SRR4302533.878636/94‑1 (MQ=60)
          TCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAA                                                                                <  SRR4302533.3785525/94‑1 (MQ=60)
           CGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAT                                                                               <  SRR4302533.6086057/94‑1 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCCCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATA                                                                              <  SRR4302533.1281521/94‑1 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGCGATCTCGGGTGGCAACGCGGCGCGAATA                                                                              >  SRR4302533.2079215/1‑94 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCCACGCGGCGCGAATA                                                                              >  SRR4302533.4843925/1‑94 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATA                                                                              <  SRR4302533.4953829/94‑1 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATA                                                                              <  SRR4302533.791713/94‑1 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTATCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACG                                                                            >  SRR4302533.6001639/1‑94 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACG                                                                            <  SRR4302533.4024076/94‑1 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACG                                                                            <  SRR4302533.1927876/94‑1 (MQ=60)
               AGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGT                                                                           <  SRR4302533.4371287/94‑1 (MQ=60)
                 TGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGA                                                                         >  SRR4302533.6017107/1‑94 (MQ=60)
                 TGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGA                                                                         >  SRR4302533.2072991/1‑94 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                       <  SRR4302533.270275/94‑1 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                       >  SRR4302533.3547225/1‑94 (MQ=60)
                    GTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTC                                                                      >  SRR4302533.244658/1‑94 (MQ=60)
                    GTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTC                                                                      >  SRR4302533.5868333/1‑94 (MQ=60)
                    GTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTC                                                                      >  SRR4302533.4310353/1‑94 (MQ=60)
                     TGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCC                                                                     >  SRR4302533.4502877/1‑94 (MQ=60)
                     TGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCC                                                                     <  SRR4302533.5174477/94‑1 (MQ=60)
                     TGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCC                                                                     >  SRR4302533.4065049/1‑94 (MQ=60)
                      GAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTATCGCCGGAAGAGAGATCTCGGATGGCAACGCGGCGCGAATACGTGACTCCA                                                                    >  SRR4302533.2821818/1‑94 (MQ=60)
                      GAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCA                                                                    >  SRR4302533.5063596/1‑94 (MQ=60)
                       AATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAG                                                                   >  SRR4302533.4812674/1‑94 (MQ=60)
                       AATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAG                                                                   <  SRR4302533.4633639/94‑1 (MQ=60)
                       AATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAG                                                                   >  SRR4302533.6202520/1‑94 (MQ=60)
                        ATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGA                                                                  <  SRR4302533.676190/94‑1 (MQ=60)
                        ATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGA                                                                  >  SRR4302533.1592102/1‑94 (MQ=60)
                         TCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAC                                                                 <  SRR4302533.5420954/94‑1 (MQ=60)
                         TCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAC                                                                 >  SRR4302533.3290081/1‑94 (MQ=60)
                          CATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACC                                                                >  SRR4302533.4794555/1‑94 (MQ=60)
                           ATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCT                                                               >  SRR4302533.3440914/1‑94 (MQ=60)
                           ATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCT                                                               >  SRR4302533.3869676/1‑94 (MQ=60)
                            TCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTA                                                              >  SRR4302533.1938634/1‑94 (MQ=60)
                             CAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAA                                                             <  SRR4302533.1335057/94‑1 (MQ=60)
                             CAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAA                                                             >  SRR4302533.3885924/1‑94 (MQ=60)
                              AAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAA                                                            <  SRR4302533.5664088/94‑1 (MQ=60)
                              AAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAA                                                            <  SRR4302533.4622001/94‑1 (MQ=60)
                              AAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAA                                                            <  SRR4302533.6191831/94‑1 (MQ=60)
                               AGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAC                                                           <  SRR4302533.3227294/94‑1 (MQ=60)
                               AGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAC                                                           <  SRR4302533.5621156/94‑1 (MQ=60)
                               AGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAC                                                           <  SRR4302533.453529/94‑1 (MQ=60)
                                GGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACG                                                          >  SRR4302533.2251326/1‑94 (MQ=60)
                                GGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACG                                                          <  SRR4302533.6259552/94‑1 (MQ=60)
                                   GGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTT                                                       <  SRR4302533.2589760/94‑1 (MQ=60)
                                    GCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTT                                                      <  SRR4302533.1552921/94‑1 (MQ=60)
                                     CATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTA                                                     >  SRR4302533.234673/1‑94 (MQ=60)
                                     CATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTA                                                     <  SRR4302533.6068775/94‑1 (MQ=60)
                                     CATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTA                                                     <  SRR4302533.4634953/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                    >  SRR4302533.2760493/1‑94 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                    <  SRR4302533.3343295/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                    <  SRR4302533.5376935/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                    <  SRR4302533.871350/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                    >  SRR4302533.1028174/1‑94 (MQ=60)
                                        TCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTC                                                  >  SRR4302533.4663901/1‑94 (MQ=60)
                                        TCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTC                                                  <  SRR4302533.4256195/94‑1 (MQ=60)
                                         CAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCC                                                 <  SRR4302533.2366016/94‑1 (MQ=60)
                                         CAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCC                                                 <  SRR4302533.1711035/94‑1 (MQ=60)
                                         CAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCC                                                 >  SRR4302533.1124467/1‑94 (MQ=60)
                                          AATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCG                                                <  SRR4302533.1879684/94‑1 (MQ=60)
                                          AATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCG                                                >  SRR4302533.4339522/1‑94 (MQ=60)
                                           ATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGG                                               >  SRR4302533.1019655/1‑94 (MQ=60)
                                           ATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGG                                               >  SRR4302533.199058/1‑94 (MQ=60)
                                            TACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGC                                              >  SRR4302533.4036911/1‑94 (MQ=60)
                                             ACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCT                                             <  SRR4302533.3098609/94‑1 (MQ=60)
                                             ACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCT                                             >  SRR4302533.4031790/1‑94 (MQ=60)
                                             ACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCT                                             <  SRR4302533.72043/94‑1 (MQ=60)
                                              CCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTA                                            >  SRR4302533.4123965/1‑94 (MQ=60)
                                              CCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTA                                            >  SRR4302533.4511901/1‑94 (MQ=60)
                                              CCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTA                                            >  SRR4302533.2567402/1‑94 (MQ=60)
                                               CCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAC                                           <  SRR4302533.4148166/94‑1 (MQ=60)
                                               CCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAC                                           <  SRR4302533.1161977/94‑1 (MQ=60)
                                               CCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAC                                           <  SRR4302533.6045375/94‑1 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                          >  SRR4302533.755619/1‑94 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                          >  SRR4302533.3534998/1‑94 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                          >  SRR4302533.2847855/1‑94 (MQ=60)
                                                 CCCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGG                                         <  SRR4302533.4442121/94‑1 (MQ=60)
                                                 ACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGG                                         <  SRR4302533.6017009/94‑1 (MQ=60)
                                                  CCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGC                                        >  SRR4302533.4768013/1‑94 (MQ=60)
                                                   CGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA                                       >  SRR4302533.950316/1‑94 (MQ=60)
                                                   CGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA                                       >  SRR4302533.6148564/1‑94 (MQ=60)
                                                     CACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                     >  SRR4302533.194564/1‑94 (MQ=60)
                                                     CACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                     >  SRR4302533.3945386/1‑94 (MQ=60)
                                                     CACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                     >  SRR4302533.1653470/1‑94 (MQ=60)
                                                     CACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                     >  SRR4302533.4220760/1‑94 (MQ=60)
                                                      ACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAACCGTTTTAGTCCGGCTGCGGCAAGA                                    >  SRR4302533.2935937/1‑94 (MQ=60)
                                                      ACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGA                                    <  SRR4302533.171045/94‑1 (MQ=60)
                                                      ACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGA                                    >  SRR4302533.2390326/1‑94 (MQ=60)
                                                       CCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAA                                   >  SRR4302533.914226/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                  >  SRR4302533.3281536/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                  >  SRR4302533.10159/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                  >  SRR4302533.827685/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                  >  SRR4302533.693036/1‑94 (MQ=60)
                                                         TTGTCCTACCGCGGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTAAGTCCGGCTACGGCAAGAATG                                 <  SRR4302533.3019413/94‑1 (MQ=60)
                                                         TTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATG                                 <  SRR4302533.4718870/94‑1 (MQ=60)
                                                         TTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATG                                 <  SRR4302533.3779273/94‑1 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                >  SRR4302533.401357/1‑94 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                <  SRR4302533.237306/94‑1 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                <  SRR4302533.3901052/94‑1 (MQ=60)
                                                           GTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGAT                               >  SRR4302533.4211230/1‑94 (MQ=60)
                                                           GTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGAT                               <  SRR4302533.4804154/94‑1 (MQ=60)
                                                           GTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGAT                               <  SRR4302533.6284732/94‑1 (MQ=60)
                                                              CTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGC                            <  SRR4302533.575079/94‑1 (MQ=60)
                                                              CTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGC                            >  SRR4302533.343449/1‑94 (MQ=60)
                                                              CTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGC                            >  SRR4302533.1143631/1‑94 (MQ=60)
                                                               TACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCA                           <  SRR4302533.2201295/94‑1 (MQ=60)
                                                                 CCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                         >  SRR4302533.5153837/1‑94 (MQ=60)
                                                                  CGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                        <  SRR4302533.4976469/94‑1 (MQ=60)
                                                                  CGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                        <  SRR4302533.1714534/94‑1 (MQ=60)
                                                                  CGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                        >  SRR4302533.5801855/1‑94 (MQ=60)
                                                                    CCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTA                      <  SRR4302533.3276104/94‑1 (MQ=60)
                                                                     CGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAT                     >  SRR4302533.386937/1‑94 (MQ=60)
                                                                     CGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAT                     <  SRR4302533.3720652/94‑1 (MQ=60)
                                                                     CGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAT                     <  SRR4302533.836309/94‑1 (MQ=60)
                                                                      GGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATT                    >  SRR4302533.5946565/1‑94 (MQ=60)
                                                                      GGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATT                    <  SRR4302533.4264898/94‑1 (MQ=60)
                                                                        AAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCG                  >  SRR4302533.1382030/1‑94 (MQ=60)
                                                                          GAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCTCC                <  SRR4302533.2704873/94‑1 (MQ=60)
                                                                          GAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC                <  SRR4302533.4332157/94‑1 (MQ=60)
                                                                          GAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC                >  SRR4302533.518460/1‑94 (MQ=60)
                                                                             AGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTT             >  SRR4302533.4339366/1‑94 (MQ=60)
                                                                               ATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAT           >  SRR4302533.2974467/1‑94 (MQ=60)
                                                                               ATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAT           >  SRR4302533.2764772/1‑94 (MQ=60)
                                                                                    GGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGT      <  SRR4302533.3553126/94‑1 (MQ=60)
                                                                                     GGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCTCCGTTATCCAGTT     <  SRR4302533.1420514/94‑1 (MQ=60)
                                                                                        GGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGC  <  SRR4302533.4874704/94‑1 (MQ=60)
                                                                                             |                                                                                        
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGC  >  NC_000913/3987729‑3987910

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