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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,473 +G intergenic (+40/‑193) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4721.G99.3% 430.1 / ‑3.8 140intergenic (+39/‑194)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (61/78);  total (62/78)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 143 total aligned reads displayed.

ACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAA‑‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGC  >  NC_000913/4008400‑4008558
                                                                               |                                                                                         
ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc                                                                            <  2:2551572/90‑1 (MQ=255)
 tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCa                                                                           >  1:1030894/3‑93 (MQ=255)
  gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                                                                          <  2:826713/92‑1 (MQ=255)
  gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                                                                          <  2:483481/92‑1 (MQ=255)
  gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                                                                          <  2:4494890/92‑1 (MQ=255)
  gggCTACAGCCTGAAGGATTGTATTGCGTTTGGCGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAg                                                                          >  2:4215796/2‑93 (MQ=255)
   ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                                                                         <  1:2506546/93‑1 (MQ=255)
   ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt                                                                         >  2:3217413/1‑93 (MQ=255)
     cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc                                                                       <  1:3340585/93‑2 (MQ=255)
     cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc                                                                       <  2:2953294/93‑2 (MQ=255)
     cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc                                                                       >  1:2035757/1‑92 (MQ=255)
      tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg                                                                      >  1:689379/1‑91 (MQ=255)
      tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg                                                                      >  2:2407127/1‑91 (MQ=255)
      tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg                                                                      >  2:2696585/1‑91 (MQ=255)
          gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc                                                                  >  1:4332377/3‑93 (MQ=255)
          gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc                                                                  >  2:4246988/3‑93 (MQ=255)
          gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcacc                                                                  >  2:4074519/3‑93 (MQ=255)
           ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCgcaccg                                                                 >  1:750969/2‑92 (MQ=255)
            cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                <  1:3892728/93‑1 (MQ=255)
            cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAg                                                                <  1:960598/93‑1 (MQ=255)
                aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                            <  1:2805984/93‑1 (MQ=255)
                aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                            <  2:3541119/93‑1 (MQ=255)
                aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                                            <  2:1727104/93‑1 (MQ=255)
                 ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                           >  2:399184/1‑93 (MQ=255)
                 ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                                           <  1:4531157/93‑1 (MQ=255)
                              tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                              >  2:1121649/5‑93 (MQ=255)
                              tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCCGCATTATGGGAAGTGCGCACCACCGTCTGAAAGGCCTTCAt                                              <  2:2832462/89‑1 (MQ=255)
                                tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                            <  2:3075612/91‑1 (MQ=255)
                                tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                            >  2:1871860/3‑93 (MQ=255)
                                tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                            <  1:1450524/91‑1 (MQ=255)
                                tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                            <  2:2683642/91‑1 (MQ=255)
                                   tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                         <  1:1791295/93‑1 (MQ=255)
                                   tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                         >  1:81779/1‑93 (MQ=255)
                                   ggACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                         <  1:2233640/92‑1 (MQ=255)
                                    gACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                        >  2:3324421/1‑93 (MQ=255)
                                     aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                                      <  2:110176/93‑1 (MQ=255)
                                       gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg                                   <  1:1271192/93‑1 (MQ=255)
                                        ggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                  >  2:3679479/1‑92 (MQ=255)
                                        ggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                  >  1:1064987/1‑92 (MQ=255)
                                        ggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg                                  >  1:4165025/1‑92 (MQ=255)
                                         gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                 <  1:3208493/93‑3 (MQ=255)
                                         gATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga                                 <  1:1244152/93‑3 (MQ=255)
                                          aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                <  2:439256/93‑4 (MQ=255)
                                          aTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa                                <  1:2648648/93‑4 (MQ=255)
                                           tGAACGACGCGGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                               <  1:3075847/93‑5 (MQ=255)
                                           tGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaag                               <  2:3992017/93‑5 (MQ=255)
                                              acgacgCCGAGATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             <  1:451557/93‑1 (MQ=255)
                                              acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             <  1:1118114/93‑1 (MQ=255)
                                              acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             >  2:2467002/1‑93 (MQ=255)
                                              acgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                             <  2:4116374/93‑1 (MQ=255)
                                               cgacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                            <  2:3949962/93‑1 (MQ=255)
                                                gccgcgGAAATGCTGTCGATGGGGGGGAAAGGCTGCAGTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  2:202191/91‑1 (MQ=255)
                                                gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  2:103893/93‑1 (MQ=255)
                                                gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           >  2:4158913/1‑93 (MQ=255)
                                                gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  2:1314234/93‑1 (MQ=255)
                                                gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           <  1:4441255/93‑1 (MQ=255)
                                                gacgCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                           >  2:2223405/1‑93 (MQ=255)
                                                 acgCCGAAATGCTCTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                          <  1:2678236/93‑1 (MQ=255)
                                                   gCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                        <  2:507341/93‑1 (MQ=255)
                                                     cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                      <  1:2277134/93‑1 (MQ=255)
                                                     cGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc                      >  2:196956/1‑93 (MQ=255)
                                                      gAAATGCTGTCGCTGGCGGGAAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCCCGCAGTGATTGGTACt                     >  2:841908/1‑93 (MQ=255)
                                                      gAAATGCTGGCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAACTTCATCCCGAGCTGGAAGGGATTGGTACt                     >  1:2236980/1‑93 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  1:890759/93‑1 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   >  1:3391435/1‑93 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  1:3768821/93‑1 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  1:3037406/93‑1 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  2:2085270/93‑1 (MQ=255)
                                                        aaTGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa                   <  2:2098657/93‑1 (MQ=255)
                                                           gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                <  2:2580693/93‑1 (MQ=255)
                                                           gCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                >  1:3024996/1‑93 (MQ=255)
                                                            cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATACCGAGCTGGAAGGGATTGGTACTTATGcc               >  1:3128615/1‑93 (MQ=255)
                                                            cTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc               <  1:1463743/93‑1 (MQ=255)
                                                             tGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg              >  1:1690261/1‑93 (MQ=255)
                                                              gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGTAGTGATTGGTACTAATGCcga             >  1:920091/1‑93 (MQ=255)
                                                              gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga             <  1:3254248/93‑1 (MQ=255)
                                                              gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga             <  1:2055819/93‑1 (MQ=255)
                                                              gTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga             <  2:1568400/93‑1 (MQ=255)
                                                               tCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac            <  1:4613932/93‑1 (MQ=255)
                                                                cGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg           <  1:378303/93‑1 (MQ=255)
                                                                 gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga          <  1:2502888/93‑1 (MQ=255)
                                                                 gATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga          <  2:2638556/93‑1 (MQ=255)
                                                                  aTGGCGGGGAAAGGCTGCATTATGGTCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac         <  2:1479920/93‑1 (MQ=255)
                                                                  aTGGCGGGGAAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac         >  2:430081/1‑93 (MQ=255)
                                                                  aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac         <  2:2255970/93‑1 (MQ=255)
                                                                  aTGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac         <  2:4315304/93‑1 (MQ=255)
                                                                     ggcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgc       <  1:3225253/89‑1 (MQ=255)
                                                                      gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg      >  1:4103495/4‑93 (MQ=255)
                                                                      gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg      <  1:1307376/90‑1 (MQ=255)
                                                                      gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg      >  1:1452411/4‑93 (MQ=255)
                                                                      gcggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAcgcg      <  2:1336706/90‑1 (MQ=255)
                                                                       cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg     <  2:4654945/91‑1 (MQ=255)
                                                                       cggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCgg     <  1:1228176/91‑1 (MQ=255)
                                                                        ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCggg    >  1:338655/2‑92 (MQ=255)
                                                                        ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    <  1:3579281/92‑1 (MQ=255)
                                                                        ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    <  2:3511676/92‑1 (MQ=255)
                                                                        ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    >  1:3730473/2‑93 (MQ=255)
                                                                        ggggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGt    <  1:3202892/92‑1 (MQ=255)
                                                                          ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc  >  2:1520282/1‑93 (MQ=255)
                                                                          ggAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGc  >  2:1678308/1‑93 (MQ=255)
                                                                               |                                                                                         
ACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGAAA‑‑GCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGC  >  NC_000913/4008400‑4008558

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGAAGGCTGGCGAACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  >  NC_000913/4008386‑4008565
                                                                                                |                                                                                             
gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                                                                                               <  SRR4302535.4233280/76‑1 (MQ=53)
gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                                                                                               <  SRR4302535.3360272/76‑1 (MQ=53)
       ccaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAT                                                                                         <  SRR4302535.2780208/83‑1 (MQ=49)
        gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                                                                                        >  SRR4302535.2628570/11‑94 (MQ=49)
         aagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTA                                                                                       >  SRR4302535.4576429/10‑94 (MQ=48)
          agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAT                                                                                      <  SRR4302535.2354039/86‑1 (MQ=48)
           gaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATG                                                                                     <  SRR4302535.1021142/87‑1 (MQ=49)
            aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                                                                                    <  SRR4302535.3344291/88‑1 (MQ=51)
            aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                                                                                    >  SRR4302535.1760338/7‑94 (MQ=51)
             AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                                                                                   <  SRR4302535.3071604/94‑1 (MQ=50)
               ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCA                                                                                 >  SRR4302535.1038215/4‑94 (MQ=57)
                tgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAG                                                                                <  SRR4302535.4002565/92‑1 (MQ=58)
                  GGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTG                                                                              <  SRR4302535.2524617/94‑1 (MQ=60)
                   GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGC                                                                             >  SRR4302535.2050421/1‑94 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                                                                            >  SRR4302535.694210/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                                                                            >  SRR4302535.1698782/1‑92 (MQ=60)
                    CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                                                                            <  SRR4302535.3363861/94‑3 (MQ=60)
                     TACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgc                                                                           >  SRR4302535.94673/1‑91 (MQ=60)
                       CAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgcac                                                                         >  SRR4302535.4034112/1‑89 (MQ=60)
                        agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC                                                                        >  SRR4302535.4362005/4‑94 (MQ=60)
                         gcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgcaccg                                                                       >  SRR4302535.756261/3‑87 (MQ=60)
                         CCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAG                                                                      >  SRR4302535.1201421/1‑94 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                     >  SRR4302535.338688/1‑94 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                     <  SRR4302535.3919415/94‑1 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                     <  SRR4302535.967423/94‑1 (MQ=60)
                           CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC                                                                     >  SRR4302535.1255636/1‑94 (MQ=60)
                            TGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCG                                                                    <  SRR4302535.2775887/94‑1 (MQ=57)
                               AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCT                                                                 <  SRR4302535.2825932/94‑1 (MQ=57)
                                GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTG                                                                <  SRR4302535.4562173/94‑1 (MQ=60)
                                 GATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGA                                                               >  SRR4302535.57736/1‑94 (MQ=60)
                                  ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA                                                              <  SRR4302535.912953/94‑1 (MQ=60)
                                  ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA                                                              <  SRR4302535.3356104/94‑1 (MQ=60)
                                    tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                            >  SRR4302535.2795015/14‑94 (MQ=60)
                                    tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                            <  SRR4302535.1707369/81‑1 (MQ=60)
                                    tgtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAG                                                            >  SRR4302535.2401960/14‑94 (MQ=60)
                                     gtattgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGA                                                           <  SRR4302535.2467044/82‑1 (MQ=60)
                                     TATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAC                                                          <  SRR4302535.582540/94‑1 (MQ=60)
                                       ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                         >  SRR4302535.237125/1‑94 (MQ=60)
                                       ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                         >  SRR4302535.2381732/1‑94 (MQ=60)
                                       ATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACC                                                         >  SRR4302535.635963/1‑94 (MQ=60)
                                        gtgcgtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302535.2084021/85‑1 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302535.2345907/94‑1 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        >  SRR4302535.2252083/1‑94 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        <  SRR4302535.2084022/94‑1 (MQ=60)
                                        TTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCT                                                        >  SRR4302535.1034242/1‑94 (MQ=60)
                                         TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                       <  SRR4302535.941555/94‑1 (MQ=60)
                                         TGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTT                                                       <  SRR4302535.992706/94‑1 (MQ=60)
                                           CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                     >  SRR4302535.1778076/1‑94 (MQ=60)
                                           CGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCA                                                     <  SRR4302535.3945934/94‑1 (MQ=60)
                                            gtttgGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAT                                                    <  SRR4302535.1862662/89‑1 (MQ=60)
                                             TG‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCC                                                 <  SRR4302535.1460901/94‑1 (MQ=60)
                                                 GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGA                                               >  SRR4302535.82362/1‑94 (MQ=60)
                                                  TGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                              <  SRR4302535.1804248/94‑1 (MQ=60)
                                                  GGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAG                                              <  SRR4302535.2249525/94‑1 (MQ=60)
                                                      GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                         >  SRR4302535.1072564/1‑94 (MQ=60)
                                                      GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                         <  SRR4302535.1280275/94‑1 (MQ=60)
                                                      GGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGG                                         >  SRR4302535.4193544/1‑94 (MQ=60)
                                                        GATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGgaa                                       <  SRR4302535.1253032/94‑4 (MQ=60)
                                                        GATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGgaa                                       <  SRR4302535.3230716/94‑4 (MQ=60)
                                                         ATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAG                                       <  SRR4302535.2667620/94‑1 (MQ=60)
                                                         ATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAG                                       >  SRR4302535.3196930/1‑94 (MQ=60)
                                                          TGAACGACGCGGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                      <  SRR4302535.3097308/94‑1 (MQ=60)
                                                          TGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                      >  SRR4302535.985089/1‑94 (MQ=60)
                                                          TGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGT                                      <  SRR4302535.3081714/94‑1 (MQ=60)
                                                            AACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGA                                    <  SRR4302535.3246209/94‑1 (MQ=60)
                                                             ACGACGCCGAGATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAT                                   <  SRR4302535.454719/94‑1 (MQ=60)
                                                             ACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAT                                   <  SRR4302535.1126090/94‑1 (MQ=60)
                                                               GACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                 >  SRR4302535.917468/1‑94 (MQ=60)
                                                               GACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTG                                 <  SRR4302535.4471707/94‑1 (MQ=60)
                                                                ACGCCGAAATGCTCTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGG                                <  SRR4302535.2697397/94‑1 (MQ=60)
                                                                   CCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                             >  SRR4302535.4202139/1‑94 (MQ=60)
                                                                   CCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAC                             >  SRR4302535.1400195/1‑94 (MQ=60)
                                                                    CGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACT                            <  SRR4302535.2293395/94‑1 (MQ=60)
                                                                     GAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTA                           >  SRR4302535.1820182/1‑94 (MQ=60)
                                                                      AAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                          >  SRR4302535.3415029/1‑94 (MQ=60)
                                                                      AAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAA                          >  SRR4302535.2344464/1‑94 (MQ=60)
                                                                       AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         >  SRR4302535.2651406/1‑94 (MQ=60)
                                                                       AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         <  SRR4302535.3794835/94‑1 (MQ=60)
                                                                       AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         <  SRR4302535.897072/94‑1 (MQ=60)
                                                                       AATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAT                         <  SRR4302535.3058682/94‑1 (MQ=60)
                                                                        ATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATG                        >  SRR4302535.2122399/1‑94 (MQ=60)
                                                                         TGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGC                       <  SRR4302535.359317/94‑1 (MQ=60)
                                                                         TGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGC                       >  SRR4302535.3046173/1‑94 (MQ=60)
                                                                          GCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGGGCACCAGCGTCTGAAAGACCTTCATACCGAGCTGGAAGGGATTGGTACTTATGCC                      >  SRR4302535.3150404/1‑94 (MQ=57)
                                                                           CTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                     <  SRR4302535.1474228/94‑1 (MQ=60)
                                                                           CTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCG                     >  SRR4302535.1702494/1‑94 (MQ=60)
                                                                            TGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGTAGTGATTGGTACTAATGCCGA                    >  SRR4302535.926616/1‑94 (MQ=60)
                                                                             GTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGAC                   <  SRR4302535.2070629/94‑1 (MQ=60)
                                                                             GTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGAC                   <  SRR4302535.3276858/94‑1 (MQ=60)
                                                                              TCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACG                  <  SRR4302535.4645279/94‑1 (MQ=60)
                                                                               CGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGA                 <  SRR4302535.380924/94‑1 (MQ=60)
                                                                                GATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGAC                <  SRR4302535.2520935/94‑1 (MQ=60)
                                                                                 ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG               <  SRR4302535.685392/94‑1 (MQ=60)
                                                                                 ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG               >  SRR4302535.811241/1‑94 (MQ=60)
                                                                                 ATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACG               <  SRR4302535.3783146/94‑1 (MQ=60)
                                                                                   GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG             >  SRR4302535.4131587/1‑94 (MQ=60)
                                                                                   GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG             >  SRR4302535.1462804/1‑94 (MQ=60)
                                                                                   GGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCG             <  SRR4302535.3247612/94‑1 (MQ=60)
                                                                                      gcggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGG            <  SRR4302535.1316714/84‑1 (MQ=60)
                                                                                       cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT           <  SRR4302535.801094/85‑1 (MQ=60)
                                                                                       cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT           <  SRR4302535.1236952/85‑1 (MQ=60)
                                                                                       cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGT           >  SRR4302535.3756306/10‑94 (MQ=60)
                                                                                       cggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGG           >  SRR4302535.341000/10‑94 (MQ=60)
                                                                                        ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG          <  SRR4302535.3225069/86‑1 (MQ=60)
                                                                                        ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG          <  SRR4302535.4492088/86‑1 (MQ=60)
                                                                                        ggggaaagGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTG          <  SRR4302535.3604091/86‑1 (MQ=60)
                                                                                        GG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCC        >  SRR4302535.3849761/1‑94 (MQ=60)
                                                                                         G‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCG       >  SRR4302535.3178349/1‑94 (MQ=60)
                                                                                            AAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302535.941217/1‑94 (MQ=60)
                                                                                            AAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCA     >  SRR4302535.3997677/1‑94 (MQ=60)
                                                                                              GGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT   <  SRR4302535.2366038/94‑1 (MQ=60)
                                                                                              GGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATT   <  SRR4302535.1830240/94‑1 (MQ=60)
                                                                                                GCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  <  SRR4302535.2522363/94‑1 (MQ=60)
                                                                                                |                                                                                             
TGGAAGGCTGGCGAACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTA  >  NC_000913/4008386‑4008565

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