Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 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.0% 315.4 / 13.0 121intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base C (9/106);  total (9/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01
*NC_0009133,957,0912.C95.0% 313.6 / 13.0 121intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base C (9/106);  total (9/112)
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.37e-01

Only 100 of 131 total aligned reads displayed.

GAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957001‑3957173
                                                                                           ||                                                                                     
gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑cc                                                                                     <  2:24706/93‑2 (MQ=255)
    aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑Ccctct                                                                                 <  1:2235990/93‑2 (MQ=255)
     ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑Ccctctc                                                                                <  1:2051478/93‑3 (MQ=255)
      gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑Ccctctcc                                                                               <  2:83347/93‑4 (MQ=255)
      gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑Ccctctcc                                                                               <  1:5512205/93‑4 (MQ=255)
       taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCCA‑‑CCCTccccc                                                                              >  2:5306005/1‑91 (MQ=255)
          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                             <  1:1490220/93‑1 (MQ=38)
           taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                            <  1:3943828/93‑1 (MQ=37)
            aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                           <  1:5874305/93‑1 (MQ=37)
            aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                           <  1:1506081/93‑1 (MQ=37)
            aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                           <  1:5127257/93‑1 (MQ=37)
            aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                           <  1:442242/93‑1 (MQ=37)
              gAAGCTGCGAGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                         <  1:614449/93‑1 (MQ=25)
              gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                         <  2:1391639/93‑1 (MQ=37)
               aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                        <  2:308817/93‑1 (MQ=37)
                aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                       <  2:308317/93‑1 (MQ=37)
                aaCTGCGCGGCTATATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                       <  2:1602838/93‑1 (MQ=25)
                   tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCCTACCCCGAGCCTCTCCCACacggcgaggg                                                                    >  1:1180540/1‑93 (MQ=255)
                      gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAgggagagggaga                                                                 <  2:5863167/93‑1 (MQ=25)
                      gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAgggagagggaga                                                                 <  1:1462075/93‑1 (MQ=25)
                       cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCCACCCTCTCCCCCcgggagagggagga                                                                >  2:5349304/1‑91 (MQ=25)
                       cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaa                                                                <  1:5147220/93‑1 (MQ=25)
                        ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGaaa                                                               <  1:2054224/93‑1 (MQ=25)
                          cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAac                                                             <  1:3113233/93‑1 (MQ=21)
                          cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAac                                                             <  2:3949230/93‑1 (MQ=21)
                          cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACACGGGGGGGGGGaaaaa                                                             >  1:6123685/1‑92 (MQ=21)
                             tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCGCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACt                                                         >  1:4305824/1‑93 (MQ=21)
                              atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                        <  2:1701736/93‑1 (MQ=21)
                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                       <  1:5177571/93‑1 (MQ=21)
                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                       <  2:1476689/93‑1 (MQ=21)
                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                       <  1:5765629/93‑1 (MQ=21)
                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                       <  1:48730/93‑1 (MQ=21)
                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                       <  2:6072941/93‑1 (MQ=21)
                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                    <  1:5565873/93‑1 (MQ=21)
                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                    <  1:4575268/93‑1 (MQ=21)
                                   aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                  <  2:2422879/93‑3 (MQ=21)
                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGATAAAACACTCAAGGCCTtc                                              <  2:2436726/93‑1 (MQ=18)
                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:3465767/93‑1 (MQ=18)
                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:4706130/93‑1 (MQ=18)
                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:927910/93‑1 (MQ=18)
                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:83547/93‑1 (MQ=18)
                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:3726613/93‑1 (MQ=18)
                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:1435796/93‑1 (MQ=18)
                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:4978611/93‑1 (MQ=18)
                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  1:5707864/93‑1 (MQ=18)
                                          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGACTTCTcc                                           <  1:1415221/93‑1 (MQ=18)
                                           aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  2:2489644/93‑1 (MQ=18)
                                           aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  1:2148729/93‑1 (MQ=18)
                                            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:2488997/93‑1 (MQ=18)
                                            aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:4609262/93‑1 (MQ=18)
                                             cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCGCACAGGGAGAGGGAGAAAACAATCAAGGCCTTCTCCTgg                                        >  2:2718552/1‑93 (MQ=18)
                                              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:2782371/93‑1 (MQ=18)
                                              gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:2265983/93‑1 (MQ=18)
                                              gTATTGCTGTTGCGGGGTAAGTGCGCGCTGATGCCCTCCCCCCCACCCTCTCCCCCAAGGCGAAGGGGGAAACACCCAAACGCTTCTccctga                                       >  2:3382122/1‑89 (MQ=12)
                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  2:3572386/93‑1 (MQ=18)
                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  2:677933/93‑1 (MQ=18)
                                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  2:3302854/93‑1 (MQ=18)
                                                 ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4612686/93‑1 (MQ=18)
                                                  tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  1:1635564/93‑1 (MQ=18)
                                                  tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  1:4438388/93‑1 (MQ=18)
                                                   tcTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  2:3497629/92‑1 (MQ=18)
                                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCCCTCCCCCATGGAGGGGGGGGGAACACACCAAGCCGTCTCCCGGAGAAGGc                                 >  1:5374225/1‑93 (MQ=12)
                                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  1:4805151/93‑1 (MQ=18)
                                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  1:4017842/93‑1 (MQ=18)
                                                      gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:3425918/93‑1 (MQ=18)
                                                      gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:4595024/93‑1 (MQ=18)
                                                       ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:1001283/93‑1 (MQ=18)
                                                       ttGAGGGTTAAGTGCGCGCTGATGCCCTCGCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:3740307/93‑1 (MQ=18)
                                                         gCGTGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:6103709/93‑1 (MQ=18)
                                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:3264256/93‑1 (MQ=18)
                                                         gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:1641988/93‑1 (MQ=18)
                                                            ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  2:1441983/93‑1 (MQ=18)
                                                            ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  2:5205762/93‑1 (MQ=18)
                                                            ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  1:1657751/93‑1 (MQ=18)
                                                                aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:5340369/93‑1 (MQ=18)
                                                                aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:2676680/93‑1 (MQ=18)
                                                                aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:5975077/93‑1 (MQ=18)
                                                                 aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:2776490/93‑1 (MQ=18)
                                                                 aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:3664536/93‑1 (MQ=18)
                                                                 aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:335793/93‑1 (MQ=18)
                                                                  gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  1:441583/93‑1 (MQ=18)
                                                                   tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  1:4714255/93‑1 (MQ=21)
                                                                   tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:1113653/93‑1 (MQ=21)
                                                                   tGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCCGGAGAAGGCCTTGCTATTAGTTGc                  <  1:1554658/93‑1 (MQ=18)
                                                                    gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  2:797858/93‑1 (MQ=21)
                                                                      gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt               <  2:2590604/93‑1 (MQ=21)
                                                                       cgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa              <  2:3129113/93‑1 (MQ=21)
                                                                           gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  1:832/93‑1 (MQ=25)
                                                                              gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc       <  1:1580149/93‑1 (MQ=25)
                                                                               cccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  2:4872583/93‑1 (MQ=25)
                                                                                ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  2:3915710/93‑1 (MQ=25)
                                                                                 cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    <  1:2350572/93‑1 (MQ=37)
                                                                                 cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCCTGCTATTAGTTGCGGTACAGCACCTTa    >  2:136447/1‑93 (MQ=25)
                                                                                  tCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:1243656/1‑93 (MQ=25)
                                                                                  tCACCCCGACCCTCTCCCACAGGGAGCGGGAGGAAACACTCAAGGCCTTCTCCTTGCGCAGGCCTTGCTATTAGTTGTGGTACCGGGCCTTaa   >  2:1513959/1‑93 (MQ=16)
                                                                                  tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:5846674/1‑93 (MQ=37)
                                                                                  tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:1102475/1‑93 (MQ=37)
                                                                                  tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:1563701/1‑93 (MQ=37)
                                                                                  ccACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  2:4925545/92‑1 (MQ=37)
                                                                                   cACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  2:4189783/93‑1 (MQ=37)
                                                                                           ||                                                                                     
GAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957001‑3957173

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ 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
AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT  >  NC_000913/3957002‑3957183
                                                                                            |                                                                                            
gcgagtcacagAAGAAACTGCTCCGCTATATGACAGATATGGAACGTATTGCTGTTGCGGGGTAAGTGCGCCCTGAGGCCCTCTCCCCGA‑‑CCCT                                                                                           <  SRR4302532.941842/83‑1 (MQ=49)
   AGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CCCTCTC                                                                                        <  SRR4302532.2310182/94‑1 (MQ=60)
    GGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCC                                                                                         <  SRR4302532.2119792/94‑1 (MQ=60)
     GTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCC                                                                                        <  SRR4302532.5703483/94‑1 (MQ=60)
       ggggAAGAAACTTCGCGGCTATATGACAGATAAGAAACGTATTGCTGTTGCGGGTTAAGTTCGCGCTGATGCCCTCCCCCCGACCCTCTCCCAC                                                                                      <  SRR4302532.5058144/90‑1 (MQ=60)
       AGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCAC                                                                                      <  SRR4302532.693345/94‑1 (MQ=60)
         GTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAG                                                                                    <  SRR4302532.1540495/94‑1 (MQ=60)
          TAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGG                                                                                   <  SRR4302532.4076723/94‑1 (MQ=60)
           AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGG                                                                                  <  SRR4302532.5303639/94‑1 (MQ=60)
           AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGG                                                                                  <  SRR4302532.1556917/94‑1 (MQ=60)
           AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGG                                                                                  <  SRR4302532.6077845/94‑1 (MQ=60)
           AAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGG                                                                                  <  SRR4302532.457539/94‑1 (MQ=60)
             GAAGCTGCGAGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAG                                                                                <  SRR4302532.635869/94‑1 (MQ=60)
                 CTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCCTACCCCGAGCCTCTCCCACACGGCGAGGG                                                                            >  SRR4302532.1219978/1‑94 (MQ=60)
                  TGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGA                                                                           <  SRR4302532.1615054/94‑1 (MQ=60)
                     GCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAA                                                                        <  SRR4302532.1511343/94‑1 (MQ=60)
                      CGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAA                                                                       <  SRR4302532.5324292/94‑1 (MQ=60)
                       GGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGAAAA                                                                      <  SRR4302532.2122631/94‑1 (MQ=60)
                       GGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCCCTCCCACAGGGGGAGGGAGAAAA                                                                      >  SRR4302532.5431759/1‑94 (MQ=60)
                        GCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAcggggggggggaaaaa                                                                     >  SRR4302532.6335004/1‑78 (MQ=60)
                         CTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACA                                                                    <  SRR4302532.3216522/94‑1 (MQ=60)
                           ATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCGCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACT                                                                  >  SRR4302532.4451411/1‑94 (MQ=60)
                             ATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCA                                                                <  SRR4302532.1529299/94‑1 (MQ=60)
                              TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                               <  SRR4302532.50451/94‑1 (MQ=60)
                              TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                               <  SRR4302532.5966259/94‑1 (MQ=60)
                              TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                               <  SRR4302532.4435097/94‑1 (MQ=60)
                              TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                               <  SRR4302532.5355871/94‑1 (MQ=60)
                               GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAG                                                             <  SRR4302532.1915310/94‑1 (MQ=60)
                               GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                             <  SRR4302532.5812098/94‑1 (MQ=60)
                                ACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGG                                                            <  SRR4302532.5759144/94‑1 (MQ=60)
                                ACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGG                                                            <  SRR4302532.4731294/94‑1 (MQ=60)
                                 CAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCGACCCGCTCCCACAGGGAGAGGGAGAAAACACTCAAGGC                                                           <  SRR4302532.3086777/94‑1 (MQ=60)
                                     TATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTC                                                        <  SRR4302532.3253024/94‑1 (MQ=60)
                                     TATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGATGGACTCACCCCGA‑‑Ccccctcccacagcgagagagcgagaaaactccagcacctc                                                      >  SRR4302532.1684079/1‑54 (MQ=54)
                                      ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                       <  SRR4302532.3822783/94‑1 (MQ=60)
                                      ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                       <  SRR4302532.959695/94‑1 (MQ=60)
                                        GAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCC                                                     <  SRR4302532.5906851/94‑1 (MQ=60)
                                         AAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGACTTCTCCT                                                    <  SRR4302532.1462704/94‑1 (MQ=55)
                                          AACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTG                                                   <  SRR4302532.2220533/94‑1 (MQ=60)
                                                TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                             <  SRR4302532.4769970/94‑1 (MQ=60)
                                                 TGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGG                                            <  SRR4302532.1690928/94‑1 (MQ=60)
                                                 TGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGG                                            <  SRR4302532.4589634/94‑1 (MQ=60)
                                                  GCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑Ccccctcccccatggagggggggggaacacaccaagccgtctcccggagaaggc                                         >  SRR4302532.5560417/1‑41 (MQ=60)
                                                    TGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCT                                         <  SRR4302532.4153232/94‑1 (MQ=60)
                                                    TGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCT                                         <  SRR4302532.4970628/94‑1 (MQ=60)
                                                     GTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                        <  SRR4302532.4751670/94‑1 (MQ=60)
                                                      TTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                       <  SRR4302532.1035368/94‑1 (MQ=60)
                                                      TTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                       <  SRR4302532.4797533/94‑1 (MQ=60)
                                                      TTGCGGGTTAAGTGCGCCCTGACGCCCCCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                       <  SRR4302532.4579156/94‑1 (MQ=55)
                                                      TTGAGGGTTAAGTGCGCGCTGATGCCCTCGCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTG                                       <  SRR4302532.3865998/94‑1 (MQ=50)
                                                        GCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCT                                     <  SRR4302532.3372812/94‑1 (MQ=60)
                                                           GGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATT                                  <  SRR4302532.1713931/94‑1 (MQ=60)
                                                            GTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCTCCCCCAGGGAGGGGGGGAAAACACTCAAGGACTTCTCCTGGAGAAGGCCTTGCTATTT                                 >  SRR4302532.6263290/1‑94 (MQ=39)
                                                               AAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT                              <  SRR4302532.2766509/94‑1 (MQ=60)
                                                                AGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                             <  SRR4302532.2869667/94‑1 (MQ=60)
                                                                AGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                             <  SRR4302532.3787605/94‑1 (MQ=60)
                                                                 GTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGC                            <  SRR4302532.456855/94‑1 (MQ=60)
                                                                  TGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCG                           <  SRR4302532.4875107/94‑1 (MQ=60)
                                                                  TGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCCGGAGAAGGCCTTGCTATTAGTTGCG                           <  SRR4302532.1607229/94‑1 (MQ=52)
                                                                         TGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGC                    <  SRR4302532.5188497/94‑1 (MQ=60)
                                                                          GATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCA                   <  SRR4302532.871/94‑1 (MQ=60)
                                                                             GCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCT                <  SRR4302532.1633654/94‑1 (MQ=60)
                                                                             GCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCT                <  SRR4302532.2221497/94‑1 (MQ=60)
                                                                                CTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAA             <  SRR4302532.2428994/94‑1 (MQ=60)
                                                                                 tccccccgcccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT            <  SRR4302532.1360927/83‑1 (MQ=60)
                                                                                   ACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA          >  SRR4302532.4383038/1‑94 (MQ=60)
                                                                                   ACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGA          >  SRR4302532.2292258/1‑94 (MQ=60)
                                                                                     CCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATG        >  SRR4302532.1800449/1‑94 (MQ=60)
                                                                                       cgaccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG      <  SRR4302532.1213051/89‑1 (MQ=60)
                                                                                       cgaccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG      >  SRR4302532.614219/6‑94 (MQ=60)
                                                                                        gaccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA     >  SRR4302532.2396082/5‑94 (MQ=60)
                                                                                         C‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT  >  SRR4302532.2093808/1‑94 (MQ=60)
                                                                                         C‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT  <  SRR4302532.433881/94‑1 (MQ=60)
                                                                                          ccCTCTCCCACAGGGAGAGGGAGATAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA   >  SRR4302532.443944/3‑94 (MQ=60)
                                                                                          ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGTGGTACAGCACCTTAATGATGTGATA   >  SRR4302532.4268279/3‑94 (MQ=60)
                                                                                          ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA   >  SRR4302532.2926230/3‑94 (MQ=60)
                                                                                          ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA   <  SRR4302532.3297974/92‑1 (MQ=60)
                                                                                          ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCCTAATGATGTGATA   <  SRR4302532.3342236/92‑1 (MQ=60)
                                                                                            |                                                                                            
AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGA‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT  >  NC_000913/3957002‑3957183

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