Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 1548 369.5 20095359 96.9% 19472402 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,890,225 C→T R345R (CGG→CGA b1808 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,890,2250CT99.5% 731.2 / ‑5.5 202R345R (CGG→CGAb1808hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (113/88);  total (114/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 205 total aligned reads displayed.

GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314
                                                                                          |                                                                                         
gAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCg                                                                                         >  1:5776347/1‑93 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  2:7587318/93‑1 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        >  2:570980/1‑93 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                        <  1:7552995/93‑1 (MQ=255)
 aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCgt                                                                                        >  1:378042/1‑93 (MQ=255)
  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       >  2:404390/1‑93 (MQ=255)
  gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                       >  2:8724111/1‑93 (MQ=255)
   cggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtgt                                                                                      <  2:7340268/93‑1 (MQ=255)
    ggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGtgtg                                                                                     <  1:2267239/93‑1 (MQ=255)
       gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                  >  2:4755469/1‑93 (MQ=255)
        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 >  1:9116621/1‑93 (MQ=255)
        tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                 >  1:6891262/1‑93 (MQ=255)
         gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                                <  1:7288458/93‑1 (MQ=255)
          agagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAc                                                                               >  1:4460296/1‑93 (MQ=255)
           gagCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACa                                                                              <  1:806768/93‑1 (MQ=255)
            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                             <  1:8452748/93‑1 (MQ=255)
            agCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACaa                                                                             <  2:3937215/93‑1 (MQ=255)
             gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                            <  2:1492702/93‑1 (MQ=255)
              ccAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACg                                                                           >  2:1661392/1‑93 (MQ=255)
               cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTCCAGTCGTGTGCGATACAACGt                                                                          <  2:6448126/93‑1 (MQ=255)
                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:6158506/93‑1 (MQ=255)
                agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                         <  2:852753/93‑1 (MQ=255)
                  agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                       >  2:8260482/1‑93 (MQ=255)
                    tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                     >  2:2085864/1‑93 (MQ=255)
                        aTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCta                                                                 >  1:2656873/1‑92 (MQ=255)
                        aTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTc                                                                 <  1:9578406/93‑1 (MQ=255)
                          ggCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaa                                                               >  1:4014159/1‑93 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:3058919/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  1:717241/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:6647653/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:4961810/93‑1 (MQ=255)
                           gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                              <  2:7379966/93‑1 (MQ=255)
                            cgcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAt                                                             <  2:3820219/93‑1 (MQ=255)
                               gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                          >  1:9346212/1‑93 (MQ=255)
                                  gTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCa                                                       <  1:5966714/93‑1 (MQ=255)
                                   tGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTTTGCGATACAACGTGGCGCGCTCAAATGCCGCAt                                                      <  2:3650967/93‑1 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     >  1:9631697/1‑93 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     <  2:3005242/93‑1 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     <  2:9342144/93‑1 (MQ=255)
                                    gCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATc                                                     >  1:5558101/1‑93 (MQ=255)
                                       ttCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAg                                                  >  2:7009157/1‑93 (MQ=255)
                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                <  2:4358212/93‑1 (MQ=255)
                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                <  2:9899023/93‑1 (MQ=255)
                                         cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                                >  2:5589815/1‑93 (MQ=255)
                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                             >  1:9776467/1‑93 (MQ=255)
                                            aCCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgg                                             >  1:5450348/1‑93 (MQ=255)
                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAgcc                                            >  1:9242203/1‑91 (MQ=255)
                                             ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                            <  1:3583627/93‑1 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  1:7439613/1‑93 (MQ=255)
                                              cAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCg                                           >  1:4255063/1‑93 (MQ=255)
                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTCGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGGTCAGATGCCGCATCCTGCATGGccg                                          >  2:7038345/1‑91 (MQ=255)
                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                          >  1:1683057/1‑93 (MQ=255)
                                               aGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCgg                                          >  2:7038343/1‑93 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  2:1115207/93‑1 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:1204015/93‑1 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        >  1:7944362/1‑93 (MQ=255)
                                                 tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGCGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:9873459/93‑1 (MQ=255)
                                                 gatcTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                        <  1:5868897/90‑1 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTTCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  1:5384769/1‑93 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  2:2403205/1‑93 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  2:7231543/93‑1 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  1:445134/93‑1 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   <  2:6118344/93‑1 (MQ=255)
                                                      gTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGccc                                   >  1:8045834/1‑93 (MQ=255)
                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGCTACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  1:193884/1‑93 (MQ=255)
                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  2:3976936/1‑93 (MQ=255)
                                                       tAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGcccc                                  >  1:2506909/1‑93 (MQ=255)
                                                        aGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCa                                 <  1:471348/93‑1 (MQ=255)
                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  1:891976/1‑93 (MQ=255)
                                                         gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                                >  1:6542592/1‑93 (MQ=255)
                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              <  1:8617629/93‑1 (MQ=255)
                                                           ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                              >  1:9474726/1‑93 (MQ=255)
                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             >  2:3048346/1‑93 (MQ=255)
                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             <  2:7731875/93‑1 (MQ=255)
                                                            cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                             >  2:4834483/1‑93 (MQ=255)
                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           >  2:9714774/1‑93 (MQ=255)
                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           <  1:79089/93‑1 (MQ=255)
                                                              gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACa                           >  1:6161279/1‑93 (MQ=255)
                                                                tgattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGt                         >  2:4222403/1‑93 (MQ=255)
                                                                 gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                        >  2:3344290/1‑93 (MQ=255)
                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       >  2:2627852/1‑93 (MQ=255)
                                                                  attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                       <  1:3752859/93‑1 (MQ=255)
                                                                     gatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGc                    >  2:4549465/1‑93 (MQ=255)
                                                                      atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                   >  2:5010265/1‑93 (MQ=255)
                                                                      atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                   >  1:3933488/1‑93 (MQ=255)
                                                                       tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                  <  1:1099454/93‑1 (MQ=255)
                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 <  2:5724801/93‑1 (MQ=255)
                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 >  2:2810538/1‑93 (MQ=255)
                                                                        ctctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAc                 <  1:3334560/93‑1 (MQ=255)
                                                                         tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACg                <  2:512795/93‑1 (MQ=255)
                                                                             ttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGGTTCGCCACGTCAt            >  2:160421/1‑93 (MQ=255)
                                                                               aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa          >  2:8849336/1‑93 (MQ=255)
                                                                               aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa          <  1:7656669/93‑1 (MQ=255)
                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         >  1:3471655/1‑93 (MQ=255)
                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         >  1:2435324/1‑93 (MQ=255)
                                                                                gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc         <  2:5767356/93‑1 (MQ=255)
                                                                                 ccGCTCCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGACCGCCCCAGTGACAGTTTCGCCACGTCATAACa        <  2:8830015/93‑1 (MQ=255)
                                                                                   gCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaaa      >  2:3415597/1‑93 (MQ=255)
                                                                                    cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAGGCCCCAGTGACAGTTTCGCCACGTTATAACAAAg     <  2:3564714/93‑1 (MQ=255)
                                                                                     ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt    <  1:1584348/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTC  >  NC_000913/1890135‑1890314

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
GCTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  >  NC_000913/1890132‑1890318
                                                                                             |                                                                                             
GCTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGT                                                                                               <  SRR4302540.7952157/94‑1 (MQ=60)
GCTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGT                                                                                               >  SRR4302540.2288546/1‑94 (MQ=60)
GCTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGT                                                                                               >  SRR4302540.1971090/1‑94 (MQ=60)
 CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTC                                                                                              >  SRR4302540.10033592/1‑94 (MQ=60)
 CTGAGCGGCATGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTC                                                                                              <  SRR4302540.1979788/94‑1 (MQ=60)
  TGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCG                                                                                             >  SRR4302540.6517604/1‑94 (MQ=60)
  TGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCG                                                                                             >  SRR4302540.5840155/1‑94 (MQ=60)
   GAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGT                                                                                            >  SRR4302540.381951/1‑94 (MQ=60)
    AGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTG                                                                                           <  SRR4302540.7637677/94‑1 (MQ=60)
      CGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTG                                                                                         >  SRR4302540.14017/1‑94 (MQ=60)
      CGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTG                                                                                         <  SRR4302540.2672764/94‑1 (MQ=60)
       GGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGC                                                                                        <  SRR4302540.2291563/94‑1 (MQ=60)
          GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     >  SRR4302540.6056151/1‑94 (MQ=60)
          GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     >  SRR4302540.6967695/1‑94 (MQ=60)
          GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     >  SRR4302540.9226203/1‑94 (MQ=60)
          GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     >  SRR4302540.9270846/1‑94 (MQ=60)
           GGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATA                                                                                    <  SRR4302540.2478341/94‑1 (MQ=60)
            GAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAC                                                                                   <  SRR4302540.7369790/94‑1 (MQ=60)
            GAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAC                                                                                   >  SRR4302540.4508649/1‑94 (MQ=60)
              GAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAA                                                                                 <  SRR4302540.815341/94‑1 (MQ=60)
               AGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAC                                                                                <  SRR4302540.4254874/94‑1 (MQ=60)
               AGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAC                                                                                <  SRR4302540.8550883/94‑1 (MQ=60)
                   AGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGG                                                                            >  SRR4302540.473667/1‑94 (MQ=60)
                    GAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGC                                                                           <  SRR4302540.9847516/94‑1 (MQ=60)
                    GAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGC                                                                           <  SRR4302540.6728654/94‑1 (MQ=60)
                      GTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGC                                                                         >  SRR4302540.9155293/1‑94 (MQ=60)
                          AATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTA                                                                     >  SRR4302540.2685706/1‑94 (MQ=60)
                           ATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCA                                                                    <  SRR4302540.9691524/94‑1 (MQ=60)
                            TGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAA                                                                   >  SRR4302540.4057400/1‑94 (MQ=60)
                              GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302540.3092482/94‑1 (MQ=60)
                              GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302540.724936/94‑1 (MQ=60)
                              GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302540.3820755/94‑1 (MQ=60)
                               CGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATG                                                                <  SRR4302540.4683057/94‑1 (MQ=60)
                                 CGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCC                                                              >  SRR4302540.9457549/1‑94 (MQ=60)
                                   AAGTTCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCGGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGC                                                            >  SRR4302540.859220/1‑94 (MQ=60)
                                     GTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCAT                                                          <  SRR4302540.6032738/94‑1 (MQ=60)
                                      TGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATC                                                         <  SRR4302540.7664719/94‑1 (MQ=60)
                                      TGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATC                                                         >  SRR4302540.5619319/1‑94 (MQ=60)
                                      TGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATC                                                         >  SRR4302540.9745191/1‑94 (MQ=60)
                                        CATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCA                                                       >  SRR4302540.9705155/1‑94 (MQ=60)
                                           TCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCA                                                    <  SRR4302540.1153540/94‑1 (MQ=60)
                                             GTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAG                                                  <  SRR4302540.5717625/94‑1 (MQ=60)
                                             GTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAG                                                  >  SRR4302540.3073681/1‑94 (MQ=60)
                                              TACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGG                                                 >  SRR4302540.5510328/1‑94 (MQ=60)
                                              TACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGG                                                 >  SRR4302540.9891068/1‑94 (MQ=60)
                                               ACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGC                                                >  SRR4302540.3980386/1‑94 (MQ=60)
                                               ACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGC                                                <  SRR4302540.8904663/94‑1 (MQ=60)
                                               ACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGCC                                                >  SRR4302540.9352896/1‑94 (MQ=60)
                                                CCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCG                                               >  SRR4302540.7522511/1‑94 (MQ=60)
                                                CCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCG                                               >  SRR4302540.4301358/1‑94 (MQ=60)
                                                CCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCG                                               <  SRR4302540.2351580/94‑1 (MQ=60)
                                                CCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCG                                               <  SRR4302540.3622769/94‑1 (MQ=60)
                                                 CAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGG                                              >  SRR4302540.1701326/1‑94 (MQ=60)
                                                   GTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAC                                            >  SRR4302540.9407588/1‑94 (MQ=60)
                                                   GTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAA                                            >  SRR4302540.8036961/1‑94 (MQ=60)
                                                    gatCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.5933775/91‑1 (MQ=60)
                                                    TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.1216384/94‑1 (MQ=60)
                                                    TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.8692575/94‑1 (MQ=60)
                                                    TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.5933776/94‑1 (MQ=60)
                                                    TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.5927970/94‑1 (MQ=60)
                                                    TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGCGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302540.9988765/94‑1 (MQ=60)
                                                       CTGTAGCCCGGCTGATTGATCTCTTTTAGACGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCC                                        >  SRR4302540.5597531/1‑94 (MQ=60)
                                                        TGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTTCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCC                                       >  SRR4302540.5443997/1‑94 (MQ=60)
                                                        TGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCC                                       >  SRR4302540.8139309/1‑94 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGCTACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      >  SRR4302540.195811/1‑94 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      <  SRR4302540.449874/94‑1 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      >  SRR4302540.3734501/1‑94 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      >  SRR4302540.2533980/1‑94 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      <  SRR4302540.449875/94‑1 (MQ=60)
                                                         GTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCC                                      <  SRR4302540.8400467/94‑1 (MQ=60)
                                                          TAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCA                                     <  SRR4302540.4389569/94‑1 (MQ=60)
                                                           AGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAG                                    <  SRR4302540.476400/94‑1 (MQ=60)
                                                           AGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAG                                    >  SRR4302540.7393273/1‑94 (MQ=60)
                                                           AGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAG                                    >  SRR4302540.901454/1‑94 (MQ=60)
                                                           AGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAG                                    >  SRR4302540.6615192/1‑94 (MQ=60)
                                                             CCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTG                                  >  SRR4302540.9587096/1‑94 (MQ=60)
                                                             CCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCAATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTG                                  >  SRR4302540.9587097/1‑94 (MQ=60)
                                                              CCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGA                                 <  SRR4302540.8717742/94‑1 (MQ=60)
                                                               CGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGAC                                >  SRR4302540.8146197/1‑94 (MQ=60)
                                                               CGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGAC                                <  SRR4302540.3603940/94‑1 (MQ=60)
                                                                GGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACA                               >  SRR4302540.6229933/1‑94 (MQ=60)
                                                                 GCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAG                              <  SRR4302540.79917/94‑1 (MQ=60)
                                                                   TGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTT                            >  SRR4302540.7338211/1‑94 (MQ=60)
                                                                    GATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTT                           >  SRR4302540.7962389/1‑94 (MQ=60)
                                                                    GATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTT                           <  SRR4302540.8005182/94‑1 (MQ=60)
                                                                     ATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTC                          <  SRR4302540.3793802/94‑1 (MQ=60)
                                                                     ATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTC                          <  SRR4302540.7061277/94‑1 (MQ=60)
                                                                       TGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGC                        <  SRR4302540.9967943/94‑1 (MQ=60)
                                                                        GATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCC                       >  SRR4302540.322024/1‑94 (MQ=60)
                                                                        GATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCC                       >  SRR4302540.3976007/1‑94 (MQ=60)
                                                                          TCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAC                     <  SRR4302540.1110776/94‑1 (MQ=60)
                                                                          TCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAC                     <  SRR4302540.540459/94‑1 (MQ=60)
                                                                           CTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACG                    <  SRR4302540.3371178/94‑1 (MQ=60)
                                                                             CTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTC                  >  SRR4302540.3913402/1‑94 (MQ=60)
                                                                                 TAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAA              >  SRR4302540.1699057/1‑94 (MQ=60)
                                                                                  AGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAC             >  SRR4302540.3509454/1‑94 (MQ=60)
                                                                                  AGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAC             >  SRR4302540.8091828/1‑94 (MQ=60)
                                                                                  AGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAC             >  SRR4302540.2461557/1‑94 (MQ=60)
                                                                                  AGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAC             <  SRR4302540.7742876/94‑1 (MQ=60)
                                                                                   GCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAACGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACA            >  SRR4302540.1767920/1‑94 (MQ=60)
                                                                                    CCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAA           >  SRR4302540.2167270/1‑94 (MQ=60)
                                                                                      GCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAG         <  SRR4302540.5447737/94‑1 (MQ=60)
                                                                                       CTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGT        >  SRR4302540.7546824/1‑94 (MQ=60)
                                                                                        TTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTT       <  SRR4302540.1601652/94‑1 (MQ=60)
                                                                                        TTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTC       >  SRR4302540.3390814/1‑94 (MQ=60)
                                                                                        TTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCACGTCATAACAAAGTTC      >  SRR4302540.9552062/1‑94 (MQ=60)
                                                                                           AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACACAGTTCCA    >  SRR4302540.1784755/1‑94 (MQ=60)
                                                                                           AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    >  SRR4302540.288393/1‑94 (MQ=60)
                                                                                           AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    >  SRR4302540.3599931/1‑94 (MQ=60)
                                                                                           AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    <  SRR4302540.7456959/94‑1 (MQ=60)
                                                                                           AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    >  SRR4302540.7872797/1‑94 (MQ=60)
                                                                                            GTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAG   >  SRR4302540.3007070/1‑94 (MQ=60)
                                                                                             TCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  <  SRR4302540.5841751/94‑1 (MQ=60)
                                                                                             |                                                                                             
GCTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  >  NC_000913/1890132‑1890318

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