Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,259,919 C→T R156H (CGC→CAC)  fruK ← fructose‑1‑phosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,9190CT99.4% 551.1 / NA 158R156H (CGC→CAC) fruKfructose‑1‑phosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (81/76);  minor base A (1/0);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

GGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTC  >  NC_000913/2259829‑2260008
                                                                                          |                                                                                         
ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGa                                                                                         >  2:119162/1‑93 (MQ=255)
ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGa                                                                                         >  1:3332333/1‑93 (MQ=255)
ggTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGa                                                                                         <  2:361062/93‑1 (MQ=255)
    tCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCa                                                                                     <  2:444122/93‑1 (MQ=255)
      aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   >  2:4538337/1‑93 (MQ=255)
      aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   >  2:1673679/1‑93 (MQ=255)
      aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   <  1:2672945/93‑1 (MQ=255)
      aCCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATc                                                                                   <  1:793065/93‑1 (MQ=255)
          gCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGt                                                                               >  1:2334544/1‑93 (MQ=255)
          gCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGt                                                                               <  1:2879025/93‑1 (MQ=255)
          gCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGt                                                                               <  2:874029/93‑1 (MQ=255)
           ccACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTc                                                                              <  2:4499904/93‑1 (MQ=255)
            cACGGAGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGCACACTGACTACGCAGGAGAGTCATCCAGTCg                                                                             >  1:3439932/1‑93 (MQ=255)
             aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGATTCATCCAGTCgg                                                                            >  2:2558193/1‑93 (MQ=255)
             aCGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCgg                                                                            <  1:4612657/93‑1 (MQ=255)
                gTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGa                                                                         >  1:4610402/1‑93 (MQ=255)
                 tGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGaa                                                                        <  2:3649444/93‑1 (MQ=255)
                 tGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGaa                                                                        >  1:1345357/1‑93 (MQ=255)
                  gCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAAc                                                                       <  2:437979/93‑1 (MQ=255)
                  gCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAAc                                                                       <  2:1959742/93‑1 (MQ=255)
                   ccGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACg                                                                      >  1:3997573/1‑93 (MQ=255)
                    cGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGc                                                                     >  1:4133883/1‑93 (MQ=255)
                    cGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGc                                                                     >  2:3711863/1‑93 (MQ=255)
                     gCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCt                                                                    <  2:3037820/93‑1 (MQ=255)
                     gCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCt                                                                    >  2:4362433/1‑93 (MQ=255)
                      cTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCtt                                                                   >  1:2027473/1‑93 (MQ=255)
                      cTTTCAAACCTGCTACTAACGCTTCACCGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATACAGTCGGTGAACGCtt                                                                   >  1:867715/1‑93 (MQ=255)
                       tttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTc                                                                  <  2:3681947/93‑1 (MQ=255)
                       tttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTc                                                                  >  1:3594361/1‑93 (MQ=255)
                        ttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTcc                                                                 >  2:2822181/1‑93 (MQ=255)
                        ttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTcc                                                                 >  2:969920/1‑93 (MQ=255)
                        ttCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTcc                                                                 >  1:2984053/1‑93 (MQ=255)
                         tCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCg                                                                <  2:4006670/93‑1 (MQ=255)
                          cAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCgg                                                               <  1:2334662/93‑1 (MQ=255)
                          cAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCgg                                                               <  2:2254558/93‑1 (MQ=255)
                               cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          >  2:4277431/1‑93 (MQ=255)
                               cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          >  1:2704120/1‑93 (MQ=255)
                               cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          <  2:190775/93‑1 (MQ=255)
                               cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          <  1:2466842/93‑1 (MQ=255)
                               cTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGa                                                          >  2:4418340/1‑93 (MQ=255)
                                tGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGAc                                                         <  2:2480923/93‑1 (MQ=255)
                                     ctaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTg                                                    <  1:1406043/93‑1 (MQ=255)
                                     ctaACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTACGTGAACGCTTCCGGGCTGACGACTg                                                    >  2:3792931/1‑93 (MQ=255)
                                      taACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGa                                                   <  2:4384121/93‑1 (MQ=255)
                                      taACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCACGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGa                                                   <  1:2509674/93‑1 (MQ=255)
                                       aaCGCTTCACGGCTGCTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAc                                                  <  1:2676181/93‑1 (MQ=255)
                                       aaCGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAc                                                  <  2:4642905/93‑1 (MQ=255)
                                       aaCGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAc                                                  <  1:4085602/93‑1 (MQ=255)
                                         cGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACgg                                                <  1:312546/93‑1 (MQ=255)
                                            ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTaa                                             >  2:631089/1‑93 (MQ=255)
                                            ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTaa                                             >  2:3203493/1‑93 (MQ=255)
                                            ttCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTaa                                             >  2:2255661/1‑93 (MQ=255)
                                             tCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAg                                            >  1:1080517/1‑93 (MQ=255)
                                               aCGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCt                                          >  2:1207993/1‑93 (MQ=255)
                                                cGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCtt                                         >  2:3380436/1‑93 (MQ=255)
                                                 ggCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTc                                        <  2:4457315/93‑1 (MQ=255)
                                                  gCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTcc                                       <  2:2195087/93‑1 (MQ=255)
                                                    tactaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGc                                     >  1:818295/1‑93 (MQ=255)
                                                     actaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCt                                    >  2:866503/1‑93 (MQ=255)
                                                     actaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCt                                    <  2:1128636/93‑1 (MQ=255)
                                                      ctaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTg                                   >  2:279248/1‑93 (MQ=255)
                                                      ctaTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTg                                   <  1:3520909/93‑1 (MQ=255)
                                                       taTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGa                                  <  1:1593089/93‑1 (MQ=255)
                                                        aTCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGac                                 >  1:241769/1‑93 (MQ=255)
                                                           aaaGATAATGCAAGCACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGAcaca                              <  2:757231/93‑1 (MQ=255)
                                                            aaGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAg                             <  1:2965973/93‑1 (MQ=255)
                                                              gATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAc                           >  2:1664173/1‑93 (MQ=255)
                                                               aTAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAcc                          >  1:4429001/1‑93 (MQ=255)
                                                               aTAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAcc                          >  1:3776990/1‑93 (MQ=255)
                                                                tAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCa                         <  2:3612248/93‑1 (MQ=255)
                                                                  aTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCata                       <  1:2198829/93‑1 (MQ=255)
                                                                  aTGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCata                       <  2:203825/93‑1 (MQ=255)
                                                                   tGCAAGGACCCTGACTACGCAGGTGAGTCATCCAGTCGGCGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCatat                      >  1:4197974/1‑93 (MQ=255)
                                                                   tGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCatct                      <  2:2957268/93‑3 (MQ=255)
                                                                   tGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCatat                      >  1:2061010/1‑93 (MQ=255)
                                                                   tGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCatat                      >  2:3453011/1‑93 (MQ=255)
                                                                    gCAAGGACACTGACTACGCAGGTGATTCATCCAGTCGGTAAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACGCAGGCCATATc                     <  2:2004692/93‑1 (MQ=255)
                                                                    gCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATc                     >  2:4156349/1‑93 (MQ=255)
                                                                      aaGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGa                   >  2:1080762/1‑93 (MQ=255)
                                                                      aaGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGa                   >  2:4561815/1‑93 (MQ=255)
                                                                        ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 >  1:3050630/1‑93 (MQ=255)
                                                                        ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 >  2:3590393/1‑93 (MQ=255)
                                                                        ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 <  2:507898/93‑1 (MQ=255)
                                                                        ggACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAc                 >  1:1259396/1‑93 (MQ=255)
                                                                            acTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc             <  1:950630/93‑1 (MQ=255)
                                                                            acTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc             >  2:294337/1‑93 (MQ=255)
                                                                            acTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAc             >  1:4182031/1‑93 (MQ=255)
                                                                               gACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGa          <  2:2980269/93‑1 (MQ=255)
                                                                                 cTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGc        <  1:1051403/93‑1 (MQ=255)
                                                                                 cTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGc        <  2:1347669/93‑1 (MQ=255)
                                                                                 cTACGCAGGTGAGTCATACAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGc        >  2:1961467/1‑93 (MQ=255)
                                                                                  tACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGcc       <  2:3933539/93‑1 (MQ=255)
                                                                                  tACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGcc       >  2:730275/1‑93 (MQ=255)
                                                                                   aCGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCa      <  1:1829796/93‑1 (MQ=255)
                                                                                   aCGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCa      >  2:466268/1‑93 (MQ=255)
                                                                                    cGCAGGTGAGTCATCCAGTCGGTGAACGCCTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAg     <  2:2358707/93‑1 (MQ=255)
                                                                                     gCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGc    <  2:3043451/93‑1 (MQ=255)
                                                                                     gCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGc    <  2:2620203/93‑1 (MQ=255)
                                                                                     gCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGc    >  1:4080685/1‑93 (MQ=255)
                                                                                       aGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTc  <  1:1391870/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTC  >  NC_000913/2259829‑2260008

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
AGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAG  >  NC_000913/2259828‑2260012
                                                                                           |                                                                                             
AGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGA                                                                                             >  SRR4302535.3355524/1‑94 (MQ=60)
       ACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCC                                                                                      <  SRR4302535.798638/94‑1 (MQ=60)
       ACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCC                                                                                      <  SRR4302535.2692068/94‑1 (MQ=60)
         CAGCCACGGTGCCGGTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAG                                                                                    <  SRR4302535.1411563/94‑1 (MQ=60)
          AGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGT                                                                                   >  SRR4302535.2351250/1‑94 (MQ=60)
          AGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGT                                                                                   >  SRR4302535.1110703/1‑94 (MQ=60)
           GCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTC                                                                                  <  SRR4302535.2899410/94‑1 (MQ=60)
            CCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCG                                                                                 <  SRR4302535.2720569/94‑1 (MQ=60)
            CCACGGAGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGCACACTGACTACGCAGGAGAGTCATCCAGTCG                                                                                 >  SRR4302535.3463852/1‑94 (MQ=60)
              ACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGT                                                                               <  SRR4302535.4643995/94‑1 (MQ=60)
                GGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGA                                                                             >  SRR4302535.4641726/1‑94 (MQ=60)
                 GTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAA                                                                            >  SRR4302535.1354963/1‑94 (MQ=60)
                   GCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACG                                                                          >  SRR4302535.4024925/1‑94 (MQ=60)
                   GCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACG                                                                          <  SRR4302535.3408236/94‑1 (MQ=60)
                    CCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGC                                                                         >  SRR4302535.4162182/1‑94 (MQ=60)
                     CGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCCAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCT                                                                        >  SRR4302535.2134259/1‑94 (MQ=60)
                      GCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTT                                                                       >  SRR4302535.2042087/1‑94 (MQ=60)
                      GCTTTCAAACCTGCTACTAACGCTTCACCGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATACAGTCGGTGAACGCTT                                                                       >  SRR4302535.873857/1‑94 (MQ=60)
                       CTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTC                                                                      <  SRR4302535.2661418/94‑1 (MQ=60)
                       CTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTC                                                                      >  SRR4302535.3619305/1‑94 (MQ=60)
                        TTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCC                                                                     >  SRR4302535.3004984/1‑94 (MQ=60)
                        TTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCC                                                                     <  SRR4302535.3599636/94‑1 (MQ=60)
                        TTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCC                                                                     <  SRR4302535.4499690/94‑1 (MQ=60)
                           CAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGG                                                                  <  SRR4302535.2351371/94‑1 (MQ=60)
                               CCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGA                                                              >  SRR4302535.2723441/1‑94 (MQ=60)
                                CTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGAC                                                             <  SRR4302535.2484624/94‑1 (MQ=60)
                                CTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGAC                                                             >  SRR4302535.1378926/1‑94 (MQ=60)
                                    TACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCT                                                         <  SRR4302535.890550/94‑1 (MQ=60)
                                      CTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGA                                                       <  SRR4302535.1416091/94‑1 (MQ=60)
                                       TAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCACGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGAC                                                      <  SRR4302535.2527767/94‑1 (MQ=60)
                                        AACGCTTCACGGCTGCTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACG                                                     <  SRR4302535.2695328/94‑1 (MQ=60)
                                        AACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACG                                                     >  SRR4302535.237128/1‑94 (MQ=60)
                                        AACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACG                                                     <  SRR4302535.4113555/94‑1 (MQ=60)
                                          CGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGT                                                   <  SRR4302535.314710/94‑1 (MQ=60)
                                             TTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAG                                                >  SRR4302535.1088203/1‑94 (MQ=60)
                                                   GCTACTATCAAAGATAAGGCAAGGACCCTGACTCCGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGGTGACGCCTGACGGTAAGCTTCCG                                          <  SRR4302535.1264201/94‑1 (MQ=60)
                                                    CTACTATCAAGGATAATGCAAGGACACTGACTTCGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGC                                         <  SRR4302535.1548781/94‑1 (MQ=60)
                                                    CTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGC                                         >  SRR4302535.824037/1‑94 (MQ=60)
                                                      ACTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTG                                       >  SRR4302535.869588/1‑94 (MQ=60)
                                                       CTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGA                                      <  SRR4302535.3000349/94‑1 (MQ=60)
                                                       CTATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGA                                      <  SRR4302535.3545362/94‑1 (MQ=60)
                                                        TATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGAC                                     <  SRR4302535.1604586/94‑1 (MQ=60)
                                                        TATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGAC                                     >  SRR4302535.243427/1‑94 (MQ=60)
                                                         ATCAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACA                                    >  SRR4302535.4127033/1‑94 (MQ=60)
                                                           CAAAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACA                                  >  SRR4302535.2524969/1‑94 (MQ=60)
                                                             AAGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGA                                <  SRR4302535.2986783/94‑1 (MQ=60)
                                                              AGATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGAC                               >  SRR4302535.2852555/1‑94 (MQ=60)
                                                               GATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACC                              >  SRR4302535.3803029/1‑94 (MQ=60)
                                                               GATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACC                              >  SRR4302535.4459342/1‑94 (MQ=60)
                                                               AATAATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACC                              <  SRR4302535.3948378/94‑1 (MQ=60)
                                                                   ATGCAAGGACCCTGACTACGCAGGTGAGTCATCCAGTCGGCGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATAT                          >  SRR4302535.4226700/1‑94 (MQ=60)
                                                                   ATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATAT                          >  SRR4302535.2075862/1‑94 (MQ=60)
                                                                   ATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATAT                          <  SRR4302535.316729/94‑1 (MQ=60)
                                                                   ATGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATAT                          <  SRR4302535.2214495/94‑1 (MQ=60)
                                                                    TGCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATC                         >  SRR4302535.620175/1‑94 (MQ=60)
                                                                     GCAAGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCG                        <  SRR4302535.2267147/94‑1 (MQ=60)
                                                                        AGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAC                     >  SRR4302535.3071999/1‑94 (MQ=60)
                                                                        AGGACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAAC                     >  SRR4302535.1268409/1‑94 (MQ=60)
                                                                           ACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGA                  >  SRR4302535.4067376/1‑94 (MQ=60)
                                                                           ACACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCAGCTGACACAGACCATATCGAACTGA                  >  SRR4302535.2002372/1‑94 (MQ=60)
                                                                            CACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGAC                 >  SRR4302535.4210646/1‑94 (MQ=60)
                                                                             ACTGACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACC                <  SRR4302535.957361/94‑1 (MQ=60)
                                                                                GACTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAG             >  SRR4302535.4233578/1‑94 (MQ=60)
                                                                                  CTACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCC           <  SRR4302535.1058881/94‑1 (MQ=60)
                                                                                    ACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAG         <  SRR4302535.1843007/94‑1 (MQ=60)
                                                                                    ACGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAG         >  SRR4302535.4177368/1‑94 (MQ=60)
                                                                                     CGCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGC        >  SRR4302535.4108598/1‑94 (MQ=60)
                                                                                      GCAGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCT       <  SRR4302535.3024742/94‑1 (MQ=60)
                                                                                        AGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCA     <  SRR4302535.1638092/94‑1 (MQ=60)
                                                                                        AGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCA     <  SRR4302535.1401804/94‑1 (MQ=60)
                                                                                        AGGTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCA     >  SRR4302535.3131026/1‑94 (MQ=60)
                                                                                          GTGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGA   >  SRR4302535.4216204/1‑94 (MQ=60)
                                                                                           TGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAG  <  SRR4302535.617390/94‑1 (MQ=60)
                                                                                           |                                                                                             
AGGTTTCACCAGCCACGGTGCCGCTTTCAAACCTGCTACTAACGCTTCACGGCTACTATCAAAGATAATGCAAGGACACTGACTACGCAGGCGAGTCATCCAGTCGGTGAACGCTTCCGGGCTGACGCCTGACGGTAAGCTTCCGCTGACACAGACCATATCGAACTGACCGAGCCAGCTCAGAG  >  NC_000913/2259828‑2260012

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