BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,717	+G	coding (3/999 nt)	ysgA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,717	1	.	G	99.2%	314.0 / ‑5.7	120	M1M (ATG→ATC) †	ysgA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base G (56/62); minor base T (0/1); total (56/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 126 total aligned reads displayed.

CCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CCTTATAAAAATGT  >  NC_000913/4013642‑4013784

cccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                                       >  2:2820229/1‑76 (MQ=255)

cccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGACGAAGAGAGACGc                                                                       >  2:1393563/1‑76 (MQ=255)

 ccTTGTGAAGGGATCGAGGAGAAGCCCGCCACAAGCACAACGTCAGGGGTATGCACGATGGCCGAAGCGAGAGGcg                                                                      >  1:4752128/1‑75 (MQ=255)

  cTTGTGAAGGGATAGAGTTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                     <  1:4573423/76‑3 (MQ=255)

  cTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                     >  1:3363499/1‑74 (MQ=255)

      tgAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgca                                                                 >  1:2648108/1‑76 (MQ=255)

       gAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcag                                                                <  1:105899/76‑1 (MQ=255)

        aaGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                               >  1:152695/1‑76 (MQ=255)

        aaGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                               >  2:1406937/1‑76 (MQ=255)

        aaGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                               <  1:2858836/76‑1 (MQ=255)

          gggATAGAGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                                             >  1:1127695/1‑76 (MQ=255)

           ggATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAgg                                                            >  1:3544040/1‑76 (MQ=255)

           ggATAGAGGTGAAGCCCGCCACAATCGCGTCGCCCGGGGCCTGAACGATGGTCGAAGCGAGAGGCGCTGCAGCAgg                                                            >  2:2402888/1‑76 (MQ=255)

             aTAGGGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                          <  1:5376981/76‑1 (MQ=255)

              gagagGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCATCAGGTGc                                                         <  2:190384/75‑1 (MQ=255)

               agagGTGAGCCCCGCCACAATCGCGTCGTCCGGGGTTTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCa                                                        <  2:876947/76‑1 (MQ=255)

               agagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGCTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCa                                                        >  2:2672866/1‑76 (MQ=255)

                gagGTGACCCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                       <  2:3480228/76‑1 (MQ=255)

                agtgTGAAGCCCACCAAAATCGCGTCGTCCGGGGTCTGAACGAGGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                       <  2:177360/73‑1 (MQ=255)

                 agGTGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                      <  1:5824231/76‑1 (MQ=255)

                  ggTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAt                                                     >  1:3491319/1‑76 (MQ=255)

                     gAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                  <  2:3941674/76‑1 (MQ=255)

                     gAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                  <  1:3812481/76‑1 (MQ=255)

                     gAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                  >  2:3682187/1‑76 (MQ=255)

                       aGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCCGa                                                >  2:4185197/1‑76 (MQ=255)

                       aGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                >  1:4659134/1‑76 (MQ=255)

                       aGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                <  1:4867630/76‑1 (MQ=255)

                       aGCCCGCCACAATCGCGGCCTCCGGGGCCTGAACGATGGTCGAAGCGAGAGGCGATGCAGGAGGTGCTAATCCAGa                                                >  2:2994195/1‑76 (MQ=255)

                         cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                              <  1:1265267/76‑1 (MQ=255)

                         cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                              <  1:3212094/76‑1 (MQ=255)

                         cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                              <  1:1921631/76‑1 (MQ=255)

                          ccGCCACAAGCGCGGCGGCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                                             <  2:141664/76‑1 (MQ=255)

                            gCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgtt                                           >  1:2681652/1‑76 (MQ=255)

                                 aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                      >  2:3967523/1‑76 (MQ=255)

                                    cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                   <  2:5645630/76‑1 (MQ=255)

                                    cgcgTCGCCCGGGGTCTGAACGAAGGTCGAAGCGAGAGGCGACGCAGCAGGAGCAAACCCAGATtgctgtgttgct                                   >  2:2633627/1‑74 (MQ=255)

                                    cacgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCCGGTGCAAATCCAGATTGTTGTgttgtt                                   >  2:4844821/3‑76 (MQ=255)

                                     gcggcgGCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                  <  2:1254515/76‑1 (MQ=255)

                                     gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                  >  1:2741899/1‑76 (MQ=255)

                                      cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                 <  1:1844425/76‑1 (MQ=255)

                                      cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                 >  2:2090763/1‑76 (MQ=255)

                                      cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                 >  2:4037348/1‑76 (MQ=255)

                                      cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                 >  2:5338145/1‑76 (MQ=255)

                                      cgtcgtCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc                                 <  1:679854/76‑1 (MQ=255)

                                       gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                >  1:952986/1‑76 (MQ=255)

                                       gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                <  2:2043156/76‑1 (MQ=255)

                                       gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                <  1:3440955/76‑1 (MQ=255)

                                        tcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                               >  1:3688166/1‑76 (MQ=255)

                                         cgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCAt                              <  2:1129303/76‑1 (MQ=255)

                                          gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                             >  1:1434335/1‑76 (MQ=255)

                                           tcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATgg                            <  1:4231604/76‑1 (MQ=255)

                                            ccGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                           <  2:975071/76‑1 (MQ=255)

                                            ccGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                           <  1:492589/76‑1 (MQ=255)

                                            ccGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                           <  1:5458388/76‑1 (MQ=255)

                                             cGGGGTCTGAGCGATGGGCGAAGCGAGAGGCGATGCAGCCGGTGCCAATCCAGCGGGTTGTGTTGTTGCCATGGTa                          <  2:2195948/76‑1 (MQ=255)

                                             cGGGGTCTGACCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa                          <  2:2191263/76‑1 (MQ=255)

                                             cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa                          >  1:3294937/1‑76 (MQ=255)

                                              ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt                         <  1:2447292/76‑1 (MQ=255)

                                              ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt                         >  1:5057484/1‑76 (MQ=255)

                                              ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt                         >  1:5798757/1‑76 (MQ=255)

                                                ggTCTGAGCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                       <  2:399414/76‑1 (MQ=255)

                                                ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                       >  1:2021037/1‑76 (MQ=255)

                                                ggTCGGAACGATGGGCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                       <  2:1901821/76‑1 (MQ=255)

                                                   cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                    <  1:2141559/76‑1 (MQ=255)

                                                   cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                    >  2:3066372/1‑76 (MQ=255)

                                                   cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                    <  1:1552640/76‑1 (MQ=255)

                                                    tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                   <  2:1140967/76‑1 (MQ=255)

                                                    tGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                   >  2:4668546/1‑76 (MQ=255)

                                                     gAACTATGTTCGCAGCGAGAGGCAATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGt                  <  1:2994195/76‑1 (MQ=255)

                                                      aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                 <  1:979279/76‑1 (MQ=255)

                                                       ccGAGGGTGGAAGCGAGAGGCGGTGCAGCAGGTGCCCAGCCAGATTGTTGTTTTGTTGCCATGGTATTCTCCGTA‑c               <  1:831402/75‑1 (MQ=255)

                                                       aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCCTGGTATTCTCCGTA‑c               >  2:121054/1‑76 (MQ=255)

                                                       aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑c               <  2:4243448/76‑1 (MQ=255)

                                                        cGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑cc              <  1:2260185/76‑1 (MQ=255)

                                                         gATGGTCGAAGCGGGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc             <  1:5698577/76‑2 (MQ=255)

                                                         gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc             <  1:3721435/76‑2 (MQ=255)

                                                         gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc             <  2:4209006/76‑2 (MQ=255)

                                                         gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccc             >  1:5342439/1‑75 (MQ=255)

                                                          aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccct            <  1:260926/76‑3 (MQ=255)

                                                           tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑ccctt           <  2:769324/76‑4 (MQ=255)

                                                           tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATGGTTGTGTTGTTGCCATGGTATTGTCCGTA‑ccctt           <  2:791744/76‑4 (MQ=255)

                                                              tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata         >  1:36462/1‑76 (MQ=255)

                                                               ggAAGCGAGAGGCTATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTACCATGGTATTCTCCGTACCCTTATaa        <  1:2587659/75‑1 (MQ=255)

                                                               cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa        >  1:1297728/1‑76 (MQ=255)

                                                                gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa       >  2:3203476/1‑76 (MQ=255)

                                                                gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa       >  2:4195047/1‑76 (MQ=255)

                                                                gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa       >  1:2791574/1‑76 (MQ=255)

                                                                gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa       <  2:1481556/76‑1 (MQ=255)

                                                                gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa       >  1:3898726/1‑76 (MQ=255)

                                                                 aaGCGAGAGGCGATGCAGCAGGTGCAAATCCCGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa      <  1:1224137/76‑1 (MQ=255)

                                                                 aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa      <  2:299599/76‑1 (MQ=255)

                                                                 aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa      >  2:747607/1‑76 (MQ=255)

                                                                 aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGCACCCTTATaaaa      >  1:141423/1‑76 (MQ=255)

                                                                  aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa     <  1:802022/76‑1 (MQ=255)

                                                                  aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa     <  2:4719208/76‑1 (MQ=255)

                                                                   gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt    >  2:5063136/1‑76 (MQ=255)

                                                                   gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt    <  2:3491319/76‑1 (MQ=255)

                                                                    cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg   >  1:346773/1‑76 (MQ=255)

                                                                    cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg   >  2:1597425/1‑76 (MQ=255)

                                                                     gagaGGCGATGCAGCATGTTCAAATCCAGATTGTTGTGTTGTTGCTATGGTATTCTCCCTACCCTTATATAAATGt  <  1:1035332/76‑1 (MQ=255)

CCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGC‑ATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTA‑CCTTATAAAAATGT  >  NC_000913/4013642‑4013784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 8 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`