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 884 257.2 12958910 96.9% 12557183 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.C93.7% 362.6 / 23.2 142intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (1/8);  new base C (10/123);  total (11/131)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
*NC_0009133,957,0912.C93.7% 363.6 / 23.2 142intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (1/8);  new base C (10/123);  total (11/131)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 158 total aligned reads displayed.

TGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957000‑3957173
                                                                                               ||                                                                                    
tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑c                                                                                     <  2:3286395/93‑1 (MQ=255)
   gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCCCCCCCCC‑‑A‑‑cccc                                                                                  >  2:5077826/1‑90 (MQ=255)
   gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccgt                                                                                  <  1:3349452/93‑4 (MQ=255)
   gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccct                                                                                 <  2:1184061/93‑1 (MQ=255)
      ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctctc                                                                              >  1:858430/1‑91 (MQ=255)
      ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctctc                                                                              <  2:3494109/93‑3 (MQ=255)
      ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctctc                                                                              <  2:610770/93‑3 (MQ=255)
      ggtagagaAGAAACTGCGCGACTATATGACAGATATGAAACGTATTGCTGTTGCGTGTTAAGTGCGCGCTGATGCCCTCGCCCCG‑‑A‑‑Ccctctc                                                                              <  2:2820609/93‑3 (MQ=37)
        taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccctctccc                                                                            <  2:5852981/93‑5 (MQ=255)
           gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcaca                                                                           <  2:4117146/93‑1 (MQ=38)
           gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcaca                                                                           <  1:4224652/93‑1 (MQ=38)
           gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcaca                                                                           <  2:3856588/93‑1 (MQ=38)
            taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAGGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACag                                                                          <  1:2387100/93‑1 (MQ=37)
             aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCcacacg                                                                         >  2:2073569/1‑91 (MQ=37)
             aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCTCCCCC‑‑A‑‑CACCCCcccacccc                                                                       >  1:386143/1‑89 (MQ=255)
              aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACaggg                                                                        <  1:2990561/93‑1 (MQ=37)
                 aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggaga                                                                     <  2:2604224/93‑1 (MQ=37)
                  aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagag                                                                    <  2:3153435/93‑1 (MQ=37)
                  aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagag                                                                    <  1:1636071/93‑1 (MQ=37)
                   cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagagg                                                                   <  1:354695/93‑1 (MQ=25)
                    tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                                  <  1:6093181/93‑1 (MQ=25)
                    tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                                  <  1:5449868/93‑1 (MQ=25)
                    tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagaggg                                                                  <  2:4080743/93‑1 (MQ=25)
                     gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCC‑‑CTCTCCCACagggagaggga                                                                 <  1:129307/93‑1 (MQ=25)
                      cgcgGCTATATGACAGATAGGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACagggagagggag                                                                <  2:1195481/93‑1 (MQ=25)
                       gctgCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAgggagagggaga                                                               <  1:1005045/90‑1 (MQ=21)
                       gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAgggagagggaga                                                               <  2:3815556/93‑1 (MQ=25)
                        cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaa                                                              <  2:1727305/93‑1 (MQ=25)
                         ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaa                                                             <  2:5410557/93‑1 (MQ=25)
                          gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGGGCGCGCTGATGCCCCTGTCACGG‑‑A‑‑CGTCTCCCACCCGGAGAGGGGGaaca                                                            >  1:2089165/1‑91 (MQ=255)
                           cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                                           >  1:3247110/1‑93 (MQ=21)
                           cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAac                                                           <  1:1351145/93‑1 (MQ=21)
                             atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                         <  1:4174166/93‑1 (MQ=21)
                              tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACt                                                        <  2:2831649/93‑1 (MQ=21)
                              tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACt                                                        <  1:2623633/93‑1 (MQ=21)
                               atGACAGATATGAAACGTATTGCTGTTGCGGGTTACGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                       <  2:4552696/93‑1 (MQ=21)
                                tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                      <  1:4823172/93‑1 (MQ=21)
                                 gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:1783949/93‑1 (MQ=21)
                                 gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:1013524/93‑1 (MQ=21)
                                 gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:1454246/93‑1 (MQ=21)
                                  aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:622719/93‑1 (MQ=21)
                                   cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                  >  2:3191891/1‑92 (MQ=21)
                                   cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                  <  1:2237195/93‑2 (MQ=21)
                                   cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCACG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                  <  2:4658949/93‑2 (MQ=21)
                                     gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  2:5481310/93‑4 (MQ=21)
                                      atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:1523584/93‑5 (MQ=21)
                                        atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:3153491/93‑1 (MQ=18)
                                        atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:3935850/93‑1 (MQ=18)
                                          gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  1:6083213/93‑1 (MQ=18)
                                          gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:1923959/93‑1 (MQ=18)
                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctca                                           <  1:2069487/93‑2 (MQ=18)
                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:4411706/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:5338660/93‑1 (MQ=18)
                                              cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  1:5554151/93‑1 (MQ=18)
                                              cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  2:6374374/93‑1 (MQ=21)
                                               gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:254129/93‑1 (MQ=18)
                                               gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  1:4455735/93‑1 (MQ=18)
                                                 aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  2:2040703/93‑1 (MQ=18)
                                                    gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  1:363040/93‑1 (MQ=18)
                                                      tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  2:1050033/93‑1 (MQ=18)
                                                       gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:4627217/93‑1 (MQ=18)
                                                        ttGCGGGTTAAGTGCGCGCTAATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  2:213704/93‑1 (MQ=18)
                                                         tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                             <  2:2910210/93‑1 (MQ=18)
                                                         tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                             <  2:3887439/93‑1 (MQ=18)
                                                          gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:4099594/93‑1 (MQ=18)
                                                           cGGGTTAAGTGCGCGCTGATGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  2:5291380/93‑1 (MQ=18)
                                                           cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  2:2130620/93‑1 (MQ=18)
                                                           cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  1:1302114/93‑1 (MQ=18)
                                                           cGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           >  2:1003920/1‑93 (MQ=18)
                                                            gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                          <  2:6110777/93‑1 (MQ=18)
                                                             ggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  2:3128559/93‑1 (MQ=18)
                                                              gTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        <  1:1165686/93‑1 (MQ=18)
                                                               ttAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  1:2583127/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:6138053/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:3732241/93‑1 (MQ=18)
                                                                 aaGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:2692296/93‑1 (MQ=18)
                                                                 aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACGCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:4512314/93‑1 (MQ=18)
                                                                  aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:5691133/93‑1 (MQ=18)
                                                                  aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:2604575/93‑1 (MQ=18)
                                                                   gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  1:3979240/93‑1 (MQ=18)
                                                                     gcgcgcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  1:3262793/93‑1 (MQ=21)
                                                                         gcTGATGCCGTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  1:1044885/93‑1 (MQ=21)
                                                                         gcTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  2:83492/93‑1 (MQ=21)
                                                                          cTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa            <  2:6451767/93‑1 (MQ=21)
                                                                           tGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  1:1848148/93‑1 (MQ=25)
                                                                           cgATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  2:3119394/92‑1 (MQ=25)
                                                                            gATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  1:555552/93‑1 (MQ=25)
                                                                             aTGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa         <  2:6324766/93‑1 (MQ=25)
                                                                              tGCCCTCACCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc        <  1:3453432/93‑1 (MQ=25)
                                                                              tGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc        <  2:760353/93‑1 (MQ=25)
                                                                               gCCCTCACCCCG‑‑ACCCTCTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc       >  2:5520604/1‑93 (MQ=25)
                                                                                 cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:500823/93‑1 (MQ=25)
                                                                                 cccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTCCTCCTGGAGAACGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:2821905/93‑1 (MQ=21)
                                                                                  ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  2:5280271/93‑1 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  2:5104223/93‑1 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  2:2698188/93‑1 (MQ=25)
                                                                                    tCACCCCG‑‑ACGCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTGGCGGTACAGCACCTTaa   <  2:3182267/93‑1 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  1:6225475/1‑93 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACAATCAAGGCCTTCTCCTGGAGAAAGCCTTGCTATTAGTTTCGGTACAGCACCTTaa   >  2:5642319/1‑93 (MQ=25)
                                                                                     cACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  2:2853495/93‑1 (MQ=37)
                                                                                               ||                                                                                    
TGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957000‑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
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  NC_000913/3957003‑3957184
                                                                                             |                                                                                              
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccgtc                                                                                            <  SRR4302534.3439134/94‑5 (MQ=60)
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CCCTC                                                                                            <  SRR4302534.5627625/94‑1 (MQ=60)
  AGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CCCTCTC                                                                                          >  SRR4302534.880103/1‑94 (MQ=60)
     TAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCA                                                                                         >  SRR4302534.6609810/1‑94 (MQ=60)
        GTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAG                                                                                      <  SRR4302534.4338528/94‑1 (MQ=60)
         TAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAGGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGG                                                                                     <  SRR4302534.2450531/94‑1 (MQ=60)
         TAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGG                                                                                     <  SRR4302534.6239112/94‑1 (MQ=60)
           AGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGA                                                                                   <  SRR4302534.3068744/94‑1 (MQ=60)
              AACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAG                                                                                <  SRR4302534.106716/94‑1 (MQ=60)
              AACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAG                                                                                <  SRR4302534.2819226/94‑1 (MQ=60)
               ACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGG                                                                               <  SRR4302534.3629060/94‑1 (MQ=60)
               ACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGG                                                                               <  SRR4302534.1678719/94‑1 (MQ=60)
               ACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGG                                                                               <  SRR4302534.4799384/94‑1 (MQ=60)
                CTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGG                                                                              <  SRR4302534.363768/94‑1 (MQ=60)
                 TGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGA                                                                             <  SRR4302534.4118529/94‑1 (MQ=60)
                 TGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGA                                                                             <  SRR4302534.5596709/94‑1 (MQ=60)
                 TGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGA                                                                             <  SRR4302534.6258873/94‑1 (MQ=60)
                  GCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAG                                                                            <  SRR4302534.132723/94‑1 (MQ=60)
                    GCTGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAA                                                                          <  SRR4302534.1030825/94‑1 (MQ=60)
                      GGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAA                                                                        <  SRR4302534.4217915/94‑1 (MQ=60)
                       GCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAAC                                                                       >  SRR4302534.3333575/1‑94 (MQ=60)
                        CTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACA                                                                      <  SRR4302534.1386771/94‑1 (MQ=60)
                          ATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACT                                                                    <  SRR4302534.4286918/94‑1 (MQ=60)
                           TATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTC                                                                   <  SRR4302534.2693038/94‑1 (MQ=60)
                            ATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCA                                                                 <  SRR4302534.4743327/94‑1 (MQ=60)
                             TGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAA                                                                <  SRR4302534.4954067/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302534.1039608/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302534.1412322/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302534.1492498/94‑1 (MQ=60)
                              GACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAG                                                              <  SRR4302534.1830277/94‑1 (MQ=60)
                                CAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGC                                                            <  SRR4302534.2296473/94‑1 (MQ=60)
                                  GATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCT                                                           <  SRR4302534.3967596/94‑1 (MQ=60)
                                  GATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCT                                                           <  SRR4302534.1756647/94‑1 (MQ=60)
                                   ATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAAAACTCAAaggctt                                                          >  SRR4302534.3363860/1‑88 (MQ=60)
                                    TATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTC                                                         <  SRR4302534.4169904/94‑1 (MQ=60)
                                    TATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACAGGGGGCGGGGGAAAACACTCAcggccttc                                                         >  SRR4302534.1040876/1‑86 (MQ=60)
                                     ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                        <  SRR4302534.4041275/94‑1 (MQ=60)
                                     ATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCT                                                        <  SRR4302534.3236534/94‑1 (MQ=60)
                                       GAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCC                                                      <  SRR4302534.6248630/94‑1 (MQ=60)
                                        AAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCT                                                     <  SRR4302534.3339940/94‑1 (MQ=55)
                                        AAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCAT                                                     <  SRR4302534.2123616/94‑1 (MQ=57)
                                          ACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGG                                                   <  SRR4302534.5482590/94‑1 (MQ=60)
                                          ACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGG                                                   <  SRR4302534.2654603/94‑1 (MQ=60)
                                           cttcTTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCCCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAATCACTCAAGGCCTTCTCCTGGA                                                  <  SRR4302534.6320036/90‑1 (MQ=32)
                                           CGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGA                                                  <  SRR4302534.5704094/94‑1 (MQ=60)
                                            GTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAG                                                 <  SRR4302534.3174935/94‑1 (MQ=60)
                                            GTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAG                                                 <  SRR4302534.4576179/94‑1 (MQ=60)
                                               TTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAG                                              <  SRR4302534.4457118/94‑1 (MQ=49)
                                                 GCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGC                                            <  SRR4302534.372300/94‑1 (MQ=60)
                                                 GCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGC                                            <  SRR4302534.1372389/94‑1 (MQ=60)
                                                   TGTTGCGGGTTAAGTGCGCGCTGCTGCCCTCGCCCCG‑‑ACCCTCTCCGCCAGGGAGAGGGAGACAACGCTCAAGGCGTTCTCCTGGAGAAGGCCT                                          <  SRR4302534.6293923/94‑1 (MQ=36)
                                                    GTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302534.3524881/94‑1 (MQ=35)
                                                    GTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTT                                         <  SRR4302534.4752764/94‑1 (MQ=60)
                                                      TGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCCCAGGGAGAGGGAGAAAAAACTCAAGGCCTTCTCCTGGAGAAGGCCTTGC                                       >  SRR4302534.3278612/1‑94 (MQ=49)
                                                       GCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCT                                      <  SRR4302534.4209813/94‑1 (MQ=60)
                                                        CGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTA                                     <  SRR4302534.5209461/94‑1 (MQ=60)
                                                        CGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTA                                     <  SRR4302534.1336455/94‑1 (MQ=60)
                                                        CGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTA                                     <  SRR4302534.643403/94‑1 (MQ=60)
                                                           GTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTA                                  <  SRR4302534.1196294/94‑1 (MQ=60)
                                                           GTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCTCTCCCACACGGAGAGGGAGAAAACACTCAAAGACTTCTCCTTGAGAAGGGCTTGCTATTA                                  >  SRR4302534.2162781/1‑94 (MQ=23)
                                                            TTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCCCACGGGGAGGGAGAAAACACTCAAAGCCTTCTCCTGGAGAAGGgcctgctatttt                                 >  SRR4302534.2187317/1‑82 (MQ=32)
                                                            TTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAG                                 <  SRR4302534.5807261/94‑1 (MQ=60)
                                                            TTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAG                                 <  SRR4302534.2651512/94‑1 (MQ=60)
                                                             TAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGT                                <  SRR4302534.6305008/94‑1 (MQ=60)
                                                             TAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGT                                <  SRR4302534.3832525/94‑1 (MQ=60)
                                                              AAGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT                               <  SRR4302534.2763435/94‑1 (MQ=60)
                                                              AAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACGCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT                               <  SRR4302534.4634527/94‑1 (MQ=49)
                                                              AAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTT                               <  SRR4302534.2729221/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302534.2341938/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302534.5845268/94‑1 (MQ=60)
                                                               AGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTG                              <  SRR4302534.2673493/94‑1 (MQ=60)
                                                                GTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGC                             <  SRR4302534.4085726/94‑1 (MQ=60)
                                                                  GCGCGCTGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGG                           <  SRR4302534.3349675/94‑1 (MQ=60)
                                                                      GCTGATGCCGTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACA                       <  SRR4302534.1071883/94‑1 (MQ=60)
                                                                        TGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGC                     <  SRR4302534.931591/94‑1 (MQ=60)
                                                                        TGATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGC                     <  SRR4302534.1896233/94‑1 (MQ=60)
                                                                         GATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCA                    <  SRR4302534.569246/94‑1 (MQ=60)
                                                                          ATGCCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAC                   <  SRR4302534.1262541/94‑1 (MQ=60)
                                                                           TGCCCTCACCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC                  <  SRR4302534.3546264/94‑1 (MQ=60)
                                                                             CCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTT                <  SRR4302534.513297/94‑1 (MQ=60)
                                                                             CCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTT                <  SRR4302534.3783942/94‑1 (MQ=60)
                                                                             CCCTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTCCTCCTGGAGAACGCCTTGCTATTAGTTGCGGTACAGCACCTT                <  SRR4302534.2896096/94‑1 (MQ=60)
                                                                               CTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAAGCCTTCTCCTGGAGAAGGGCTTGCTATTAGTTGCGGTACAGCACCTTAA              >  SRR4302534.6290108/1‑94 (MQ=60)
                                                                               CTCACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAA              >  SRR4302534.6394501/1‑94 (MQ=60)
                                                                                 CACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATG            >  SRR4302534.4058229/1‑94 (MQ=60)
                                                                                 CACCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATG            >  SRR4302534.1075222/1‑94 (MQ=60)
                                                                                   CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAT          <  SRR4302534.1644849/94‑1 (MQ=60)
                                                                                   CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAT          <  SRR4302534.4736628/94‑1 (MQ=60)
                                                                                   CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAT          >  SRR4302534.3133339/1‑94 (MQ=60)
                                                                                     CCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGT        <  SRR4302534.4086855/94‑1 (MQ=60)
                                                                                     CCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCCGGACAGCACCTTAATGATGT        >  SRR4302534.5661005/1‑94 (MQ=60)
                                                                                       ccgcccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGT        <  SRR4302534.2311525/88‑1 (MQ=60)
                                                                                        cgaccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTG       <  SRR4302534.4257880/89‑1 (MQ=60)
                                                                                         gaccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGA      <  SRR4302534.1871668/90‑1 (MQ=60)
                                                                                          accCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGAT     >  SRR4302534.312800/4‑94 (MQ=60)
                                                                                          C‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT   <  SRR4302534.6392646/94‑1 (MQ=60)
                                                                                          C‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT   >  SRR4302534.5777950/1‑94 (MQ=60)
                                                                                          C‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATAT   >  SRR4302534.2981984/1‑94 (MQ=60)
                                                                                           ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA    <  SRR4302534.4855107/92‑1 (MQ=60)
                                                                                           ccCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATA    >  SRR4302534.2650714/3‑94 (MQ=60)
                                                                                             CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  <  SRR4302534.6413475/94‑1 (MQ=60)
                                                                                             CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  SRR4302534.5725023/1‑94 (MQ=60)
                                                                                             |                                                                                              
GCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATC  >  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: