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 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,092 +CC intergenic (+26/+58) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,0911.C95.4% 401.8 / 16.1 153intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base C (5/141);  total (5/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01
*NC_0009133,957,0912.C95.4% 402.6 / 16.1 153intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base C (5/141);  total (5/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 167 total aligned reads displayed.

TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173
                                                                                               ||                                                                                       
 tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑c                                                                                        <  2:3254375/93‑1 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:1551755/93‑3 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:3112562/93‑3 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:5164218/93‑3 (MQ=255)
    gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc‑‑‑t                                                                                  <  1:703850/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  1:1021298/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  2:809937/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  1:1555842/93‑1 (MQ=255)
            gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCcaca                                                                            <  2:408603/93‑1 (MQ=38)
             taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACag                                                                           <  1:2864458/93‑1 (MQ=37)
              aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagg                                                                          <  1:2420264/93‑1 (MQ=37)
               aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggg                                                                         <  1:3184175/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  1:3927307/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:6031079/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:5567110/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:4771747/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:2121284/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagag                                                                     <  1:80414/93‑1 (MQ=37)
                    cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagg                                                                    <  2:1259763/93‑1 (MQ=25)
                    cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagg                                                                    <  1:4563853/93‑1 (MQ=25)
                     tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagaggg                                                                   <  1:152831/93‑1 (MQ=25)
                      gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagaggga                                                                  <  1:177466/93‑1 (MQ=25)
                       cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagggag                                                                 <  2:3759527/93‑1 (MQ=25)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGaaaa                                                             <  1:5135331/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCCC‑‑‑TCCCACAGGGAGCGGGAGaaaa                                                             >  2:3279794/1‑93 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAaca                                                           <  1:5849000/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCC‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAaca                                                           >  2:5469641/1‑93 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCAACCCGACA‑‑CTC‑‑‑TCCCCCCAGGGGAGGGGGAAAaca                                                           >  1:3343323/1‑93 (MQ=25)
                              atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAacac                                                          <  2:2783006/93‑1 (MQ=21)
                                atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTc                                                        <  2:1174069/93‑1 (MQ=21)
                                atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTc                                                        <  2:2793655/93‑1 (MQ=21)
                                  gACAGATATGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑CGC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:3294301/93‑1 (MQ=18)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:1806617/93‑1 (MQ=21)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:2508504/93‑1 (MQ=21)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:1904400/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:1973059/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:397507/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:3033432/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:675265/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:5979940/93‑1 (MQ=21)
                                    cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                  <  1:5549775/93‑2 (MQ=21)
                                     aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                 <  2:1443408/93‑3 (MQ=21)
                                      gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  1:9338/93‑4 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:5607446/93‑5 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAAGGAGAAAACACTCAAggcct                                               <  2:6649291/93‑5 (MQ=21)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:2687964/93‑1 (MQ=18)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:2281234/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:392906/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:4103406/93‑1 (MQ=18)
                                             accGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  1:4575859/91‑1 (MQ=18)
                                             aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  2:4929641/93‑1 (MQ=18)
                                              aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:3706130/93‑1 (MQ=18)
                                              aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:3890355/93‑1 (MQ=21)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  1:2062873/93‑1 (MQ=18)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  1:5716829/93‑1 (MQ=18)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  2:624573/93‑1 (MQ=18)
                                                gTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:399765/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCTCACCGCGAGAGGGAGAAAAAACGCAAAGCCTTCTCCTGGAGaa                                    >  2:4808822/1‑93 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:3736839/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:5425655/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:2692190/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4010176/93‑1 (MQ=18)
                                                    tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  2:888500/93‑1 (MQ=18)
                                                      cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                 <  2:5464697/93‑1 (MQ=18)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  1:1270206/93‑1 (MQ=18)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  2:4740879/93‑1 (MQ=18)
                                                        tttGCGGGTTAAGTTCGCGCTGATCCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:552755/92‑1 (MQ=18)
                                                        gTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:6429211/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:5530893/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:4337576/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:4950168/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  2:5130233/93‑1 (MQ=18)
                                                          tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                             <  2:1958098/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:657926/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:3748690/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:5063741/93‑1 (MQ=18)
                                                            cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  2:3034061/93‑1 (MQ=18)
                                                             gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                          <  1:4431487/93‑1 (MQ=18)
                                                              ggTTAAGTGCGCGCTGCTGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  1:398457/93‑1 (MQ=18)
                                                              ggTTAAGTGCGCGCTGATGCCCCCCCCCCC‑‑C‑‑CCCCCTTCCCCCGGGGGGGGGGGGAAAAAACCCAAGCCCTCCCCTGGGGAAGGCCCTGCttt                         >  2:4506482/1‑91 (MQ=255)
                                                               gTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        <  2:2273293/93‑1 (MQ=18)
                                                                 tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  2:5187312/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:3226459/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:1487841/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:2729588/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:5071597/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:5640282/93‑1 (MQ=18)
                                                                    gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  2:1970743/93‑1 (MQ=18)
                                                                     tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:4398858/93‑1 (MQ=21)
                                                                      gcgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  1:2529276/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  2:6278852/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  1:6403099/93‑1 (MQ=18)
                                                                            tGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  1:2980254/93‑1 (MQ=25)
                                                                             gATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:3058079/93‑1 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGGGAGGGAGAAAACACTCAAGGCCTTTTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     >  2:3315625/1‑93 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  1:5960311/93‑1 (MQ=25)
                                                                                   cTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    >  1:6513542/1‑93 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:1019747/93‑1 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:2986405/93‑1 (MQ=37)
                                                                                     cACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:79103/93‑1 (MQ=37)
                                                                                               ||                                                                                       
TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173

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
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  NC_000913/3957003‑3957184
                                                                                             |                                                                                                
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑C‑‑CCTC                                                                                            <  SRR4302547.716857/94‑1 (MQ=60)
 CAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctct                                                                                             <  SRR4302547.207446/94‑6 (MQ=60)
 CAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctct                                                                                             <  SRR4302547.1910553/94‑6 (MQ=60)
 CAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctct                                                                                             <  SRR4302547.1041322/94‑6 (MQ=60)
 CAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctct                                                                                             <  SRR4302547.1586798/94‑6 (MQ=60)
         TAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGG                                                                                     <  SRR4302547.2921675/94‑1 (MQ=60)
         TAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGG                                                                                     <  SRR4302547.3115930/94‑1 (MQ=60)
          AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGG                                                                                    <  SRR4302547.2244094/94‑1 (MQ=60)
          AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGG                                                                                    <  SRR4302547.2469256/94‑1 (MQ=60)
          AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGG                                                                                    <  SRR4302547.5842267/94‑1 (MQ=60)
           AGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGA                                                                                   <  SRR4302547.6431814/94‑1 (MQ=60)
           AGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGA                                                                                   <  SRR4302547.3247628/94‑1 (MQ=60)
            GAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCC‑‑TCTCCCACAGGGAG                                                                                  <  SRR4302547.4007185/94‑1 (MQ=60)
            GAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAG                                                                                  <  SRR4302547.5590625/94‑1 (MQ=60)
               ACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGG                                                                               <  SRR4302547.81916/94‑1 (MQ=60)
                CTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGG                                                                              <  SRR4302547.4658274/94‑1 (MQ=60)
                CTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGG                                                                              <  SRR4302547.5163113/94‑1 (MQ=60)
                 TGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGA                                                                             <  SRR4302547.155619/94‑1 (MQ=60)
                  GCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAG                                                                            <  SRR4302547.180765/94‑1 (MQ=60)
                    GCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAA                                                                          <  SRR4302547.2579981/94‑1 (MQ=60)
                       GCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAAC                                                                       <  SRR4302547.5241135/94‑1 (MQ=60)
                         TATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACAC                                                                     <  SRR4302547.5969985/94‑1 (MQ=60)
                           TATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTC                                                                   <  SRR4302547.6019613/94‑1 (MQ=60)
                             TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                                <  SRR4302547.4305938/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302547.2559219/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302547.1942077/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302547.5667212/94‑1 (MQ=60)
                               ACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGG                                                             <  SRR4302547.2012231/94‑1 (MQ=60)
                               ACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGG                                                             <  SRR4302547.3093757/94‑1 (MQ=60)
                               ACAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGG                                                             <  SRR4302547.6103982/94‑1 (MQ=60)
                                CAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGC                                                            <  SRR4302547.5664284/94‑1 (MQ=60)
                                CAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccccctcccacacggagggggggaaaaaaaccaaggc                                                              >  SRR4302547.2224033/1‑58 (MQ=60)
                                  GATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCT                                                           <  SRR4302547.9510/94‑1 (MQ=60)
                                   ATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTT                                                          <  SRR4302547.3312236/94‑1 (MQ=60)
                                   ATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTT                                                          <  SRR4302547.5003634/94‑1 (MQ=60)
                                    TATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTC                                                         <  SRR4302547.1003472/94‑1 (MQ=60)
                                    TATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTC                                                         <  SRR4302547.1509883/94‑1 (MQ=60)
                                     atgaactgcattgatgtcggcaGTTAAGTGCGCGCTGATTCCCTGCCCCCG‑‑C‑‑CCGTCGCCCACAGGGAGAGGGAGAAAACACTCAAGgccttat                                                       <  SRR4302547.6259867/72‑8 (MQ=11)
                                     ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                        <  SRR4302547.2741917/94‑1 (MQ=60)
                                     ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                        <  SRR4302547.1328040/94‑1 (MQ=60)
                                        AAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCT                                                     <  SRR4302547.3877077/94‑1 (MQ=60)
                                         ACCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTG                                                    <  SRR4302547.4670574/94‑1 (MQ=42)
                                           CGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGA                                                  <  SRR4302547.5834913/94‑1 (MQ=60)
                                           CGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGA                                                  <  SRR4302547.2104034/94‑1 (MQ=60)
                                               TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                              <  SRR4302547.2746234/94‑1 (MQ=60)
                                               TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                              <  SRR4302547.4091688/94‑1 (MQ=60)
                                               TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                              <  SRR4302547.5537858/94‑1 (MQ=60)
                                               TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                              <  SRR4302547.3813136/94‑1 (MQ=60)
                                                   TGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCT                                          <  SRR4302547.1295783/94‑1 (MQ=60)
                                                    TTTGCGGGTTAAGTTCGCGCTGATCCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.562715/94‑1 (MQ=49)
                                                    GTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.6563239/94‑1 (MQ=49)
                                                    GTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.5645052/94‑1 (MQ=60)
                                                    GTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.3725099/94‑1 (MQ=60)
                                                    GTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.2586554/94‑1 (MQ=60)
                                                    GTCGCGGGTTAAGTGCGCCCTGATGCCCTCCCCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.1102190/94‑1 (MQ=38)
                                                    CTTGCGGGTTAAGTGCGCGCTGATGCTCTCACCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302547.6546601/94‑1 (MQ=45)
                                                     TTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                        <  SRR4302547.5052664/94‑1 (MQ=60)
                                                     TTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                        <  SRR4302547.6703651/94‑1 (MQ=60)
                                                       GCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCT                                      <  SRR4302547.3825212/94‑1 (MQ=60)
                                                       GCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCC‑‑TCTCCCCCCGGGGGAGGGGGAAAACACCCAAAGCCTTCTCCTCGAGAAGGGCCTGCT                                      >  SRR4302547.2248669/1‑94 (MQ=15)
                                                         GGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCC‑‑C‑‑TCTCCCCCAGGGGGAGGGAGAAAACACTCcaaagcttctcctctaggagggcctgctat                                     >  SRR4302547.2255933/1‑64 (MQ=0)
                                                         GGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAT                                    <  SRR4302547.4523089/94‑1 (MQ=60)
                                                          GGTTAAGTGCGCGCTGCTGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATT                                   <  SRR4302547.405868/94‑1 (MQ=60)
                                                          GGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATT                                   <  SRR4302547.895682/94‑1 (MQ=60)
                                                            TTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAG                                 <  SRR4302547.4559902/94‑1 (MQ=60)
                                                             TAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACGGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGT                                <  SRR4302547.4027872/94‑1 (MQ=60)
                                                              AAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT                               <  SRR4302547.3290853/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302547.5756713/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302547.5663915/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302547.5176235/94‑1 (MQ=60)
                                                                GTGCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGC                             <  SRR4302547.731518/94‑1 (MQ=60)
                                                                  GCGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGG                           <  SRR4302547.2580394/94‑1 (MQ=60)
                                                                   CGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGT                          <  SRR4302547.5435538/94‑1 (MQ=60)
                                                                   CGCGCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGT                          <  SRR4302547.6536655/94‑1 (MQ=56)
                                                                      GCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACA                       <  SRR4302547.2958970/94‑1 (MQ=60)
                                                                      GCTGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACA                       <  SRR4302547.1793472/94‑1 (MQ=60)
                                                                        TGATGCCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGC                     <  SRR4302547.3039555/94‑1 (MQ=60)
                                                                             CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTT                <  SRR4302547.282125/94‑1 (MQ=60)
                                                                              CCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTA               >  SRR4302547.6649136/1‑94 (MQ=60)
                                                                              CCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTA               <  SRR4302547.6083903/94‑1 (MQ=60)
                                                                               CTCACCCCG‑‑ACCC‑‑TCTCCCCCAGGGGGAGGGAGAAAACACTCAAGGCCTTCTTCTGGAGAAAGCCTTGTTAATAGTTGCcgcaaagacccctta              >  SRR4302547.1553319/1‑79 (MQ=20)
                                                                                TCACCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT             <  SRR4302547.2674732/94‑1 (MQ=60)
                                                                                TCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT             <  SRR4302547.3045826/94‑1 (MQ=60)
                                                                                TCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT             <  SRR4302547.1039732/94‑1 (MQ=60)
                                                                                 CACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTGGTTGCGGTACAGCACCTTAATG            >  SRR4302547.2089307/1‑94 (MQ=60)
                                                                                 CACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATG            <  SRR4302547.80582/94‑1 (MQ=60)
                                                                                  ACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA           >  SRR4302547.4329731/1‑94 (MQ=60)
                                                                                  ACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA           >  SRR4302547.3970247/1‑94 (MQ=60)
                                                                                    cccccgaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA           <  SRR4302547.6273102/85‑1 (MQ=60)
                                                                                    cccccgaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA           <  SRR4302547.6488127/85‑1 (MQ=60)
                                                                                    CCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATG         >  SRR4302547.6295771/1‑94 (MQ=60)
                                                                                    CCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATG         >  SRR4302547.5238099/1‑94 (MQ=60)
                                                                                     CCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGT        <  SRR4302547.4753994/94‑1 (MQ=60)
                                                                                        cgaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG       <  SRR4302547.5862786/89‑1 (MQ=60)
                                                                                        cgaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG       <  SRR4302547.3942351/89‑1 (MQ=60)
                                                                                        cgaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG       <  SRR4302547.4000347/89‑1 (MQ=60)
                                                                                         gaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      <  SRR4302547.1039347/90‑1 (MQ=60)
                                                                                         gaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      <  SRR4302547.642216/90‑1 (MQ=60)
                                                                                         gaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      >  SRR4302547.6349048/5‑94 (MQ=60)
                                                                                         gaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      <  SRR4302547.3791917/90‑1 (MQ=60)
                                                                                         gaccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      <  SRR4302547.2105085/90‑1 (MQ=60)
                                                                                          accC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGAT     <  SRR4302547.6687241/91‑1 (MQ=60)
                                                                                          C‑‑C‑‑TCTCCCACCGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGACAATGCCTAGCTATTAGTTGCGGTACTGCACCTTAATGATGTGATAT   >  SRR4302547.1416436/1‑94 (MQ=60)
                                                                                          C‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT   >  SRR4302547.2650681/1‑94 (MQ=60)
                                                                                           ccC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA    >  SRR4302547.5147919/3‑94 (MQ=60)
                                                                                             C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  SRR4302547.4015008/1‑94 (MQ=60)
                                                                                             C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  SRR4302547.1288827/1‑94 (MQ=60)
                                                                                             C‑‑TCTCCAACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  SRR4302547.4293631/1‑94 (MQ=60)
                                                                                             |                                                                                                
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  NC_000913/3957003‑3957184

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