Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,896,324 +G coding (698/711 nt) b3718 ← putative isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,896,3241.G100.0% 871.0 / NA 241V233A (GTG→GCG) b3718putative isomerase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (134/107);  total (134/107)

Only 100 of 256 total aligned reads displayed.

TTCGTAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410
                                                                                              |                                                                                      
ttCGTAG‑TAATTAAATTAGGAGGACCGGGCCGGAGCGTATTTCACCAATCCGGCGTTA‑TTCCTAAGCTTAACCGAGTGCTAATTCTGCGGCag                                                                                        <  2:646915/93‑2 (MQ=255)
  cGTAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagg                                                                                       >  2:1308461/1‑91 (MQ=255)
  cGTAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagc                                                                                       <  2:2821280/93‑3 (MQ=255)
    tAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagcgg                                                                                     >  1:3329308/1‑89 (MQ=255)
    tAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCagcgg                                                                                     >  1:4666662/1‑89 (MQ=255)
      gTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCt                                                                                   >  2:5346095/1‑93 (MQ=255)
       ttATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCtt                                                                                  >  2:5917253/1‑93 (MQ=255)
       ttATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCtt                                                                                  >  1:3807712/1‑93 (MQ=255)
        tATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCttt                                                                                 >  1:750153/1‑93 (MQ=255)
         aTTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTa                                                                                >  2:5673140/1‑93 (MQ=255)
          tttAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTAt                                                                               >  2:4676475/1‑93 (MQ=255)
          tttAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTAt                                                                               >  1:1098701/1‑93 (MQ=255)
           ttAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATc                                                                              <  2:753084/93‑1 (MQ=255)
           ttAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATc                                                                              <  1:4006422/93‑1 (MQ=255)
            tAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCa                                                                             >  1:2278835/1‑93 (MQ=255)
            tAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCa                                                                             >  1:3466415/1‑93 (MQ=255)
             aaaTAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAg                                                                            <  2:3085155/93‑1 (MQ=255)
              aataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGc                                                                           >  1:599914/1‑93 (MQ=255)
              aataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGc                                                                           <  1:5001693/93‑1 (MQ=255)
              aataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGc                                                                           <  2:82705/93‑1 (MQ=255)
               ataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCg                                                                          >  1:3038595/1‑93 (MQ=255)
               ataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCg                                                                          >  2:1505271/1‑93 (MQ=255)
               ataatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGCCAGCGGCTTTATCAGCg                                                                          >  2:4998536/1‑93 (MQ=255)
                taatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGa                                                                         >  1:4852512/1‑93 (MQ=255)
                 aatGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTCTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAt                                                                        >  1:803075/1‑93 (MQ=255)
                   tGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTa                                                                      >  2:1586332/1‑93 (MQ=255)
                   tGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTa                                                                      >  1:5780937/1‑93 (MQ=255)
                     a‑gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAcc                                                                    >  2:5792290/1‑93 (MQ=255)
                       gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCa                                                                   >  1:3953601/1‑93 (MQ=255)
                       gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCa                                                                   >  1:2209845/1‑93 (MQ=255)
                       gaACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCa                                                                   >  1:5748510/1‑93 (MQ=255)
                        aaCAGGCCGGAGCGTAATTCACACCTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCAt                                                                  <  1:5475705/93‑1 (MQ=255)
                          cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCa                                                                <  1:2168400/93‑1 (MQ=255)
                          cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCa                                                                <  1:677780/93‑1 (MQ=255)
                          cAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCa                                                                >  1:242698/1‑93 (MQ=255)
                           aGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAg                                                               <  1:1473652/93‑1 (MQ=255)
                           aGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAg                                                               <  1:5366390/93‑1 (MQ=255)
                            ggCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGc                                                              >  1:3483823/1‑93 (MQ=255)
                            ggCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGc                                                              >  1:4921107/1‑93 (MQ=255)
                            ggCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGc                                                              >  1:668551/1‑93 (MQ=255)
                              ccGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGa                                                            <  2:5525544/93‑1 (MQ=255)
                               cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                           <  2:365180/93‑1 (MQ=255)
                               cGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGaa                                                           >  2:4846512/1‑93 (MQ=255)
                                ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  2:1674439/1‑93 (MQ=255)
                                ggAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAg                                                          >  1:4122574/1‑93 (MQ=255)
                                 gAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGTGATTACCATCAGCGAAgg                                                         >  1:1335542/1‑93 (MQ=255)
                                 gAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAgg                                                         <  1:1608136/93‑1 (MQ=255)
                                 gAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAgg                                                         >  2:5067068/1‑93 (MQ=255)
                                   gCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAt                                                       <  2:4192291/93‑1 (MQ=255)
                                   gCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAt                                                       >  1:2900933/1‑93 (MQ=255)
                                    cGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATg                                                      <  1:4431309/93‑1 (MQ=255)
                                      tAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCa                                                    >  2:4451008/1‑93 (MQ=255)
                                       aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   >  1:100415/1‑93 (MQ=255)
                                       aaTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAg                                                   <  2:4429651/93‑1 (MQ=255)
                                        aTTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGc                                                  <  2:3338907/93‑1 (MQ=255)
                                          tCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTg                                                <  2:2233183/93‑1 (MQ=255)
                                           cacacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               <  1:5150310/93‑1 (MQ=255)
                                           cacacaCCCGGCCTTATTGCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGc                                               <  2:1279391/93‑1 (MQ=255)
                                            acacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCa                                              <  1:5250667/93‑1 (MQ=255)
                                             cacaTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCaa                                             <  1:570962/93‑1 (MQ=255)
                                               caTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                           >  2:1412631/1‑93 (MQ=255)
                                               caTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca                                           >  1:699233/1‑93 (MQ=255)
                                                aTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAacac                                          >  1:938260/1‑93 (MQ=255)
                                                   cGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGa                                       >  2:4554954/1‑93 (MQ=255)
                                                      ccTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGc                                    >  2:2374695/1‑93 (MQ=255)
                                                       cTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCt                                   <  2:6042958/93‑1 (MQ=255)
                                                         tATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                 >  2:1897274/1‑93 (MQ=255)
                                                         tATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTgg                                 >  2:5258374/1‑93 (MQ=255)
                                                          aTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGa                                <  1:6100460/93‑1 (MQ=255)
                                                           tttCTTAAGCTTAACCGAGTGTTAATTCCGCGGCAGCGGCTTTATCAGCTATTACCATCAGCGAAGGATGCAGCTGCAAAACGGAAGCTGGaa                               >  1:2265086/1‑93 (MQ=255)
                                                             tCTTGAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACa                             <  1:2127493/93‑1 (MQ=255)
                                                              cTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACAt                            <  2:2613219/93‑1 (MQ=255)
                                                                 aaGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCtt                         <  2:1524031/93‑1 (MQ=255)
                                                                 aaGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCtt                         <  1:968494/93‑1 (MQ=255)
                                                                  aGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  1:4207066/93‑1 (MQ=255)
                                                                  aGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTc                        <  1:271212/93‑1 (MQ=255)
                                                                   gCTTAACCGAGTGCTAATTCTGCGGCAGCTGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       <  2:3007988/93‑1 (MQ=255)
                                                                   gCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCg                       >  1:1945966/1‑93 (MQ=255)
                                                                    cTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCgg                      <  1:4311387/93‑1 (MQ=255)
                                                                     ttAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGCACATCTTCGGt                     >  1:4246938/1‑93 (MQ=255)
                                                                         ccGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcc                 <  1:5398534/93‑1 (MQ=255)
                                                                         ccGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAcc                 <  1:4412589/93‑1 (MQ=255)
                                                                          cGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCg                >  2:3800492/1‑93 (MQ=255)
                                                                           gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  1:2536657/1‑93 (MQ=255)
                                                                           gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  1:3394198/1‑93 (MQ=255)
                                                                           gAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCgg               >  1:1730207/1‑93 (MQ=255)
                                                                            aGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGa              <  2:4651139/93‑1 (MQ=255)
                                                                            aGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGa              >  2:4106082/1‑93 (MQ=255)
                                                                              tGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAcc            >  2:2101231/1‑93 (MQ=255)
                                                                                ccaaTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  2:6081543/91‑1 (MQ=255)
                                                                                cTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCtt          <  1:1412491/93‑1 (MQ=255)
                                                                                 tAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTg         >  2:4359574/1‑93 (MQ=255)
                                                                                 tAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTg         <  1:5587498/93‑1 (MQ=255)
                                                                                  aaTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGc        >  2:1961520/1‑93 (MQ=255)
                                                                                   aTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCa       >  2:4263167/1‑93 (MQ=255)
                                                                                   aTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCa       <  2:1991142/93‑1 (MQ=255)
                                                                                   aTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCa       <  1:1241202/93‑1 (MQ=255)
                                                                                    ttCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAg      >  1:4586767/1‑93 (MQ=255)
                                                                                       tGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAAc   >  1:1882824/1‑93 (MQ=255)
                                                                                        gCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACa  <  2:1169599/93‑1 (MQ=255)
                                                                                              |                                                                                      
TTCGTAGTTATTTAAATAATGA‑GAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACA  >  NC_000913/3896232‑3896410

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCA  >  NC_000913/3896231‑3896417
                                                                                             |                                                                                              
TTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACTCATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA                                                                                                <  SRR4302532.3894215/94‑1 (MQ=60)
TTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA                                                                                                <  SRR4302532.1634607/94‑1 (MQ=60)
TTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA                                                                                                <  SRR4302532.773741/94‑1 (MQ=60)
 TTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑G                                                                                              >  SRR4302532.115418/1‑94 (MQ=60)
   CGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAgcg                                                                                             >  SRR4302532.164701/1‑91 (MQ=60)
    GTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGG                                                                                            >  SRR4302532.3440323/1‑94 (MQ=60)
    GTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGG                                                                                            >  SRR4302532.4825768/1‑94 (MQ=60)
    GTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGG                                                                                            >  SRR4302532.4203470/1‑94 (MQ=60)
      AGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCT                                                                                          <  SRR4302532.2500340/94‑1 (MQ=60)
       GTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTT                                                                                         >  SRR4302532.3936263/1‑94 (MQ=60)
        TTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTT                                                                                        <  SRR4302532.1953549/94‑1 (MQ=60)
        TTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTT                                                                                        <  SRR4302532.5875433/94‑1 (MQ=60)
        TTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTT                                                                                        >  SRR4302532.775844/1‑94 (MQ=60)
         TATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTA                                                                                       >  SRR4302532.722239/1‑94 (MQ=60)
          ATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTAT                                                                                      >  SRR4302532.1135429/1‑94 (MQ=60)
           TTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATC                                                                                     >  SRR4302532.2253685/1‑94 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    >  SRR4302532.2354467/1‑94 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    >  SRR4302532.3582512/1‑94 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    <  SRR4302532.638934/94‑1 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    <  SRR4302532.1721652/94‑1 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    <  SRR4302532.638935/94‑1 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    <  SRR4302532.2607649/94‑1 (MQ=60)
            TTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCA                                                                                    <  SRR4302532.4141422/94‑1 (MQ=60)
             TAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAG                                                                                   <  SRR4302532.4397832/94‑1 (MQ=60)
             TAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAG                                                                                   <  SRR4302532.981549/94‑1 (MQ=60)
              AAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGC                                                                                  >  SRR4302532.620859/1‑94 (MQ=60)
               AATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCG                                                                                 >  SRR4302532.3138854/1‑94 (MQ=60)
               AATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCG                                                                                 >  SRR4302532.3372100/1‑94 (MQ=60)
               AATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCG                                                                                 <  SRR4302532.5173146/94‑1 (MQ=60)
                ATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGA                                                                                >  SRR4302532.5019490/1‑94 (MQ=60)
                 TAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTCTTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAT                                                                               >  SRR4302532.830503/1‑94 (MQ=60)
                 TAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGAT                                                                               >  SRR4302532.3994872/1‑94 (MQ=60)
                  AATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATT                                                                              >  SRR4302532.265423/1‑94 (MQ=60)
                  AATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATT                                                                              >  SRR4302532.2378516/1‑94 (MQ=60)
                   ATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTA                                                                             >  SRR4302532.473162/1‑94 (MQ=60)
                   ATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTA                                                                             >  SRR4302532.1903704/1‑94 (MQ=60)
                   ATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTA                                                                             >  SRR4302532.5982013/1‑94 (MQ=60)
                   ATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCAATTA                                                                             >  SRR4302532.567155/1‑94 (MQ=60)
                   ATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGGGCTAATTCTGCGGCAGCCGCTTTATCAACGATTA                                                                             >  SRR4302532.5926765/1‑94 (MQ=60)
                    TGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAC                                                                            >  SRR4302532.577879/1‑94 (MQ=60)
                      AGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCA                                                                          >  SRR4302532.1121461/1‑94 (MQ=60)
                      AGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCA                                                                          >  SRR4302532.5948608/1‑94 (MQ=60)
                      AGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCA                                                                          >  SRR4302532.2283290/1‑94 (MQ=60)
                      AGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCA                                                                          >  SRR4302532.4086833/1‑94 (MQ=60)
                        AACAGGCCGGAGCGTAATTCACACCTCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATC                                                                        <  SRR4302532.5665565/94‑1 (MQ=60)
                        AACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATC                                                                        >  SRR4302532.5092545/1‑94 (MQ=60)
                         ACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCA                                                                       >  SRR4302532.251086/1‑94 (MQ=60)
                         ACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCA                                                                       <  SRR4302532.5724330/94‑1 (MQ=60)
                         ACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCA                                                                       <  SRR4302532.1428347/94‑1 (MQ=60)
                          CAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAG                                                                      >  SRR4302532.960832/1‑94 (MQ=60)
                          CAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAG                                                                      <  SRR4302532.2240913/94‑1 (MQ=60)
                          CAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAG                                                                      <  SRR4302532.701086/94‑1 (MQ=60)
                          CAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAG                                                                      <  SRR4302532.701088/94‑1 (MQ=60)
                           AGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGC                                                                     >  SRR4302532.3600523/1‑94 (MQ=60)
                           AGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGC                                                                     >  SRR4302532.691618/1‑94 (MQ=60)
                           AGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGC                                                                     <  SRR4302532.1523302/94‑1 (MQ=60)
                           AGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGC                                                                     >  SRR4302532.5090146/1‑94 (MQ=60)
                           AGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGC                                                                     <  SRR4302532.5552392/94‑1 (MQ=60)
                            GGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCG                                                                    >  SRR4302532.5085940/1‑94 (MQ=60)
                              CCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAA                                                                  >  SRR4302532.5980929/1‑94 (MQ=60)
                               CGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAG                                                                 >  SRR4302532.4260935/1‑94 (MQ=60)
                               CGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAG                                                                 <  SRR4302532.3538240/94‑1 (MQ=60)
                                GGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGTGATTACCATCAGCGAAGG                                                                >  SRR4302532.1380215/1‑94 (MQ=60)
                                 GAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGA                                                               >  SRR4302532.695403/1‑94 (MQ=60)
                                 GAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGA                                                               <  SRR4302532.1662576/94‑1 (MQ=60)
                                  AGCGTAATTCACAGATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAT                                                              <  SRR4302532.564223/94‑1 (MQ=60)
                                  AGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAT                                                              <  SRR4302532.5419862/94‑1 (MQ=60)
                                  AGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGAT                                                              >  SRR4302532.2997732/1‑94 (MQ=60)
                                  AGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAAGAT                                                              >  SRR4302532.147260/1‑94 (MQ=60)
                                   GCGTAATTCACACATCCGTCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATG                                                             <  SRR4302532.3526792/94‑1 (MQ=60)
                                   GCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATG                                                             <  SRR4302532.5462807/94‑1 (MQ=60)
                                    CGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGC                                                            >  SRR4302532.4706191/1‑94 (MQ=60)
                                    CGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGC                                                            <  SRR4302532.4582295/94‑1 (MQ=60)
                                     GTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCA                                                           >  SRR4302532.4071936/1‑94 (MQ=60)
                                     GTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCA                                                           <  SRR4302532.4820857/94‑1 (MQ=60)
                                      TAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAG                                                          >  SRR4302532.103902/1‑94 (MQ=60)
                                        ATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT                                                        >  SRR4302532.3256469/1‑94 (MQ=60)
                                        ATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT                                                        <  SRR4302532.933874/94‑1 (MQ=60)
                                        ATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT                                                        <  SRR4302532.1948444/94‑1 (MQ=60)
                                        ATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCT                                                        <  SRR4302532.2980735/94‑1 (MQ=60)
                                         TTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTG                                                       <  SRR4302532.6291246/94‑1 (MQ=60)
                                           CACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCA                                                     <  SRR4302532.5327489/94‑1 (MQ=60)
                                           CACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCA                                                     <  SRR4302532.5524159/94‑1 (MQ=60)
                                            ACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAA                                                    <  SRR4302532.5431674/94‑1 (MQ=60)
                                             CACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAAC                                                   <  SRR4302532.590883/94‑1 (MQ=60)
                                              ACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACC                                                  >  SRR4302532.4808340/1‑94 (MQ=60)
                                              ACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACA                                                  >  SRR4302532.723159/1‑94 (MQ=60)
                                              ACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACA                                                  >  SRR4302532.4267494/1‑94 (MQ=60)
                                               CATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACAC                                                 >  SRR4302532.970416/1‑94 (MQ=60)
                                                ATCCGGCCTTATTTCTTAAGCTTAACCGGGCGCTAATTCTGCGGCAGCGGCTTTATCAGCGCTTACCATCAGCGAAGGATGCAGCTGCAACACG                                                >  SRR4302532.1376174/1‑94 (MQ=60)
                                                ATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCCACACG                                                >  SRR4302532.785971/1‑94 (MQ=60)
                                                ATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG                                                <  SRR4302532.2490554/94‑1 (MQ=60)
                                                  CCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGA                                              <  SRR4302532.2750033/94‑1 (MQ=60)
                                                  CCGGCCTTATTTCTTAAGCTCAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGA                                              >  SRR4302532.6048502/1‑94 (MQ=60)
                                                   CGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAA                                             <  SRR4302532.4136617/94‑1 (MQ=60)
                                                    GGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAG                                            <  SRR4302532.285631/94‑1 (MQ=60)
                                                      CCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCT                                          <  SRR4302532.187399/94‑1 (MQ=60)
                                                      CCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCT                                          >  SRR4302532.5923220/1‑94 (MQ=60)
                                                      CCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCT                                          >  SRR4302532.2154713/1‑94 (MQ=60)
                                                       CTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTG                                         >  SRR4302532.4494390/1‑94 (MQ=60)
                                                       CTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTG                                         >  SRR4302532.4534036/1‑94 (MQ=60)
                                                       CTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTG                                         >  SRR4302532.4369242/1‑94 (MQ=60)
                                                       CTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTG                                         >  SRR4302532.4369244/1‑94 (MQ=60)
                                                        TTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGG                                        <  SRR4302532.3787676/94‑1 (MQ=60)
                                                        TTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGG                                        <  SRR4302532.1760049/94‑1 (MQ=60)
                                                          ATTTCTTAAGCTTAACCGAGTGTTAATTCCGCGGCAGCGGCTTTATCAGCTATTACCATCAGCGAAGGATGCAGCTGCAAAACGGAAGCTGGAA                                      >  SRR4302532.2340259/1‑94 (MQ=60)
                                                          ATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAA                                      >  SRR4302532.2340258/1‑94 (MQ=60)
                                                          ATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAA                                      <  SRR4302532.6311005/94‑1 (MQ=60)
                                                           TTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAAC                                     <  SRR4302532.3500468/94‑1 (MQ=60)
                                                           TTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAAC                                     <  SRR4302532.2351976/94‑1 (MQ=60)
                                                             TCTTGAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACAT                                   <  SRR4302532.2198575/94‑1 (MQ=60)
                                                               TTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCT                                 <  SRR4302532.3595309/94‑1 (MQ=60)
                                                                TAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTT                                >  SRR4302532.1620789/1‑94 (MQ=60)
                                                                 AAGCTTAACCTAGAACTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCTAAGGAAGCAGGTGCACCACGGAAGCTGGAACATCTTC                               >  SRR4302532.4886306/1‑94 (MQ=57)
                                                                 AAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTC                               <  SRR4302532.1001541/94‑1 (MQ=60)
                                                                  AGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCG                              <  SRR4302532.4348283/94‑1 (MQ=60)
                                                                  AGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCG                              <  SRR4302532.747301/94‑1 (MQ=60)
                                                                  AGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCG                              <  SRR4302532.280605/94‑1 (MQ=60)
                                                                  AGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCG                              >  SRR4302532.2010667/1‑94 (MQ=60)
                                                                   GCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGG                             >  SRR4302532.4162543/1‑94 (MQ=60)
                                                                    CTTAACCGAGTGCTATTTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGT                            <  SRR4302532.3636261/94‑1 (MQ=60)
                                                                    CTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGCACATCTTCGGT                            >  SRR4302532.4389631/1‑94 (MQ=60)
                                                                    CTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGT                            <  SRR4302532.4457149/94‑1 (MQ=60)
                                                                       AACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAAC                         >  SRR4302532.5284000/1‑94 (MQ=60)
                                                                        ACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTCCCATCAGCGAAGGATGCAGCTGCAACCCGGTAGCTGGAACATCTTCGGTAATC                        >  SRR4302532.773880/1‑94 (MQ=60)
                                                                         CCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCG                       <  SRR4302532.5585594/94‑1 (MQ=60)
                                                                         CCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCG                       <  SRR4302532.4562903/94‑1 (MQ=60)
                                                                         CCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCG                       >  SRR4302532.5658984/1‑94 (MQ=60)
                                                                          CGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGG                      >  SRR4302532.5086504/1‑94 (MQ=60)
                                                                          CGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGG                      >  SRR4302532.2621284/1‑94 (MQ=60)
                                                                          CGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGG                      >  SRR4302532.3007878/1‑94 (MQ=60)
                                                                          CGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGG                      >  SRR4302532.1788445/1‑94 (MQ=60)
                                                                          CGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGG                      >  SRR4302532.3507748/1‑94 (MQ=60)
                                                                           GAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGA                     >  SRR4302532.6160063/1‑94 (MQ=60)
                                                                            AGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAC                    >  SRR4302532.5559972/1‑94 (MQ=60)
                                                                            AGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCAGCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGAC                    <  SRR4302532.5188582/94‑1 (MQ=60)
                                                                             GTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACC                   >  SRR4302532.1172365/1‑94 (MQ=60)
                                                                                CTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTG                <  SRR4302532.2233621/94‑1 (MQ=60)
                                                                                CTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTG                <  SRR4302532.2336341/94‑1 (MQ=60)
                                                                                CTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTG                <  SRR4302532.1459858/94‑1 (MQ=60)
                                                                                CTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTG                <  SRR4302532.2233620/94‑1 (MQ=60)
                                                                                CTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTG                <  SRR4302532.1459857/94‑1 (MQ=60)
                                                                                 TAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGC               <  SRR4302532.5781466/94‑1 (MQ=60)
                                                                                 TAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGC               <  SRR4302532.4983949/94‑1 (MQ=60)
                                                                                   ATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAG             <  SRR4302532.1282727/94‑1 (MQ=60)
                                                                                   ATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAG             >  SRR4302532.4743159/1‑94 (MQ=60)
                                                                                      CTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAAC          >  SRR4302532.1945251/1‑94 (MQ=60)
                                                                                        GCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACAT        <  SRR4302532.4569692/94‑1 (MQ=60)
                                                                                         CGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATT       <  SRR4302532.3603835/94‑1 (MQ=60)
                                                                                          GGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTT      >  SRR4302532.5236176/1‑94 (MQ=60)
                                                                                          GGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTT      <  SRR4302532.714680/94‑1 (MQ=60)
                                                                                           GCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTT     <  SRR4302532.2645933/94‑1 (MQ=60)
                                                                                             G‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTTA  >  SRR4302532.5808141/1‑94 (MQ=60)
                                                                                             |                                                                                              
TTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCA‑CGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACGGAAGCTGGAACATCTTCGGTAACCGGACCTTGCAGAACATTTTTCA  >  NC_000913/3896231‑3896417

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